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ABSTRACT 



This dissertation examines the phenomenology of statistical hadronization at ul- 
trarelativistic energies. We start with an overview of current experimental and 
theoretical issues in Relativistic heavy ion physics. We then introduce statistical 
hadronization, and show how it gives a description of particle abundances and spec- 
tra through relativistic covariance and entropy maximization. 

We argue that several statistical hadronization models are possible; In par- 
ticular, a distinction can be made between equilibrated staged freeze-out in which 
post-formation hadron interactions play an important role in determining final-state 
observables, and non-equilibrium sudden freeze-out where spectra and abundances 
get determined at the same time and further interactions are negligible. 

We attempt to falsify sudden freeze-out by examining whether particle 
abundances and spectra can be described using the same formation temperature. 
This is done both in the chemical equilibrium framework, and using a chemical non- 
equilibrium ansatz. Our fits to experimental data suggest that the sudden freeze-out 
model explains both the particle abundances and spectra. 

We then try to extract the particle formation temperature, and quantify 
post-freeze-out hadronic interactions using experimentally observable resonances. 
We discuss observed resonances and suggest further measurements that have the 
potential to distinguish between the possible freeze-out scenarios experimentally. 
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Finally, we infer from experimental data how particle formation proceeds 
in spacetime, in particular whether freeze-out dynamics agrees with the sudden 
freeze-out expectation. We examine particle spectra, and show that they are not 
sensitive enough to pick the right freeze-out dynamics. We suggest resonances and 
azimuthal anisotropy as experimental probes for this task. 
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CHAPTER 1 
The strong interaction and quark-gluon plasma 

1.1 Elementary and collective phenomena 

Physics research within the last fifty years can be divided into two classes of prob- 
lems: figuring out the fundamental rules which govern physical processes at the 
microscopic level, and examining how these fundamental rules give rise to the phe- 
nomena we observe macroscopically Particle physics is concerned with finding the 
fundamental rules. Condensed matter physics studies how these rules generate 
complex structure. 

This division is somewhat simplistic. Quantum mechanics, together with 
relativity, imply that every interacting system will always have an infinite number of 
degrees of freedom. Even the vacuum state will always contain particle-antiparticle 
pairs popping in and out of it at all times. This means that the vacuum, just like 
any material studied in condensed matter physics, may exhibit collective phenom- 
ena. If the interaction is strong enough, the collective phenomena will dominate 
"elementary" physics, and the observed degrees of freedom might be, in their nature 
and behavior, different from the fundamental ones. In other words, to understand 
the physics of one particle, we need to understand the collective medium through 
which this particle will propagate. 
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As an arena to understand general collective systems, condensed matter 
physics has some limitations. While the structures and properties of condensed 
matter systems are enormous in their variety, the "elementary rules" they are based 
on are limited and understood: electromagnetism and quantum mechanics. It is not 
at all certain that by introducing a different set of fundamental rules we will not 
arrive at a qualitatively novel macroscopic system. The fact that an analogue of 
strong interaction confinement has never been found in a condensed matter system 
reinforces the need to investigate how the "macroscopic phenomena" change when 
a qualitatively different set of "microscopic rules" is introduced. 

The purpose of heavy ion physics is to do just that: to create and study a 
"condensed matter system" using a different fundamental interaction than electro- 
magnetism. 

1.2 The strong interaction 
1.2.1 QuEirks 

According to the current theory of the strong interaction, all hadrons are composed 
of fermionic constituents called quarks. The non-interacting quark Lagrangian 
therefore comprises several massive fermion fields, denoted by a quantum number, 
called flavor (/), conserved in strong interactions 



Quarks are arranged in three flavor doublets, of which the top member has electro- 
magnetic charge -1-2/3 and the bottom member has charge —1/3 (in units where 
the proton has charge -|-1). 
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The quarks in the first doublet, called up and down {u and d from now on), have 
a mass that is negligible with respect to the mass of the observed hadrons, and are 
therefore often called "light quarks" . These quarks are the only ones stable against 
flavor-changing weak decays. 

The quarks forming the second doublet, called strange and charm (s and c) 
are more massive. The strange quark is of particular interest to the subject of this 
dissertation, since, as we will see, its mass (50 to 150 MeV) is comparable to both 
the mass of the hghtest hadron and the energy scale of the QCD phase transition. 
The charm quark is much more massive (~ 1 GeV), though it may be produced in 
the initial processes of a heavy ion collision. The third doublet, top and bottom, is 
too massive to be produced in present heavy ion collision experiments (though the 
~ 4 GeV b quarks could be studied in future heavy ion collisions). 

1.2.2 Color and gluons 

The strong interaction, responsible for holding hadrons and nuclei together, is funda- 
mentally very similar to QED in that interactions maintain the Lagrangian invariant 
under a local phase transformation of the quark wave function 



However, we postulate that quarks transform as vectors in a 3-dimensional space 
called color space, while anti-quarks are co-vectors 



(1.2) 



UU+ = 1. 



(1.3) 
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Hence, the most general form of Eq. (1.2) allows C/ to be a special unitary matrix. 
{UU^ = 1 , = 1). Such matrices form a group, called SU{2i) ^. It can be shown 
that for any U eight real numbers aj can be chosen such that 



[lA) 



where the generators of the SU (3) group tj can be represented by the Gell-Mann 
matrices 
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Just as in the case of electromagnetism, imposing the phase symmetry locally re- 
quires the introduction of a covariant derivative to subtract the effect of a varying 
U on neighboring points. This covariant derivative, however, now contains eight 
terms corresponding to each Gell-Mann matrix 



(1.6) 



Physically, substituting into Eq (1.1) means that quarks can interact through 
eight bosons called gluons. However, unlike QED, each gluon also carries a color 



^Thc results described in this section can be generalized to SU{N) invariance, corresponding 
to physically N colors and A''^ — 1 generators. Historically, the first example of this theory to be 
examined was the Yang-Mills Lagrangian, corresponding to SU{2). 
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charge, which, in the Gell-Mann representation Eq. (1.5) corresponds to 

rb + bf rb — br rr — bb r'g + gf rg — gr bg + gb b'g — gb rf + 66 — 2g'g 
^ ' ' V2 ' V2 ' i\/2 ' V2 ' i\/2 ' 2^3 

for ^1. 8 respectively. 



:i.7) 



Gauge invariance also allows for a dimension four gluon field contribution 
to the Lagrangian, equivalent to a curvature tensor in gauge space 

8 ^ 8 

L ^ L + ^ IV D^] [D^ D^]) ^L+-Y^ F,^,F^\ (1.8) 

A qualitative difference with respect to QED emerges due to the fact that the matri- 
ces ii_8 do not commute. Their commutators, represented by the SU{3) structure 
constants fj,ji,j2 

8 
i=l 

will therefore appear in the QCD field strength since 

8 8 

[D,D,] = ... + [tntjM^Al = ... + i Yl fn.nhKK (l-^O) 

jlj2 = l j,Jlj2=l 

Hence, gluons can also interact, and their interaction manifests itself in a physical 
gauge- invariant way through a term in the gluon field strength. 



In summary the Lagrangian of interacting QCD is 

3 



j=l f c,c=l 

with 



3=1 



i^cj (1.11) 



p,^ = - d-^A^ - g Y fmn^A- (1-12) 

Here, the index c runs over the color quantum number, while j runs over the gluon 
type. The gluon matrices tis contract the anti-quark co-vector with the quark 
vector. 
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While the symmetry principles used in deriving the QCD Lagrangian follow 
the QED case closely, the self-interaction term has profound consequences not seen 
in QED. This becomes evident when the theory is quantized. 

1.2.3 Quantization and the running coupling constant 



Expanding Eq. (1.11) in terms of A^j^j leads to the following Feynman rules'^ 



P ,J1 



/ MO 



Jl 




PJ3 T,J4 



(1.13) 
(1.14) 

(1.15) 



+fnn3'U.Mfi9"'9"'-9^9''n 

Where the indices correspond to color for quarks {a,b — 1..3), type for gluons 
{ji.A — 1--8) and Greek letters for gluon Lorentz indices. 

While the first term is similar to the QED case, the next two, which come 
out of the /afec^fe contribution to F^'' are new. Their profound effect can be seen 
when one attempts to calculate quantum corrections to the coupling constant g {/i 
is the momentum transfer at which the coupling constant is measured) 

9\f^') 



9\<f) 




- 1 - [n(g2) - n(/.2)] 




(1.16) 



(1.17) 



^We will ignore the additional "ghost" particle needed in some gauges. See (Peskin and 
Schroeder 1995) for more details. 



23 

The explicit calculation of this ll{q^) is a somewhat involved endeavor beyond the 
scope of this thesis (Peskin and Schroeder 1995). However, one can see that since 
boson loops have the opposite sign from fermion loops, the contribution of the last 
two diagrams will in general have an effect on the running of the coupling constant 
opposite to that of the first term. Intuitively, while the first term screens the charge, 
the gluon loops increase the field strength while maintaining the attractive potential, 
resulting in a larger effective charge than the bare one (anti-screening). Combining 
all the diagrams, one can see that in a general Yang-Mills theory with N colors and 
n light flavors, the coupling constant will be, to a leading order correction (Peskin 
and Schroeder 1995)^ 

Provided a scale /x^ can be found where g"^ « An this is a self-consistent result. 

In the QED case {N — 0,n = 1 for the light electron) a{k'^) increases 
at increased momentum transfer, and it is found that in low energy experiments 
= « = « 1. Hence, 

««-(*^) = i-i4.lM^v<) 

with the theory remaining perturbative until a very high energy scale (~ e^°°me) 
at which QED is not expected to hold anyways. 

In the QCD case, however {N — 3,n = 3 for three light quarks) the sit- 
uation is reversed, decreases with increasing momentum transfer, and even- 
tually becomes perturbative at a large scale A ( the experimentally measured 
^higher order corrections are generally not thought to change the direction of the running. 
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g^(lOOGcV) ^Qi2) 
47r ■ / 

The highly energetic quarks produced at this energy scale become bundles of approx- 
imately coUinear fast particles known as jets, with each jet having the momentum 
of the original quark. 

However, at smaller increases to the point where perturbation theory 
ceases to be a useful tool. Thus, while high-energy scattering cross-sections can be 
computed through perturbative QCD, low-energy effects such as vacuum structure 
and infrared corrections to particle propagators are governed by non-perturbative 
physics. 



1.3 The structure of the QCD vacuum at zero temperature 

The non-perturbative nature of the strong interaction in the low-energy limit should 
not come as a surprise: the fundamental degrees of freedom of QCD (quarks, gluons) 
are very different from the fundamental degrees of freedom observed in the strong 
interactions at low energies (the 100-1- hadrons in the particle data book). For 
instance, every hadron observed so far is not charged under color, and does not 
interact manifestly through gluon exchange. It must be the case, therefore, that 
somehow only color singlets survive as physical states. Moreover, QCD predicts 
gluons to be massless, and high energy scattering processes have determined light 
quarks to be extremely light (~ a few MeV). Yet the observed hadrons are much 
heavier than that (~ GeV). Both of these phenomena are not as yet rigorously 
understood. However, there are phenomenological ways to see how they arise. 
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1.3.1 Confinement 

A simple picture of what happens when two color charges (for instance, a quark- 
antiquark pair qq) become separated is provided by relativistic quantum mechanics 
(Gribov 1999): As the separation of the qq pair increases, so does the effective 
charge. At a certain point it becomes energetically convenient to create a new qq 
pair so that a separated qq becomes two tightly bound qq. 

An analogous situation exists in atomic physics: Dirac energy levels of an 
electron orbiting a Z > 137 point charge become complex (Rafelski et al. 1978) 
as it becomes energetically possible to create a real e'^e" pair out of the vacuum ( 
Az —>■ Az + e"*" + e~ ). The positron escapes to infinity, and the energy contained 
in the {Az + e~)bound + e"^ state is less than that of the field of the isolated Az ion. 
The anti-screening of QCD means that the field strength necessary for this vacuum 
instability will be reached when any two color charges become sufficiently separated 
(Gribov 1999) . 

This picture is reinforced by evidence from lattice gauge theory (Wilson 
1974), a calculational technique valid when asymptotic freedom applies. In this 
approach, the continuous spacetime is replaced by a discrete lattice of points, as it 
is assumed that smaller-scale degrees of freedom become perturbative. It is then 
possible to represent the pure-gauge action using links from one lattice point to the 
next 

C/(xi,a;2)„~e'^^-T'^^''^^*'^ (1.21) 

and calculate physical observables in euclidean spacetime {t —>■ it) by evaluating 
a trace over the lattice points (which approaches the functional integral in the 
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continuum limit) using Monte-Carlo techniques. For instance, the potential energy 
between two static color sources 



a=l a=l 

corresponds to the expectation value of a loop in Euclidean spacetime spanning the 
particle's separation r (Wilson 1974; Bali 2001). 



< V{r) >= lim -log(< W{r) >] 

i— >oo t 



< W{r,t) >-- 



J dAexp 




J dA exp 







[1.22) 



1.23) 



In the strong coupling limit and neglecting quark loops ^ it can be proven (Creutz 

T 




-1 
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Figure 1.1: Potential calculated with full three flavor QCD , rescaled by a constant 
ri at which confinement effects become significant. 



1978) that 



V{r) oc r. 



:i.24) 



^This is known as the quenched approximation. Without it, the potential wiU saturate at the 
value necessary to create a quark- antiquark pair, as per the Gribov mechanism. 
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A numerical calculation with a more realistic picture, with three quark flavors, does 
not change this conclusion; as Fig. 1.1 (Bernard et al. 2000) shows, the effective 
QCD potential between heavy quark charges can be fltted by a function increasing 
linearly with r. Heavy quark spectoscopy (cc J/^ states and bb T states) has ex- 
perimentally confirmed that the QCD potential can indeed be effectively modeled 
by a coulomb and a hnear term (Richardson 1979; Brambilla et al. 2000). 

The Gribov argument and lattice evidence point to a coherent qualita- 
tive picture of what happens as the strong interaction coupling constant increases 
beyond the perturbative limit. Instead of decreasing at large distances, as the elec- 
tromagnetic potential does, the strong interaction potential continues to increase 
monotonically until the vacuum decays into more quarks and gluons in such a way 
as to make every state ultimately observed a color singlet. 

This phenomenon, called confinement, has not been fully understood. 
Nothing similar has, as yet, been observed within many-body systems interact- 
ing via electromagnetism. Thus, the study of how confinement arises is a good 
example of the problems outlined in section 1.1. 

1.3.2 Chiral symmetry breciking 

A considerable number of hadron masses (p,nucleons,A,A,E,E*,S,S*,Q) can be fitted 
reasonably well by assuming hadrons to be made up of "constituent" quarks, whose 
number is given by the hadron 's fiavor content. If the constituent quark mass is set 
to ~ 300 MeV for u, d and ~ 500 MeV for s, observed masses for these hadrons 
can be described in terms of the constituent quark masses plus spin couplings and 
electromagnetic effects. 
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To see how such constituent quark states might arise, consider decomposing 
the hght quark part of the Lagrangian into left-handed and right-handed compo- 
nents 



d 

\ ' J 



[1.25) 



V'L = ^(1-75)V' (1-26) 
V^i?=^(l + 75)V' (1.27) 

The quark part of the Lagrangian then becomes 

Lq = ^i^Di^il^L + ^i^Di^.i/jR - m^.{iJ^i{jR + ^r^l)- (1-28) 

i=u,d,s 

We see that the only mixing term is provided by the quark mass. Therefore, in 
"bare" QCD, the left and right light quark currents should be separated with only 
a small correction (~ lOMeV for u, d, ~ lOOMeV for s). 

However, if the QCD vacuum exhibits quark condensation 

< vac\'il)'il)\vac >—< vac\%l)Li^R -\- ipRipilvac >^ 

the low energy effective Lagrangian will exhibit an additional quark mass, since, as 
implemented in effective models such as Nambu-Jona-Lasinio (Nambu and Jona- 
Lasinio 1961; Nebauer and Aichelin 2002), interactions between quarks will generate 
an effective mass term. 

We can see that this is the case in QCD since the Lagrangian in Eq. (1.28) 
is, in the mu,d,s ~^ limit, invariant under two separate symmetries 
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where / is the identity matrix, Uj are the Pauh matrices and (y.Li,Rhl3L,R are real 
numbers. Using the transformations of Eq. (1.29), corresponding to the 

C/(3)l X [/(3)i? = SU{'i)L X 5t/(3)i? X U{1)l x U{1)r 

= SU{'i)v x SU{'i)A X C/(l)y X C/(1)a 

symmetry group, it should be possible to transform any hadron into a hadron of 
inverse parity (e.g., a negative-parity vector meson, with a current given by ip^f^ip, 
into a positive parity pseudo- vector ip^fj^^y^ip) . Hence, the mass difference between 
such states (e.g. p and /i(980), or the A and the A(1405)) should be of the order 
of the quark mass. In fact, these differences are ~ SOOMeV, the same order of 
magnitude as the constituent quark mass. 

Hence, it appears that U{3)l x U{3)r is spontaneously broken to a lower 
subgroup. Identifying the lower subgroup as SU{3)v x f/(l)v (IsospinxBaryon 
number)^ we are left with eight Goldstone modes, the 7r,K,r) mesons (which are 
all considerably lighter than the constituent quark prediction, although the large 
strange quark bare mass makes the K and rj much heavier than the tt). 

The nature of the chiral symmetry breaking, and its relationship with con- 
finement, are under intense study. However, numerical results described in the next 
section suggest the two phenomena are related. 

1.4 The QCD phase transition 

The previous two sections make it clear that in empty space light colored particles 

become confined in massive (at least 140 MeV) color-neutral composite particles. 

^U{1)a, is broken through quantum corrections ('t Hooft 1976), which is why no hght isospin=0 
pseudo-scalar meson exists, the rj' having a mass of IGeV. 
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The question which immediately follows is: What happens if temperature or quark- 
density increases to such an extent that inter-hadron separation becomes compara- 
ble to the separation of quarks within a hadron? Intuitively, hadrons should then 
break down, and at the very high temperature or density limit the system should 
approach the asymptotically free limit of a relativistic ideal gas. However, does 
this occur as a phase transition (ie, a discontinuity in some order parameter) or is 
it rather a continuous cross-over, much like the formation of the electromagnetic 
plasma? At what temperature do the quark and gluon degrees of freedom start 
manifesting themselves? What does a strongly interacting QGP look hke away 
from the infinite temperature limit? 

Studying these questions, at both theoretical and experimental levels, is 
important for a variety of reasons, some of which are discussed in section 1.1. The 
Universe underwent this phase transition shortly {t ~ 10~® seconds) after the big 
bang, and evolved as a quark-gluon plasma before that time. Hence, the QCD 
phase transition might have had consequences for the evolution and structure of 
the Universe. Moreover, there are astronomical objects today, such as neutron 
stars, where a QGP might be expected to exist naturally. Finally, to reiterate 
the point made at the start of this chapter, the QCD phase transition (if it is 
indeed a phase transition) is the only vacuum phase transition which is at energies 
experimentally accessible for the foreseeable future (It is very unlikely, for instance, 
that we will be able to create any-time soon a thermalized system hot enough 
for electroweak symmetry restoration). The study of the QCD phase transition, 
therefore, is invaluable for our understanding of quantum field theory in general, 
especially in the regime where standard perturbation theory methods fail. 
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1.4.1 A rough calculation (which will prove useful later) 

We shall proceed to calculate the QCD phase transition parameters from the MIT 
bag model (Johnson 1975), an ansatz which can fit reasonably well the spectra of the 
observed non-Goldstone hadrons. The zero-temperature QCD vacuum is assumed 
to have a "bag constant" term due to deconfinement/chiral symmetry breaking, i.e. 
a vacuum energy with positive energy density (e) but negative pressure (P). The 
energy-momentum tensor is then 

T^'-TZtter + Bg^" (1.30) 

with 

e = ^matter + B (1-31) 

P — Pmatter ~ B (1-32) 

Hadrons arise as "bags" filled with free or weakly interacting quarks, trapped in 
volumes where the pressure exerted by the quarks in the bag balances the negative 
pressure exerted by the bag. 

Through solving the Dirac equation, one can find the quark wave-functions 
inside the bag and the hadron masses in terms of the bag constant. Fits to the p 
meson, nucleon and A yield a bag constant of B^^^ — ISOMeV (Chodos et al. 1974). 
It can be seen that in this model a first order phase transition may occur when a gas 
of hadrons has the same pressure as a quark-gluon plasma (Rafelski et al. 1996). 
We calculate the pressures of the particles using the grand canonical ensemble for 
a non-interacting ideal gas (see the next chapter for details) 
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where ± applies, respectively, to fermions and bosons and gi is the degeneracy. 
For the QGP phase we used 2 massless and one massive quark (fermions) with 
a degeneracy of 6 (spin x color) and eight gluons (bosons) with a degeneracy of 2 
(polarizations). For the hadron gas, in addition to the vacuum pressure given by 
the bag constant, we used all particles in the particle data book with a mass lighter 
than 2 GeV (Hagiwara et al. 2002). Light flavor chemical equihbrium was assumed 
{/iq — —/iq) and the strangeness chemical potential was put to (no strange quarks) . 



Fig. 1.2 shows the phase diagram calculated in this ansatz. As can be seen, 
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Figure 1.2: The bag model phase diagram . 

our estimate predicts a critical temperature roughly equivalent to the mass of the 
lightest hadron. While many assumptions used in this estimate can be considered 
naive at first glance, calculations on the lattice give very similar results. 
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1.4.2 Results from lattice QCD 

A phase diagram calculation can also be performed from first principles QCD 
through lattice techniques, introduced in section 1.3.1. The fact that lattice field 
theory can be used to describe finite temperature systems is evident from the fact 
that the QFT partition function 



^QFT 



J B^DAe'^^^'^^^'^'^^'"^^ (1.34) 

is related to the thermodynamic partition function 

Zth = J DijDAe'^^^^'^^'^'''/^ (1.35) 

by the substitution t ^ ij;. This substitution, together with the fact that — 
^A+2m^ means that a quantum field theory with a periodic boundary conditions in 
the time direction (anti-periodic for fermions) (Matsubara 1955; Kapusta 1979) 

^(^'^^^) (^-^^^ 
A(^x,t+'^^^A{x,t) (1.37) 

will effectively model a finite-temperature quantum field in equilibrium with a heat 
bath. 

In such a formulation the Wilson loop definition given in Eq. (1.23) has 
to be modified, since no t — > oo limit is possible. Instead, the relevant variable 
becomes the expectation value of the Polyakov Loop (Polyakov 1978; McLerran 
and Svetitsky 1981), where the quark propagator closes on itself in the periodic 
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time: ^ 



< Wt^T^t) >-- 



J dAexp 


ij.Af^-L) dtd'x 


J dA exp 


— ^^'q^ Ldtd'^x 





:i.38) 



Its far from obvious, however, if this calculation will be relevant to the 
physical world. Eq. (1.35) refers to a system in perfect equilibrium, something 
which did not apply either to the expanding early Universe or to the conditions in 
which quark-gluon plasma is created. Nevertheless, a lattice calculation can serve 
as an indication of what we can expect the quark-gluon plasma phase transition to 
be like. 
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Figure 1.3: Left: The Polyakov loop expectation value as a function of distance at 
several temperatures, a refers to the best fit zero temperature flux tube tension. 
Right: The order of the deconfinement phase transition as a function of the number 
of flavors and quark masses (physical QCD is somewhere in between the extrema 
of the diagram). 



Fig. 1.3 (Karsch and Laermann 2003) shows that the lattice does indeed 
exhibit something which looks like a phase transition. As the temperature increases, 

the distance dependence of the Polyakov loop decreases, to vanish at a critical 

^Further evidence that the Polyakov loop is the order parameter for deconfinement is the fact 
that it is not gauge-invariant for color non-singlets (McLerran and Svetitsky 1981). Hence, unless 
its expectation value becomes independent of quark separation (which it does above a certain 
temperature as shown in Fig. 1.3 ), only color singlets survive as physical states above the lattice 
scale. 
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value as expected in a deconfining phase. We still, however, do not know what 
order characterizes this phase transition, or even if a phase transition is present 
at all. As the right panel of Fig. 1.3 shows, the order of the phase transition is 
strongly dependent on theory input, such as the number of flavors and the values 
of the masses of the light and strange quark. Since performing calculations using 
realistic inputs for the light and strange quark masses (the realistic values for these 
quantities are somewhere in the middle of those in the panel) is beyond the range of 
today's computing power, the order of the phase transition is still unknown. Fig. 1.4 



200 r 




1.0 1.5 2.0 2.5 3.0 3.5 4, 



200 400 

^q(MeV) 



600 



Figure 1.4: Left: Energy density over (should be fiat for a relativistic ideal 
gas), calculated on the lattice as a function of temperature, scaled with the critical 
temperature (where the sharp step occurs). Right: The phase diagram on the 
lattice. The diagonal dotted line shows the limit beyond which the curvature of 
p/T^ can not be fitted. The diamond refers to the estimated critical point. 



(right) (Ejiri et al. 2003) shows that the phase diagram mirrors very closely that 
we obtained in the bag model picture of Fig. 1.2. To compound this. Fig. 1.4 (left) 
shows that, above the rather sharp critical temperature jump, the thermodynamics 
of the strongly interacting system seems to be close to that of a relativistic ideal 
gas with perturbative corrections (Letessier and Rafelski 2003). 

In fact, as calculated in (Karsch et al. 2003), the hadron gas phase seems 
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to be modeled closely by the resonance gas ansatz used in this thesis, while pertur- 
bation theory to a few orders can describe the QGP just above critical temperature 
reasonably well (Letessier and Rafelski 2003). Fig. 1.5 (Karsch 1998) can perhaps 
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5.40 



Figure 1.5: Polyakov loop susceptibility (left) and i/ji/j susceptibility vs (~ T) 
for various quark masses (expressed in terms of the inverse of the lattice size) . 

provide an indication of why this is the case. As can be seen, the discontinuity in 
the Polyakov loop susceptibility xl and the ipip susceptibility Xm seem occur at the 
same temperature for a range of quark masses. Hence, the onset of deconfinement 
also means that the quarks lose their dynamically acquired "constituent" mass and 
become nearly massless. This, together with the color screening in a deconfined 
medium, means that the QGP can be described by the ideal ultra-relativistic gas 
ansatz at temperatures much lower than originally expected. 
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1.5 How to study the quEirk-gluon plcisma experimentally 

We now proceed to give an overview of the main experimental issues raised in the 
study of QGP. The only way we can think of to make QGP in a laboratory is 
to collide two energetic heavy ions (ideally large enough to make an equilibrated 
system described by the grand canonical ensemble). Hopefully, the collisions will 
transform some of the energy of the heavy ion into heat, resulting in a hot hadronic 
system in the center of mass frame. If enough energy is produced, the system will 
then undergo a phase transition, evolve as quark-gluon plasma, change back into 
hadrons, and we will be able to tell that a phase transition has taken place and to 
study the properties of deconfined matter through a careful analysis of the decay 
products. 

It is easy to see that in practice telling whether quark-gluon plasma has 
been produced is not a simple matter. Indeed, the two claims of quark-gluon plasma 
production ( (Heinz and Jacob 2000) and (Gyulassy 2004) ) are both based on an 
assessment of different kinds of signatures. (And, it should be said, the kinds of 
signatures claimed by (Heinz and Jacob 2000) differ significantly from those used 
in (Gyulassy 2004)). Table 1.1 summarizes the experimental heavy ion collision 
programs energies and nuclei. Some accelerators are geared toward fixed-target 
experiments and others are colliders. The advantage of the collider is that much 
more energy is available to make a thermalized system. The disadvantage, aside 
from the cost, is that it is a lot more difficult to have an acceptance covering a large 
solid angle (rather than some region, for example mid-rapidity). Table 1.2 provides 
a list of some heavy ion experiments. In the next subsections, we will proceed to 
give a summary of the issues involved in each experiment. 
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Table 1.1: Heavy ion Experimental program. 



Accelerator 


Type 


System 
under study 


Energy (GeV/A) 
(Center of Mass) 


Place 


SIS 


Fixed target 


Various 


2-4 


GSI, Darmstadt 


AGS 


Fixed target 


Au-Au, S-Au 


2-11 


BNL, NY 


SPS 


Fixed target 


C-C, S-S, Pb-Pb 


40, 80, 158 
(8.73, 12.3, 17.3 ) 


CERN, Geneva 


RHIC 


Collider 


Au-Au 


19.5, 130, 200 


BNL,NY 


LHC (2007) 


Collider 


From p-p to Pb-Pb 


7000 


BNL,NY 


GSI (200...) 


Collider 


High //^/density 


In development 


GSI, Darmstadt 



Table 1.2: Selected heavy ion Experiments. "Telescope" refers to a detector narrow 
in both angle and rapidity, "mid-rapidity" to a detector covering rapidity region 
where the center of mass momentum vanishes. "Wide acceptance" means an area 
outside the mid-rapidity is covered. 



Name 


Location 


Predominant phase space 


What it measures 


KaoS, FOPI 


SIS 


Wide acceptance 


Hadrons, Strangeness 


E864,E878,E886 


AGS 


Wide acceptance 


Hadrons, Strangeness 


NA57 


SPS 


Telescope 


Strangeness 


NA49 


SPS 


Wide acceptance 


All hadronic signatures 


NA50, NA60 


SPS 


Telescope 


leptons, J/^ 


WA98 


SPS 


Telescope 


7 s, hadrons 


CERES/NA45 


SPS 


Mid-rapidity 


7 s, leptons 


STAR 


RHIC 


Mid-rapidity 


Hadrons, jets, strangeness 


PHENIX 


RHIC 


Mid-rapidity and telescope 


Hadrons, jets, leptons, 7 s 


PHOBOS 


RHIC 


Wide, mid-rapidity parts 


Hadrons, jets 


BRAHMS 


RHIC 


Tc4osc()])e. wide rapidit}' 


Hruh'oiis. jots, strangeness 


ALICE 


LHC 


Mid-rapidity 


Everything 


CMS 


LHC 


Mid-rapidity 


J/^, T, 7 s, jets, leptons 


ATLAS 


LHC 


Mid-rapidity 


T, jets 
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1.5.1 Compcirison benchmcirks 

The first issue to think about is a comparison benchmark. To test for quahtatively 
new physics, the system under consideration has to be compared to a system in 
which the physics can be understood in terms of phenomena encountered previously. 
In our case, this means a system with low energy, or very few hadrons, in which 
we know tat quark-gluon plasma can not be produced, and the dynamics can be 
understood either by perturbative QCD or as a superposition of hadron-hadron 
collisions. The latter can be analyzed in terms of elementary collisions measured in 
control experiments, either through bulk analytical calculations (e.g. the wounded 
nucleon model (Ftacnik et al. 1987)) or microscopic kinetic models (e.g. Quantum 
molecular dynamics, or uRQMD (Bass et al. 1998) or hadronic string dynamics, or 
HSD (Ehehalt and Gassing 1996)). 

The most obvious choice for a benchmark experiment is proton-proton col- 
lisions. Any effects pecuhar to Nucleus-Nucleus {A — A or A — B) coUisions can be 
parametrized through the nuclear modification factor (normally referred to as Rab, 
sometimes as "enhancement" or "suppression" , depending on what the measured or 
expected physics looks like): the ratio of a quantity observed in the nuclear collision 
to the same quantity observed in a p — p collision normalized by the product of the 
two nuclei mass numbers A x B, or by the number of participants. 

Proton-proton collisions might not always be a good benchmark, since there 
can be large volume or many-body effects which are not associated with the phase 
transition but still have to be accounted for. Hence, proton-nucleus, deuteron- 
nucleus or collisions between light nuclei are also studied. In these, there is a large 
volume of excited matter and many collisions, but (presumably) no phase transition. 
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More generally, it is possible to analyze observables in terms of collision centrality 




Figure 1.6: Left: The definition of impact parameter in a colliding system. Right: 
How energy is distributed after the collision. "Stopping" refers to the dissipation 
of the system's initial longitudinal momentum. 

(number of participants) or collision energy. This way, one can investigate if there 
is a critical energy (roughly correlated with temperature) or number of participants 
(roughly correlated with reaction volume and evolution time) at which qualitatively 
new behavior (e.g., a phase transition) occurs. A measure of the collision centrality 
is provided by the number of participating nuclei, related to the impact parameter, 
6, Fig. (1.6 Left), b can be found either directly, by having a detector placed 
downstream from the collision point to measure how many nucleons escaped the 
primary collision, or indirectly, by simply measuring the total multiplicity (which 
is related to the volume of the thermal system) . 

It should be noted, however, that increasing the coUision energy leads to 
some non-trivial consequences which should be taken into account. First of all, 
higher collision energy increases the role of perturbative QCD within the system. 
In the non-perturbative limit, the particles will undergo repeated soft interactions 
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after the primary collision. These interactions will stop the particles in the collision 
region and transform their initial colliding energy into thermal energy. Hence, the 
size of the colliding system, and ultimately the total number of particles produced, 
will scale as the number of participants 



As the collision energy increases, perturbative interactions will play a bigger role. 
This means that collisions, rather than participants, will play a larger role in stop- 
ping the system, and the scaling becomes (parametrized by q;,/3) 



Due to the large center of mass momentum and asymptotic freedom, stop- 
ping in primary collisions will decrease with increasing energy. Two highly energetic 
heavy ions will pass each-other nearly transparently. Their energy will be released 
into the medium due to confinement, as the color fields generated by the initial 
collision interactions become stronger and ultimately melt into incoherent pairs 
(see Fig. 1.6 (right). The two highly Lorentz-contracted nuclei can be thought of as 
capacitors of the color field). For this reason, heavy ion systems colliding at high 
energy also tend to have a lower baryon density and chemical potential. 

These considerations also control the likely initial conditions for hydrody- 
namic evolution. A system with high baryon stopping will evolve from the pancake- 
shaped collision region, localized in rapidity space around the center of mass frame. 
This system has been originally analyzed as a starting point for hydrodynamic evo- 
lution by Landau (Landau 1953). In the infinite time limit, it should give rise to an 




,s 



(1.39) 



prod ^-^participants ~l~ P -^collisions 
collisions ~ ^^participants' 




(1.41) 
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elliptical fireball with comparable transverse and longitudinal flow proflles, which 
at mid-rapidity can be well approximated by a spherically symmetric freeze-out. 

By contrast, the thermalized system formed when the two very energetic 
nuclei pass through each other will have a high degree of boost-invariance. As 
argued in (Bjorken 1983), two highly Lorentz-contracted "pancakes" passing nearly 
transparently through each-other will look the same in a boosted reference frame. 
Hence, the dynamics of this collision should exhibit a rapidity plateau. As shown 
in (Bjorken 1983) this initial condition is hydro-dynamically stable. 

Fig. 1.7 shows that the two limits can be used as an ansatz for, respectively, 
SPS and RHIC energies. At the SPS (Appelshauser et al. 1999), the pseudorapidity 
distribution is approximately Gaussian, corresponding to a sharp longitudinal struc- 
ture. Hence, we have used a spherical flow proflle to model SPS mid-rapidity data 
in chapter 3. As the right panel of Fig. 1.7 shows, the situation is perhaps different 
at RHIC. We have used the Monte- Carlo described in detail in chapter 4 to verify 
that the flat pseudorapidity distribution observed by PHOBOS (Olszewski et al. 
2002) can be reproduced by a boost-invariant source. Such a source is described in 
detail, and used to model RHIC data, in chapter 4^. 

1.5.2 Jet quenching 

Jets, or energetic streams of hadrons resulting from a perturbative QCD interaction, 
are a promising test of whether the dense system formed in a heavy ion collision 

cxriil)its liadroiiic or ])artoiiic degrees of freedom. A parton propagating through 

■^Recent data from BRAHMS (Bearden 2004) may have thrown the vaHdity of this picture 
in doubt, as the rapidity distribution of identified tt is approximately a Gaussian, similar to 
(Appelshauser et al. 1999) and predicted by the Landau model. 
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Figure 1.7: Left: SPS 19.4 GeV/A charged particles pseudorapidity distribu- 
tion measured by the NA49 collaboration Right: RHIC 130 GeV/A charged 
particle pseudorapidity distribution measured by PHOBOS, and compared to a 
boost-invariant statistical model (see chapter 4). 

a quark-gluon plasma can lose a lot of energy quickly due to repeated interactions 
with soft gluons. This energy loss is especially pronounced due to the quantum 
interference of interactions between many coherent gluons (this is known as the 
Landau- Pomeranchuk-Migdal effect) (Wang et al. 1995). The quark-gluon plasma 
is therefore expected to be much more opaque to jets than ordinary nuclear matter. 
One way to quantify such opacity is to measure jet azimuthal correlation, and 
compare it to elementary (proton-proton or proton- nucleus) collisions. As Fig. 1.8 
shows, in an opaque system one of a pair of jets should become quenched, and hence 
unobservable. For this reason, jets should lose the expected azimuthal correlation. 
In case of non-central collisions, this effect will also acquire an angular dependence 
(Gyulassy et al. 2001). 

As Fig. 1.9 (left) shows, this effect has received a remarkable experimental 
confirmation (Adams et al. 2003a). The observed loss of azimuthal correlation is 
qualitatively unique to Au-Au collisions, the jets in both p-p and d-Au collisions 
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Figure 1.8: How an opaque system leads to jet correlations . 

being perfectly correlated. As Fig. 1.9 (right) shows, the loss of azimuthal corre- 
lation is accompanied by the loss of hard (jet-energy) particles, with respect to an 
extrapolation of hard particles produced in p-p collisions. This is also unique to 
Au-Au, hard particles in d-Au being enhanced. 

While some features of the data have recently been modeled within a 
hadronic scenario with strong in-medium modifications (Capella et al. 2004), the 
full high pt dataset from p-p, Au-Au and d-Au collisions has only been convinc- 
ingly described in a model which assumes that the matter at the center of Au-Au 
collisions exhibits colored degrees of freedom and is 100 times the density of normal 
nuclear matter. Furthermore, the non-decrease of the jet suppression with px could 
only be modeled by taking the Landau- Pomeranchuck-Migdal effect into account 
(Gyulassy et al. 2003). 

While the evidence described above makes jet quenching widely regarded 
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Figure 1.9: Left: jet azimuthal correlation (number of jet pairs emitted at an an- 
gular separation A0) normalized by the number of triggered jets, as measured by 
the STAR collaboration in p-p, d-Au and Au-Au collisions and 200 GeV. Right: 
The nuclear modification factor (as defined in section 1.5.1) for d-Au and Au-Au 
collisions, measured by STAR. The curves refer to calculations done within a model 
where partonic thermalized degrees of freedom are assumed. Minimum bias means 
that no centrality selection was performed. 

as the definitive proof of QGP formation (Gyulassy 2004), jet quenching as a QGP 
diagnostic has limitations. The information it provides about QGP equilibration, 
kinetics and equation of state is rather limited. Jets will not evolve with the bulk of 
the matter, and their hadronization, due to asymptotic freedom, will not impact the 
observed energy-momentum distribution significantly. While jets can sample some 
extensive quantities, such as the density of matter they traverse, soft physics is 
needed to understand if the system having those extensive quantities is equihbrated 
in a particular phase, or if it is in a different phase from the usual hadronic matter 
(Wang et al. 1995). Soft physics is also needed as a complement to perturbative 
QCD to account for the missing jet energy. Modeling quantitatively how the jet 
energy is distributed among the soft degrees of freedom, and what the effect of the 
deposited energy is, is an interesting theoretical problem potentially rich in insights 
into fundamental physics. However, an understanding of the soft degrees of freedom 
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is essential to solve it. 

Hence, while jet quenching might convince us that QGP is there, it is not a 
very good tool for measuring its properties and how it changes into normal matter. 
For this, probes more dependent on soft physics are required. 

1.5.3 Direct photon and dilepton production 

Historically, the standard way to measure the temperature of a hot medium has been 
to measure the thermal spectrum of its emitted light. This can be expected to be 
true for a QGP as well, given the abundance of scattering reactions with quarks and 
gluons (Kapusta et al. 1991). Of course, the processes under consideration here are 
energetic enough to allow the emitted photons to produce lepton-antilepton pairs 
(Kajantie et al. 1986). Photons (7) and leptons antilepton pairs (11) will be emitted 
in reactions such as 



The particularly attractive feature of this signature is that photons and dileptons 
always emerge from the system without undergoing further not-well understood 
interactions (energy loss, hadronization etc.) However, the enormous background 
due to the many short-lived mesons that decay into photons {p,u!,r], ...) makes ex- 
tricating the signal from the background experimentally very difficult. The most 
interesting result which has come out of this approach so far is CERES 's report 
of observing a dilepton excess (Lenkeit et al. 1999) (see Fig. 1.10 ). This result 
was presented by CERN as part of the evidence for its announcement of having 
found deconfined matter (Heinz and Jacob 2000). However, it is unclear whether 




(1.42) 



47 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

m^^ (GeV/c^) 



Figure 1.10: The di-electron invariant mass (mge) distribution observed at CERES, 
per unit of pseudorapidity (rj) and scaled with the total pseudorapidity density of 
observed charged particles dNJi/dr]. The dotted, dashed and dot-dashed lines rep- 
resent the expected distribution from electromagnetic decays of known resonances. 
An excess is observed at rUee ~ 0.7. 

it signals a dilepton excess due to QGP or, rather, a shift in the p peak due to 
partial chiral symmetry restoration. Measurement of photons/dileptons is part of 
the planned ALICE program, where it can be a very useful probe if correlated 
with jet-quenching (Gale et al. 2004). The energetic reactions in Eq. (1.42) could 
produce a detectable photon/dilepton pair together with either a detectable, par- 
tially quenched, or totally quenched jet. The correlation will then be a powerful 
tomography tool. 
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Figure 1.11: Nuclear modification factor (See sectionl.5.1) of cliarmonium particles 
(detected through the J/'^ ^ fi^fi^ decay) It is plotted against the energy density 
calculated within the Bjorken model (oc R^dN/dy, where Rq, the overlapping nu- 
clear area, depends on centrality) . Minimum bias means that no centrality selection 
was performed. 

1.5.4 Chcirmonium suppression 

One of the original ways to look for deconfinement is through Charmonium suppres- 
sion (Matsui and Satz 1986). Charmonium (bound cc states, such as J/^,^', .-•) 
can be produced in the initial reactions of an intermediate energy heavy ion colli- 
sion. Normally, these particles are very narrow (tightly bound) states. However, 
in a deconfined medium, charmonium pairs will melt due to the color screening 
provided by the free quarks and gluons. Moreover, different bound states will melt 
at different temperatures/densities, thus providing an effective QGP thermometer. 

The NA50 experiment has indeed detected such a suppression (Abreu et al. 
2000) (see Fig. 1.11) which is absent in proton- proton collisions and can not be 
explained through standard nuclear absorption models. This absorption was one of 
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the pieces of evidence CERN put forth in its announcement of quark-gluon plasma 
discovery (Heinz and Jacob 2000). 

However, the NA50 experiment cannot detect total charm, but only J/^ — > 
/i/i reactions. Hence, it is possible to adjust nuclear absorption models through 
in-medium mass modification rather than deconfinement. J/"^ production is par- 
ticularly sensitive to such mass shifts since a small modification of its mass will 
make the J/^ — >^ DD decay possible. In fact, a nuclear absorption model with J/^ 
mass modification does manage to explain the observed J/^' suppression (Chaud- 
huri 2002). The story is not completely over. Other than the fact that the J/^ 
in-medium modification, though reasonable, has never been observed, the nuclear 
absorption model still can not account for the suppression of the 

The NA60 experiment, which will also measure open charm, might shed 
some hght on what exactly causes J/^ absorption at SPS energies. Charm data 
is forthcoming from RHIC and, when it turns on, the LHC (both of which also 
measure open charm). 

As energy increases, we will reach a point where several charm pairs will 
be produced in each coUision. This might change things considerably, since a QGP 
phase will dissolve the existing cc pairs, but will also allow J/^ s to form from 
initially uncorrelated c and c at hadronization (Thews et al. 2001). This means 
that the number of J/"^ will become 

Nj/^^AN,^ + BN^^ (1.43) 

{B — ii only one cc pair is expected per collision). If enough cc s are produced, 
we might observe an enhancement of charmonium produced in heavy ion collisions 
instead of a suppression (Thews et al. 2001). This is especially true if charm and 
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bottom quark states (such as the B^) are also considered, since production of such 
a state requires initially uncorrelated quark pairs to fuse (Schroedter et al. 2000). 



1.5.5 Strangeness enhancement 

Strangeness enhancement has long been considered one of the most promising sig- 
natures of QGP formation (Rafelski and MuUer 1982), as well as a useful tool to 
study soft matter produced in a heavy ion coUision (Rafelski and Letessier 2003a). 
The basic idea is that ss pairs should form much more readily in a QGP than in 
nuclear collisions through reactions such as 



Since gluons are massless and the mass of the light quarks is much less than the 
QGP temperature, the threshold for forming strange quarks is 2ms ~ 100— 300MeV, 
boosted by chiral symmetry restoration and the abundance of qq pairs. In a hadron 
gas, on the other hand, the threshold energy for the leading strange producing 

processes 



is, respectively, 600 and 800 MeV (note that the first process is also suppressed in 
a baryon-poor environment). 

(Rafelski and MuUer 1982; Koch et al. 1986) has shown (Fig. 1.12) that 
in a thermalized gluon-rich perturbative QGP the ss production rate is as much 
as an order of magnitude greater than the strangeness production rate in a hadron 
gas at a similar temperature, with gluon-gluon fusion reactions accounting for a 




(1.44) 



9 




(1.45) 
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large majority of the production. Due to the faster equihbration time and greater 



Figure 1.12: Density of strange quarks with respect to time calculated in a QGP 
(left, normalized to baryon density with ms=150MeV) and a hadron gas (mid- 
dle) kinetic production models. As can be seen, equilibration time is significantly 
shorter for the QGP phase. In addition, as the panel on the right shows, equilib- 
rium strangeness quark content will be stronger in the QGP phase at all chemical 
potentials. 

equilibrium strange quark density, therefore, a system which has undergone a QGP 
phase transition will exhibit an enhancement of strange quarks with respect to a 
system which has not. This enhancement will translate into an even greater en- 
hancement of strange hadrons after hadronization. Multi-strange hadrons (0, H, Q) 
and their anti-particles will be particularly enhanced, since their production in a 
hadron gas has a particularly high energy threshold {pp — > QJl, requiring pro- 
tons with momentum of 700MeV each) or a sequence of many reactions (e.g. 
7rp —>■ KA, ttA —>■ KE, ttS KH.) . In a hadronizing strangeness-rich QGP, on 
the other hand, multi-strange baryons form by recombination of strange quarks, 
and should not be as suppressed according to most hadronization models (Fig. 1.13 



Hence, to test for QGP one has to look for strange particles in a nucleus- 
nucleus collision and compare with p-A or p-p. A clear enhancement, significantly 




left). 
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Figure 1.13: Left: Production of multi-strange hadrons through combination of 
strange quarks from a QGP. Right: Detection of strange particles from a heavy ion 
collision due to decay topology. 

raising with particle strangeness content, would constitute evidence for deconfine- 
ment. Experimentally, this is facilitated by the fact that strange particles decay 
weakly, with a lifetime comparable to the time of flight. Hence, they can be recon- 
structed through decay topology (Fig. 1.13 right). 
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Figure 1.14: The nuclear modification factor (see section 1.5.1) for the production 
of hyperons as measured by the NA57 experiment, normalized by the number of 
participants. Particle yields for p-p and Pb-Pb collisions are compared to p-Be, 
plotted against the number of participants, calculated within the Wounded nucleon 
model (see section 1.5.1) . 
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As Fig. 1.14 shows, a very clear experimental enhancement has indeed been 
detected by the WA97 and NA57 collaborations (Andersen et al. 1999) (whose 
narrow-acceptance telescope is particularly suited to look for multi-strange par- 
ticles). The enhancement in all strange particles, and its rise with strangeness 
content, constituted a crucial part of the supporting evidence for CERN's claim 
to have produced deconfined matter (Heinz and Jacob 2000). Microscopic nuclear 
interaction models have so far been unable to even come close to reproducing the 
observed enhancement (Antinori et al. 1999). 



1000 - 




Figure 1.15: Nuclear modification factor for A (right) and A,S,r2 (left) production 
calculated within a canonical statistical hadronization model. The result is plotted 
against the number of participants Npart (left) and center of mass collision energy 
per nucleon (right). 

One alternative explanation which would account for the observed enhance- 
ment pattern is canonical suppression (Tounsi et al. 2003). Postulating that col- 
lisions in p-p and p-A systems achieve chemical equilibration, the small volume 
of these systems necessitates that the canonical (rather than grand-canonical) en- 
semble be used to calculate strange particle yields. This means that, rather than 
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introducing a chemical potential for strangeness, only outgoing states with exact 
strangeness conservation are counted in the statistical averaging. As Fig. 1.15 shows, 
this leads to a sharp drop in strangeness as system volume decreases below a critical 
value, with the effect being greater as strangeness content increases. 

This explanation, however, has several problems. For a start, it is not at 
all clear that equilibrium statistical mechanics can be used to describe a p-p and 
p-A system. NA57 p-p and p-Pb data are well described through the uRQMD 
nuclear microscopic model (Antinori et al. 1999), where chemical equilibration is 
not assumed or indeed predicted. 

In addition, as Fig. 1.15 (left) shows, the pattern of enhancement is not pre- 
dicted correctly by the canonical picture: the canonical picture has a sharp suppres- 
sion at small volume and a long plateau as volume becomes large and the canonical 
and grand-canonical pictures coincide. Experimentally, however, strangeness scales 
linearly with number of participants in all (p-p,p-Be,p-Pb,Pb-Pb) collisions. The 
scaling in the Pb-Pb system is steeper, but still linear with centrality, in accordance 
with a higher ss production per unit volume per unit time. This is expected from 
QGP-based kinetic production (Rafelski and Letessier 2002). 

A final distinguishing test between canonical suppression and QGP produc- 
tion is to lower the collision energy. As Fig. 1.15 (right) shows, the canonical model 
predicts an increase of enhancement (going to oo as the collision energy approaches 
the S and Q threshold). The QGP model predicts a decrease, with a sharp disconti- 
nuity at the energy where the QGP is not formed anymore. Lower energy SPS runs 
will shortly measure the energy dependence of enhancement and definitely clarify 
the situation. 
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1.5.6 Collectivity 

The argument used in the derivation of strangeness enhancement can also be ap- 
phed to argue that in a quark-gluon plasma system local thermal equilibrium will 
be reached on a much shorter timescale than in a system with hadronic degrees of 
freedom. Hence, the fact that most soft observables can be described through col- 
lective dynamics (thermo and hydrodynamics) is a good indication that QGP has 
formed. Are the particle yields described by a temperature and chemical potential? 
Are particle momentum distributions characterized by one temperature and collec- 
tive flow? How well does a hydrodynamic approach describe the system? A good 
part of the subsequent chapters will be devoted to addressing these topics, so we 
will not dwell on them in detail here beyond mentioning some motivational issues. 

Any system, including a hadronic gas, will evolve collectively if given enough 
time. However, a quark-gluon plasma, with its light colored degrees of freedom, 
should be much more efficient as an entropy generator than a hadron gas. Hence, 
we should expect a much stronger collective signal in a QGP than in a collection of 
hadrons. 

In this respect, we shall mention that microscopic hadronic simulations 
(Bratkovskaya et al. 2004) show that hadronic systems at SPS and RHIC energies 
do not have the time to develop a significant amount of collective flow, and hence 
mass dependence of the apparent temperature (inverse slope of the logarithm of 
the transverse mass distribution) should not be large. This is true for p-p but not 
A- A data (see Fig. 1.16). This is even more true in the case of anisotropic flow 
(discussed in detail in chapter 6). The evolution of anisotropy is described so well 
by the hydrodynamic picture that a very early system thermalization is required 
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Figure 1.16: Mass dependence of the apparent temperature (inverse slope in a 
logarithmic graph) as calculated by uRQMD and HSD models (see section 1.5.1 for 
references) . 

(Heinz and Kolb 2002). It is very difficult to see how a hadronic model could achieve 
such a short thermalization timescale. 

It should also be mentioned that statistical production from an entropy-rich 
QGP is likely to look quite different from freeze-out of an interacting hadron gas. 
In case of a first-order phase transition or sharp cross-over, the system might evolve 
out of equilibrium at the transition, and post-phase transition interactions can po- 
tentially destroy evidence for the earlier equilibration. Given this, it is remarkable 
that several experiments found particles and antiparticles to have exactly the same 
inverse slopes (Bruno et al. 2003) (see Fig. 1.17). In an interacting hadron gas, 
anti-particles should annihilate at a strongly px dependent cross-section. 

Because of the large amount of entropy expected to be generated through 
QGP equilibration, it has also been suggested that to find the phase transition we 
should look for a jump in the entropy per strangeness, entropy per baryon number. 
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Figure 1.17: Transverse mass distributions for particle and anti-particles as mea- 
sured by the WA97 experiment . 

or entropy per energy (Letessier et al. 1992). A recent systematic experimental 
study (Gazdzicki 2004) uses this reasoning to cite the kink observed in vr (~ entropy) 
production at ~ 4GeV (see Fig. 1.18 (Left)) as well as the flattening of the K/ir 
ratio (see Fig. 1.18 (right)) as evidence for deconfinement. 



1.5.7 QGP evolution 

Concluding this overview of signatures, it is necessary to underline one of the most 
remarkable, and problematic aspects of the study of QGP in heavy ion collisions: the 
range of physical approaches which needs to be implemented to study the system. At 
each stage of the system's evolution, the relevant physics and useful approximations 
are in general completely different from that of the other phases. In addition, 
different signatures probe different stages of this evolution. 
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Figure 1.18: Left: tt production per wounded-nucleus a function of en- 
ergy for A-A vs p-p collisions, plotted against Fermi's energy measure 
F = {^snn — IrriNY^'^ / y/s^N^^'^, where y/s^N is the average of the center of mass 
energy per nucleon-nucleon collision. Right: K/vr ratio in nucleus-nucleus collisions 
collisions as a function of yj/vw- 

We shall give a chronological layout of the possible evolution of a QGP 
fireball, together with the methods which can be used to describe it. 



Initial collision The short and energetic initial collision is usually approximated 
through a parton cascade model (Geiger 1993). In this picture, the collid- 
ing partons are described as flat "pancakes" of quarks and gluons, with the 
structure functions extrapolated from nucleus deep inelastic scattering. The 
collision dynamics is then governed by pQCD. 

Recently, an alternative/complementary approach has emerged: That of the 
Color Glass Condensate (lancu et al. 2002). From saturation physics one can 
argue that at high energy heavy ion collisions the gluon occupation number 
is significantly larger than one. Hence, the colliding system can be modeled 
as an incoherent classical SU (3) field. 
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Hydrodynamic evolution Somehow (the details are not entirely clear), the ini- 
tial system reaches local thermal equilibrium in a (hopefully!) deconfined 
phase. In this case, thermahzation times will be very fast (Rafelski and MuUer 
1982) and hence the mean free path very small. The bulk evolution of the 
system can then be optimally described by relativistic hydrodynamics. 

As we shall see later in the thesis, however, thermal and chemical equilibrium 
can be two very different things. Chemical equilibration is therefore not nec- 
essarily as fast, and in fact does not necessarily occur for all quantities. This 
aspect of QGP evolution is best modeled by transport theory, with the colli- 
sion terms given by quantum field theory (Rafelski and MuUer 1982),(Kapusta 
et al. 1991), (Wang et al. 1995). 

Hadronization This is the big unknown. Quite simply, there is no rigorous phys- 
ical way to approach it. It is non-perturbative, so Feynman diagrams will not 
work. It is also likely to be a far from equilibrium process, so it is doubtful 
that lattice QCD will work as more than a qualitative indication of what's 
going on. 

All we are left are effective models, based on basic physics such as thermo- 
dynamics, relativity and conservation laws. One way to make quantitative 
predictions is entropy maximization. This leads to statistical hadronization, 
the subject of this thesis. Another recently proposed ansatz is the coales- 
cence/recombination of partons (Pries et al. 2003). This picture has been 
proven to explain a range of phenomena at intermediate pr- However, it can- 
not as yet be formulated self-consistently in a soft regime, where particle mass 
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is non-negligible, as it violates energy-momentum conservation (as well as en- 
tropy non- decrease). Perhaps, in the future, it will be possible to combine 
this approach with statistical hadronization (Biro and MuUer 2004). 

It should be noted that some observables described in this section do not have 
to undergo hadronization. Photons and dileptons leave the QGP with no 
further interactions. The impact of hadronization on jets is also dynamically 
not relevant. Since hadronization is a soft process, it will not change the 
momentum distribution of the jet significantly, but will simply "dress" the 
outgoing quark in a hadronic coating. 

Soft probes (strangeness, charm, collectivity), on the other hand, will undergo 
a potentially non-trivial hadronization before reaching the observer. This 
means understanding hadronization is a crucial step in understanding soft 
physics. It also means, however, that soft probes are tools in understanding 
the phase transition from a quark-gluon plasma to a hadron gas. 

Post- hadronization evolution It is far from obvious that hadrons stop interact- 
ing after they form. In fact, one of the most important questions one can ask 
about hadronization is the extent of interactions which follow it. 

One thing which is certain is that hadronic phase space is very different from 
quark-gluon plasma phase space. Hence, an equilibrated quark-gluon plasma 
will probably produce an out-of-equilibrium hadron gas. Equilibration times 
for a hadron gas are much slower than for a QGP, since the hghtest color- 
neutral boson has a mass of 137 MeV. Hence, the best way to analyze such 
a system is through transport theory. Collision terms are measured experi- 
mentally in elementary coUisions or calculated through effective field theory 
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techniques. Quantum molecular dynamics models (Bass et al. 1998; Ehehalt 
and Gassing 1996) are based on this approach. 

Final freeze-out At this stage, particles decouple from the system and reach the 
observer. It is sometimes fashionable to talk about a "chemical freeze-out" 
(when inelastic collisions cease) and "thermal freeze-out" (when particles stop 
interacting altogether). A look at particle reaction data shows that this de- 
scription is inappropriate: Inelastic reactions, such as strange quark exchange, 
can be just as low-threshold as elastic hadronic reactions. It is better to talk 
about hadronization (when hadrons become appropriate degrees of freedom) 
and freeze-out. 

1.6 The scope of this thesis 

This thesis, therefore, explores the effect of statistical hadronization on the observed 
soft particle abundances and spectra in heavy ion collision. We shall concentrate on 
hyperons, given their relevance to strangeness enhancement. However, non-strange 
hadrons shall also be examined. 

We shall use experimental data (from SPS and RHIC) to constrain the 
statistical hadronization picture, and differentiate between freeze-out scenarios. We 
will test these scenario's ability to constrain both hadron yields and spectra, as well 
as direct detection of unstable short-lived resonances. 

Some of the questions which we would like to ask are: 
• What are the temperature and chemical potentials characterizing the system? 
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• Is there evidence of transverse flow and other coUective effects? 

• Does statistical hadronization happen in chemical equilibrium, or are non- 
equilibrium effects important? 

• What is the dynamics of statistical hadronization? How does particle emission 
proceed in spacetime? 

• What is the significance of the post-hadronization interacting hadron gas 
phase? What is its impact on observables? 

• What is the experimental data's sensitivity to these observables? What is the 
significance of fit results, and how can they be further confirmed? 

In the second chapter, we shall review relativistic statistical mechanics and how it 
applies to hadronization. We shall also describe how a variety of different freeze- 
out scenarios arise out of statistical hadronization. In chapters 3 and 4 we will use 
statistical hadronization to describe particle spectra measured, respectively, in SPS 
and RHIC experiments. In chapter 5, we will use yields of short-lived resonances 
as a way to differentiate between the different freeze-out models and to constrain 
freeze-out dynamics. In chapter 6 we shall use azimuthal asymmetry in less than 
central collisions as a tool to study freeze-out. Finally, we will give a few conclusions 
and an outlook on unresolved issues. 
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CHAPTER 2 



Statistical hadronization: an overview 



2.1 Introduction 

The use of statistical mechanics to describe production of strongly interacting parti- 
cles was pioneered by Fermi (Fermi 1950) and developed by Pomeranchuk (Pomer- 
anchuk 1951), Landau (Landau 1953) and Hagedorn (Hagedorn 1965). Indeed, 
Hagedorn's observation (Hagedorn 1965) that an infinite number of hadrons of in- 
creasing mass leads to an exponential density of states, and hence to a critical 
temperature where the canonical partition function diverges, provided a piece of 
evidence for a QCD phase transition before even the discovery of quarks. 

The basic idea of this approach is that in a high-energy process driven by 
strong interaction a large number of particles is hkely to be produced (Fermi 1950). 
In these circumstances, the dynamic part of the reaction cross-section will average 
out, and be considered a "volume" constant. The distribution of each particle 
(1...N) will then be given by the corresponding phase space weight ^ 



It can be shown (Fermi 1950; Hagedorn 1965; Rischke 2002) that in the many-body 

infinite energy limit this distribution approaches an exponential, with the slope 
^See appendix C 




(2.1) 
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related to the energy per particle (canonical limit). If the phase space in Eq. (2.1) 
also provides for quantum number conservation, it becomes possible to derive a 
grand-canonical limit. 

The statistical model has been applied to every experimentally studied 



2003). It is not obvious, however, that particles described by such a statistical model 
are effectively thermalized. In particular, a quark-gluon plasma, should exhibit a 
further degree of equilibration: The fact that its composed of colored, massless 
degrees of freedom means it should achieve local thermalization in a timescale con- 
siderably shorter than its evolution. It is therefore expected that the QGP evolves 
as a continuous fluid. 

2.2 Thermodynamics and hydrodynamics from Kinetic theory 

Thermodynamics and hydrodynamics also arise as a r ^ oo limit to the Boltzmann 

transport equation governing many-particle distributions 



Where C[f] is the Entropy-generating particle collision term. We can describe the 
collision term semi-classically, by taking the cross-sections calculated in quantum 
field theory a{P,P' — > p,p') but assuming that the distributions of the particles 

are uncorrelatcd beyond Quantum statistics requirements. For a theory invariant 



strongly interacting system, from e^e 



hadrons to heavy ions (Becattini et al. 




(2.2) 
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under time reversal, and considering just two body interactions^ 

c[f]{x,p)^ J d'[x\X,X',p,P,P']S\P + P'-p-p')a{P,P'^p,p') (2.3) 

{[f{X,P)fiX',P') + FfiX,P')fiX',P)] - [fix,p)f{x',p') + Ffix,p')f{x',p)]} 
where F ensures the Fermion-Boson anti-symmetrization requirement 

— 1 Fermions 

1 Bosons (2.4) 
Boltzmann{Distinguishable) 

The Boltzmann H-theorem, generahzed to quantum statistics (Lifshitz and 
Pitaevskii 1981) states that Entropy 

S^j d'[x', X, X',p, P, P'] {-f(x,p) Hf(x,p)] - F{1 + Ff(x,p)) ln[l + Ff{x,p)]) 

(2.5) 

will always monotonically increase, until detailed balance is reached, C[f] will go do 
and f reduces to the entropy-maximizing Fermi-Dirac, Bose Einstein or Maxwell- 
Boltzmann distribution 

^ oo 

W, x) = rnn^ — = ^"-^"e""''^'''/^ (2.6) 

1 n=l 

This will happen after "many" collisions will have taken place 

^oo » ^^^^N^ (2-7) 

< (7 >< y > 

(too is actually reached quite quickly in a strongly interacting system with high 

particle density) 

■^In general, the frequency of many-body interactions goes as a power of the mean-free-path, 
< cr >niady (< N /V >)~^'""^^ . Hence, C[f] is actually a perturbative series and Eq. (2.3) is 
the first term. Both a QGP (Rafclski and Mullcr 1982) and a hadron gas (Bass ct al. 1998) 
have been described using just the two-body term. However, this approach will break down in 
the non-perturbative limit, or equivalently when long-range correlations become too large for the 
semi-classical approximation. A QGP hadronization can not therefore be described this way. 
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/o in Eq. (2.6) is characterized by temperature T and fugacity A, as well 
as u^{x) = 7(1, v), the 4- velocity vector representing the motion of the volume 
element. If the interactions will be strong enough that the mean free path of the 
particles is negligible with respect to the collective motion motion of the system 

» — (2.8 

dT Too 

the system will always have local thermal equilibrium and evolve as a continuous 
fluid 

9^7^^^ = (2.9) 
d^n^" = (2.10) 

where T^" is the energy-momentum tensor and d^^n^ is the number current (for any 
conserved current such as isospin, strangeness, etc.) We can express these quantities 
in terms of pressure (P), energy density (e) and number density 

^ J ^/o(pM^ ^) ^ p^M- + (p + e)u''u'' (2.11) 

n^ = J n^fo{p^,x) (2.12) 

Hence, if we know the equation of state of our fluid we will be able to close our 
system of equations and evolve the system from any set of initial conditions. 

2.3 Statistical hadronization and the Cooper-Frye formula 

QGP, with its high number density of strongly interacting particles, is a good can- 
didate for an equilibrated fluid. However, this situation will change in the QGP-HG 
transition, in which the 10 colored massless partons become over 200 color-neutral 
hadrons. 
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The Boltzmann formalism of Eq. (2.2) can not, as far as we can see, be 
applied to hadronization, since the assumptions that go in it (weakly correlated 
particles scattering locally) do not apply to a non-perturbative quantum phase 
transition. However, we know from lattice QCD that hadronization is a relatively 
fast process, and obviously hadronization of a large thermally equilibrated system 
can not decrease entropy. For this reason, and the fact that color-neutral hadrons are 
more weakly interacting than quarks, it is very likely that an equilibrated entropy- 
rich QGP will emit hadrons according to phase space as described in section 2.1. 

If the system is large, energy and quantum number conservation can be 
accomplished by Lagrange multipliers, so, in the rest-frame with respect to the 
collective flow the hadrons will be distributed according to /o as given in Eq. (2.6). 

the post hadronization hadron current will be given by 



where the Hadron fugacity is given by the product of the fugacities of its constituent 
quarks quantum (not necessarily the QGP's chemical potentials, as we'll see later) 



Using the "fast hadronization" assumption again, we can define a "Preeze-out hy- 
persurface" defining a locus in space-time from which these statistically hadronizing 
particles are emitted, labeled by a 4- vector E'*. Since T,^ is a 3-surface in Minkowski 
4-space, it can always be expressed as a function of three parameters {u,v,w). Its 
element can then be given in in a Lorentz-covariant way using Stokes's theorem 




(2.13) 



hadron 



(2.14) 



diy aE'* aE^ 

du dv dw 



(2.15) 
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where e^,^Q,^ is the Levi-Civita symbol, the number of particles produced in such a 
volume element is then Lorentz-invariant, and computable in the volume element's 
rest frame as 

ji^d^j:" = dN (2.16) 

combining Eq. (2.13) and Eq. (2.16) we obtain the famous Lorentz-invariant Cooper- 
Prye formula (Cooper and Prye 1974) 

= I d'j:^p^f(p,u^,T,X) (2.17) 

If hadronization of the full volume takes a lot of time, the emitted particle can 
find itself in the QGP again if p^T^^ < 0. It is still unclear how to handle this 
in a rigorous way, but typically one just truncates the distribution to exclude this 
unphysical region (perhaps adjusting the limits of integration to conserve energy, 
entropy etc.) (Bugaev 1996),(Anderlik et al. 1999) 

= j d'^,p^f{pX, T, A)e(S^p^) (2.18) 

where Q{x) is the step function. 

Since all hadrons apart from the pion are considerably heavier than the 
typical hadronization temperature {p^u^ > m >> T), the sum in Eq. (2.6) can be 
truncated at n=0, corresponding to the Boltzmann approximation 

f{p,) = Ae-^"'^"/^ (2.19) 

2.3.1 Total particle yields 

If the thermodynamic parameters T, A do not vary within the hadronizing volume, 
the total number of particles will be independent of E'* and u^. To see this, we 
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integrate the Cooper-Prye formula over momentum space 

JdN^ J ^p^^d'E^XeP^-'/^ (2.20) 
and insert a U/^u^ {— 1) in the integrand. The two integrals then decouple 



N 



(2.21) 



The first integral is just a normalization constant. Since p^u^ — Erest, we are left 
with 

N^Vg J d^pAe-V^^/^ (2.22) 

where g is the particle degeneracy (For a colorless particle, it will be equal to 25*+ 1 
where S is the spin) ^. Using the Bessel function definition (Gradshteyn and Ryzhik 
1989) 

(2.23) 



we can find N analytically 



(2n)!" 



X^ 



An 



WXTK, 



(2.24) 



The energy density in a rest frame with respect to the collective flow also follows 

. = / ^'V^'Xe-V^'^ ^ A, (^) . A, (^)) 

(2.25) 

Comparing Eq. (2.19) and Eq. (2.6) it can be seen that these formulae can be 
generalized to the Fermi- Dirac and Bose- Einstein distribution by the substitution 

^™'(f)-E(±i)"*'^^"^'(t^) 

n=l 

■^It is perhaps not immediately intuitive how a non interacting gas dispersion relation can be 
used to describe a strongly interacting system. The "trick" , here, is that we will generate all the 
strongly excited states (resonances) in the data book, and the strong interaction excited states can 
be considered, to a good approximation, as narrow resonances. As was explicitly shown (Dashen 
and Rajaraman 1974) if resonances are narrow enough, an ideal gas of resonances is an appropriate 
description. As this section shows, our approximation can be improved through finite width eflfects 



70 

If the particle has a finite width, the Bessel function will be further in- 
tegrated over the range of masses to take the mass spread into account. For a 
resonance with width, we obtain 

^iK) ^TfJ2 ni{M)Fr{M, r,(6,_,-, M))dM (2.27) 



where 



• bi^j is the decays branching ratio 

• TTT'threshoid IS the threshold mass for the decay, rrij 

• Fr(M, Fi) = ^M-m^^+r'^/A Breit-Wigner formula 

• ri{bi-^j, M) is the energy-dependent width for the decay under consideration. 
In general, this is a non-trivial particle dependent function. Here, we only 
consider the dominant energy dependence of the width, namely the decay 
threshold energy phase space factor. The explicit form for decays with low 
relative angular momentum has been studied through corresponding reverse 
production cross-sections, and found to be (Terazawa 1995; Gounaris and 
Sakurai 1968; Pratt and Bauer 2003) 

(2.28) 

where I is the relative angular momentum of the decay and F* is the energy- 
independent constant found in the particle data book (Hagiwara et al. 2002) 

A^o is the Breit-Wigner and phase space normalization 

poo 

N= Fhre^t-u,^aneriM,Ti{b^^J,M))dM (2.29) 

'mthreshold 



F,(6,_,-,M)=6, 



^ _ f ^threshold \ ' 
\ m / 
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f^i /-i imT (b ^) 

= / F,^eit-wigner{M,Ti{bi^j,M))dM + / tl^=l^^J-—^dz 



2.4 Treatment of resonance decays 



Most of the 200 particles produced in statistical hadronization will be short-lived 
resonances, whose decay is undetectable through an analysis of the particle trajec- 
tories as it happens after the particle was produced. More complicated cascades, 
with sequential decays, are also possible. While the number of particles made from 
these decays is suppressed by the high mass of the resonance, they are also enhanced 
by resonance's typically high spin degeneracy and the fact that many particles are 
produced in a typical decay (Fig. (2.1) left). 

The expected yield of observed "hght" particles will then have to include 
the products of all resonance decays of short-lived products. As Fig. 2.1 (right) 
shows, there is no way that this contribution can be neglected. However, on the 
bright side, resonances are an important tool for QGP diagnostics: Since their quark 
numbers are the same as the base particles, their relative abundance is controlled 
by the temperature only. Their abundance is also sensitive to post-hadronization 
dynamics (see chapter 5). In short, resonances will be a very sensitive freeze-out 
probe, one which will be explored from different angles in the next chapters of this 
thesis. 

^This change of variables, ^ = z, used within our programs, is absolutely essential for the 
numerical evaluation of the integral. Otherwise, the result will be plagued by F-dependent system- 
atic errors, since N converges slowly as the upper limit of the numerical integration oo. On the 
other hand, the integral with changed variables is also prone to error if the width is too small, since 
the Gaussian integration might miss the peak. The solution, which our calculations implement 
(Torrieri et al. 2004), is to split the integral in two parts: f°° = f'"°+^^ -|- f°° The 

first part can be done through conventional integration, the second through variable change 
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Figure 2.1: Left: Typical resonance decay patterns. Right: Resonance contribution 
to observed particle yield as a function of freeze-out temperature . 

To include resonances in particle yields and ratios, it is enough to sum the 
resonance decay products to the base yields 



Vi^?;+... 



(2.30) 



where Ni is given as in Eq. (2.24) and bj^i^^^^ is the branching ratio. The only 
technical difficulty in this calculation is that decay products in the particle data 
book (Hagiwara et al. 2002) are given up to Clebsh-Gordan coefficient factor of 
Isospin. In case of two-body decays, we calculated C-G coefficients automatically 
using the algorithms in CERNLIB (Various 1999). Three body decays, are a bit 
more problematic, since there are several possible ways of combining partial isospin 
sums. We averaged over these, giving ( J, m a and refer to the total and z 
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isospin components of the resonance and decay products respectively) 

bj^i2 = ^'o-.i2(< Jomo\jij2mim2 >f (2.31) 

1 X 

bj^l23 = ^0^123 3E(< Ji2mi2\jij2mim2 >< Jo"^o|ii2i3^"i2"^3 >)^ (2.32) 

312 

+ X](< Jismisljij^mims >< Jo"^o|ii3j2"^i3"^2 >)^ 

jl3 

+ 5Z(< J23m23\j2j3m2m3 >< Jomo\j23jim23mi >f] 

323 

2.4.1 Spectra 

The contribution of resonances to particle spectra is non-trivial, since the full kine- 
matics of the decay has to be understood. However, the resonance contribution 
to spectra is where the data is really sensitive to freeze-out dynamics since post- 
hadronization interactions might re-thermalize resonance decay products. Hence, 
searching for a bona-fide resonance contribution which could be reconstructed by 
invariant mass is a stringent test of fast statistical freeze-out. 

In analogy with Eq. (2.30) we assume that 

where the first term is given by the Cooper-Frye formula and the second gives the 
distribution of decay products in a cascade. Assuming j is also produced through 
statistical hadronization, all is left is to find -E'l^^ in terms of Ej^^ (which will 
again be given by the Cooper-Prye formula. The relations we will develop will 
themselves be recursive, to handle the case when j has a cascade component). 

We shall assume that all the short-lived resonances have decayed. We shall 
also assume that lots of particles are emitted in all directions, thereby averaging 
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over any matrix-element dependence on factors such as spin. The distribution of 
the resonance decay product momenta will then be given simply by the integral 
over the available region of the Lorentz- invariant phase space ^. 

-g^ tJ E,-—L (2.34) 



n ,^ N N 

= n ^^^T.P^ - P'^^^H E^-E,- E,) (2.35) 

1=2 * i=2 1=2 

B is a normalization factor to make sure that 

= h^.2..n (2.36) 



Two body decay 

For a two-body decay, j — > 12, we can go to the resonance's rest frame (denoted by 
the *) where 

El^^{ml + m]-ml) (2.37) 



p\ = ^El-m\ = -p\ (2.38) 
and take advantage of energy conservation and Lorentz-invariance 

=P>r=^>,- (2.39) 
hence, the integral in Eq. (2.34) for just two bodies reduces to 

After expanding the delta function argument explicitly we get a condition in terms 
of the transverse mass and momentum {jnTiVT) and rapidity {xff 

rriTimTj cosh(|/i - yj) - ptiPtj cos(6^i - 9j) 



rrij 



El = (2.41) 



^See appendix C 
^See appendix A 
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Rearranging this becomes 

(eW-2/i)2 _ 2Ae^.-!/i + 1 = (2.42) 

^ ^ mjEl+p-riPTjCosiOi - gj) 
niTimTj 

Which introduces constraints hnking rapidity, transverse mass and relative angle 

Vj = yi + ln(^ ± VA'^ - 1) A > 1 (2.44) 

Eq. (2.44) can be used to eliminate rapidity from the integration and to introduce 

integration constraints on uit'- 

mjEl+pTiPTjCos{9i-9j) ^ ^ 
rriTimTj 

together with | cos(^i — 9j)\ < 1 leads to 

rrijEi — pTiPTj > rriTimTj (2-46) 

with the limiting cases 



rrijEl —pTi^m^j - = rriTimTj (2.47) 
Rearranging and squaring, we get 

PTi(^Ti ~ '"^j) = {'rrijEi — iriTimTjY (2.48) 
putting this in the quadratic form and remembering that m^-^ — pf.^ = m\ 

mlrri^j - 2mTjmTimjEl + m]{p%^ + Ef) = (2.49) 

Using the quadratic formula, we get two solutions (a factor of 2 cancels out from 
top and bottom) 

mTirrijEl ± w {rriTimjy — m\m?AEl -\- p^i) 
TTlTj = ^ 5 (2.50) 
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which simphfy since 



= nij y p\^Ef - p\^ml = rrijPTiPl (2.51) 
If one puts this into Eq. (2.50) one gets 
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= -^{rriTiEt ± p^ipD (2.52) 
reassuringly, this is always physical since ttitiEI > priP* 

In a similar vein, Eq. (2.44) yields the constraints for the relative angle 9j. 
Putting in the form of A (Eq. 2.43) into Eq. (2.44) and manipulating one gets 

-^^•^^*+"^^^^^<cos(^,-^0 (2.53) 



PTlPTj 

f —rriiE* + rriTimT^ ^ 



- arccos ( ZlIhEJ^^I}!^] < < 6, + arccos f -^^El + mr^rurA 

PTlPTj J V PTlPTj 



Putting everything together, we finally have an analytically solvable expression for 
the decay products distribution in terms of the resonance distribution 

= B j mTjdmTj J d9, J dy,5{y, = 2/1 + \n{A ± VI^)) [E,^] 

(2.55) 

^ _ mjEl+ PTlPTj cos{ei-ej) 
mTiiriTj 

e^ = e,± arccos ( -^^^l + ^TirurA 

\ PTlPTj ) 

Integrating over Eq. (2.56) to find Ni/Nj we discover that the normalization has to 
be 

(2.58) 
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An important simplification arises in the rapidity- invariant limit (Bjorken 1983): 
In this case, the integral over the S function in Eq. (2.55) has no effect, and the 
different Fourier components of the distributions 

= vio + vu cos{9) + vu cos{29)... (2.59) 

dN- 

Ej-z^ = vjo + Vji cos{9) + Vj2 cos(2^)... (2.60) 
a pj 

do not mix (ie vjo has no effect on vim^i) (Baran et al. 2003) 

If the original distribution function is azimuthally invariant (no 6 depen- 
dence), Eq. (2.55) simphfies to (SoUfrank et al. 1990) 

^ „ ^ = - — - / dy. / dmLj ^ „ ^ , 2.61 
dm^^dyi ^npt Jy_ Jmr- dm^jdyj 

J = , =, (2.62) 

JptjPti - (iT^REt - mTjiriTi cosh{yj - yi)^) 

y± = y,±smh-^(^] (2.63) 



ElrriTi cosh(|/j - yi) ± pri \ pf - m^Ti sinh {yj - yi) 

= rnj 2 • (2-64) 

mj.^ smh [yj — yi) + mf 



Decays of three bodies and more 

The strategy to pursue if there are more than 3 bodies in the decay is to change 
variables from pi,p2, ...,Pn to pi,P2, ■■■,Pn-i, S2..n,P2..n where 

n n 

S2..n^{Y.^^^'-^Y.P^y (2-65) 

1=2 i=2 

is the invariant mass of "all the other" decay products (Byckling and Kajantie 1969) 
and p2..n is their combined momentum. The invariant phase-space in Eq. (2.34) 
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becomes, applying the 5 function in Eq. (2.40) 

(9(pi,P2, ■■■,Pn-l,S2..n:P2..n) 



El-rr- = B / dS2..n 



d^Pj X f PmPi 



d{pi,...,Pr, 



(2.66) 



rrii 



El / K 



' d^Pj 



El.n = sjpl.n + S2.. 



Ei^^pi.n + mi 



here T„ is given as in Eq. (2.34), the Jacobian is 



d{pi,P2, ■■■,Pn-l,S2..n,P2..n) 



d{pi: ...,Pn) 



^Jruj - (g2..n + m^f^jruj - (g2..n - "m^)' 



(2.67) 



(2.68) 



and the hmits of integration are given by energy conservation 

n 

s- = ^f s+ = (m,-mi)2 (2.69) 

1=2 

Thus, the n-body decay can be expressed as a convolution of a 2-body and an n-1 
body decay, and further convoluted to n 2-body decays. 

Unsurprisingly, the decay products of the cascades such as 

r^' — > TyTT — > TTTTTTTr (2.70) 

(Fig. 2.1 left) will have a very similar distribution: Ej^^ of the intermediate state 
(77 in Eq. (2.70) ) will be fed back into a formula such as Eq. (2.34), and the only 
difference from Eq. (2.66) is that the intermediate mass is a constant 

d{pi,P2,...,Pn~l,S2..n,P2..n 



d{Pu...:Pn) 



S(^S ^^intermediate) 



(2.71) 



or a Lorentzian if F 7^ 0. 



If the resonance distribution is azimuthally invariant Eq. (2.66) for the 3- 
body case simplifies to (SoUfrank et al. 1990) 



d^Pi 



-D I -C/s^23 



d^p. 



s-»23 



(2.72) 
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B 



2'KQ[{mj + mi)2, {m.j - mi)^, {m.2 - m^Y, {m.2 + ma)^] 



(2.73) 



Q[a, b, c, d\ 




— \/{a — x){b — x){x — c){x — d) 



(2.74) 



where ( E. 




) 



is given by Eq. (2.61). 



Prom the above discussion, it is clear that in general an n-body decay will 
involve 2n + l numerical integrals. Cascades of n steps with a rii decay at each step 



resonance distribution. 

It is clear that in general Monte-Carlo integration becomes the only way to 
evolve decays beyond a simple topology. Algorithms such as MAMBO (Kleiss and 
Stirling 1992) can generate points with a uniform density in n-body phase space. 
In the subsequent chapters we will use MAMBO to perform integrals within the 
Monte-Carlo program, and multi-dimensional Gaussian integration (Various 1999) 
in our fits. The latter, therefore, can only accommodate a restricted number of 
resonances where the decay topology is not too complicated. 

2.5 What kind of statistical hadronization? 

The formation of QGP should be accompanied by statistical hadronization at some 
point: This, in a sense, true by definition, since phase transitions characterize 
equilibrated systems. However, by itself, statistical hadronization is not a proof of 
quark gluon plasma formation, since any system will thermalize after enough time 
as discussed in section 2.2. 

While, as discussed in sections 1.5.5 and 1.5.6, the equilibration time of a 



will require a total of Yl^=i integrals to find the final decay products from the 
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system of quarks should be larger by an order of magnitude than that of a system 
of hadrons, hence evidence of early thermalization can be construed as a QGP 
signature, some thought should be put into whether statistical hadronization of a 
QGP can be distinguished from freeze-out of an equilibrated hadronic system. 

An obvious signal to look for is whether the hadron formation temperature 
approaches the temperature predicted by lattice for deconfinement. However, as 
explained in chapter 1, lattice results have still to converge on a definite value for the 
phase transition temperature. In addition, if the phase transition is sharp enough, 
phenomena such as super-cooling can lower the temperature at which hadrons are 
emitted. 

A phase transition from a QGP to a hadron gas should also considerably 
increase the system's mean free path. Due to this, the Cooper-Frye prescription 
should be particularly appropriate. While this has certain consequences, which the 
next two chapters will explore, it is also hardly a definite distinction, due to the 
Cooper-Frye formula's fiexibility in the choice of d^T, and u'^ (see chapter 4) 

Finally, it should be mentioned that the phase space available to QGP is 
very different to that available to a HG: The hghtest degree of freedom in a (pertur- 
bative) QGP is massless, with a degeneracy of 16. The hghtest degree of freedom 
in a Hadron Gas, on the other hand, has a mass comparable to the hadronization 
temperature and a degeneracy of 3. This suggests that, if the phase transition is 
sharp, statistical emission of hadrons from a QGP will not be based on chemical 
equilibration at all, but rather on production of hadrons from quarks with thermal 
weights. Section 2.4.3 describes a statistical hadronization model based on these 
assumptions. 
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These considerations suggest that "statistical hadronization" actually en- 
compasses a variety of models, different in both physical interpretation and ob- 
servational consequences. This is indeed the case. The following two sub-sections 
describe the two "camps" into which statistical models can be categorized. 

2.5.1 Long freeze-out 

The most naive model for freeze-out assumes particle formation at the deconfine- 
ment temperature, a long-interacting hadron gas phase, and thermal freeze-out at 
~ 100 MeV. Superficially, such a picture has experimental support. It is possible 
to fit hadrons produced at RHIC (Braun-Munzinger et al. 2001) and SPS (Braun- 
Munzinger et al. 1999) (Fig. (2.2)) using a model based on chemical equilibration 
Xq = and Eqns. (2.24) and (2.30) , and get a temperature of the order of 
~ 170 MeV. The chemical potential decreases with higher energy, as expected (see 
Chapter 1). In this model hadrons will continue interacting, so particle spectra will 
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Figure 2.2: Equilibration fits for SPS (left) and RHIC (right) energies . The 
xVDoF = 38/16 at the SPS and 5.7/7 at RHIC. 

have to be modeled with a different set of thermodynamic parameters. Experiments 
have done these fits, neglecting resonance abundances, and have found a spectra 
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freeze-out temperature of lOOMeV for all particles except S s and Q, consistent 
with a staged thermal freeze-out (Adams et al. 2003b), (Afanasiev et al. 2003) 
where particles with fewer interactions with pions (particularly the as there is 
no Qtt resonance) freeze-out earlier. 

T,^ in the Cooper-Prye formula will then not be a good physical description 
of particle emission, since particles will not be emitted from a definite point in 
spacetime, but will emerge from a "freeze-out hyper-volume" (a continuous limit 
of fewer and fewer interactions). However, the Cooper-Prye description can still be 
used as a prescription to switch from a continuous medium to a gas of hadrons (not 
necessarily at a phase transition, but simply the space-time region the mean free 
path becomes "large") (Teaney et al. 2001) 

This picture, however, has problems both at a fundamental level and 
in terms of experimental agreement. The fits in Fig. 2.2 actually have a very 
small statistical significance. The experimental fits performed in (Adams et al. 
2003b), (Afanasiev et al. 2003) are demonstrably meaningless, since they do not 
include resonances which have already been experimentally detected (see chapter 
5). On a fundamental level, It is hard to see how entropy gets conserved in a QGP- 
HG chemically equilibrated phase transition unless the transition is slow enough to 
permit the system's volume to increase significantly (something not corroborated 
by the lattice). 

It is also difficult to see why particle spectra should be described by a 
hydrodynamical model at all, given that an interacting hadron gas at T = 120 — 170 
should include many out-of-equilibrium inelastic interactions. In particular, the 
possibility for baryon-antibaryon annihilation at any cross-section makes it puzzling 
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that the uit slopes for particle and anti-particle as virtually identical, as shown 
in chapter one. (Fig. 1.16, (Bruno et al. 2003)). This feature is predicted by 
hydrodynamic models with a few simple reactions (quark- level kinetics, gg <-> qq), 
but certainly not in an interacting hadrons gas picture with many particle-specific 
channels (A -|- A pions,p -|- A Pions and Kaons,...) . 

Another failure of the low temperature scenario is Hanbury-Brown-Twiss 
interferometry (HBT) applied to tt and K. This technique, which relies on bosonic 
two-particle correlations assuming no quantum mechanics before emission, has been 
used to estimate the spacetime shape of the system when interactions stop (Wiede- 
mann and Heinz 1999). It has been found that "naive" hydrodynamics with early 
freeze-out fails to explain the observed correlations (Heinz and Kolb 2002) unless 
fast emission and a small freeze-out radius (which contradict the model) are intro- 
duced by hand. 

While the effectiveness of HBT as a description has been questioned in light 
of its assumptions (For instance, how to properly take short lived resonances (Bolz 
et al. 1993) and final state interactions (Gastineau and Aichelin 2000; Wong 2004) 
into account?), the need of hydrodynamic models to introduce fast freeze-out Ad 
Hoc to even qualitatively explain the data, together with its other failings, is a 
strong motivation to look for a scenario where fast freeze-out is part of the model, 
and particles stop interacting soon after formation. 

2.5.2 Explosive freeze-out 

One framework which would explain the rapid freeze-out and lack of post- 
hadronization interactions is explosive hadronization of a super-cooled plasma 
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(Csorgo and Csernai 1994),(Rafelski and Letessier 2000). If the QGP-HG phase 
transition is first order, at a certain critical temperature/density the two phases will 
coexist. We can estimate this critical density through the Bag model, described in 
the previous chapter (section 1.4.1) 



(2.75) 



where B is the vacuum pressure/bag constant. If the Quark gluon plasma exhibits 
strong collective expansion, ie 



(2.76) 



it will expand past the phase of coexistence. The QGP will then experience negative 
pressure, from the outside vacuum. What happens in this situation is described 
in hydrodynamics as "viscous fingering" (Sonin 1991): A mechanical instabihty 
develops on the surface where the two phases meet, and the lower pressure QGP 
"fingers" into the higher-pressure vacuum. Hadron emission will then occur through 
"bubbhng" at that surface, as the mechanical instabihty tears the fireball apart 
(Fig. 2.3). The hadronization freeze-out surface will then acquire a physical meaning 
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V 



Flow 




Hadrons 
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Figure 2.3: How expansion plus vacuum pressure leads to viscous fingering . 
as the surface where the two vacua come into contact and viscous fingers develop. 
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Indeed, the whole freeze-out dynamics is modeled through the choice of T,^ and the 
Matching conditions (Csernai et al. 2003) 

d'^,{T^GP - T^g) = (2.77) 
d'j:,(n^Qap - <g) = (2.78) 
d'^.iS^HG - S'qgp) > (2.79) 

Where n'^ refers to any conserved quantum number current (flavor, strangeness, 
charge etc.). It becomes clear that the last equation is the most problematic. While 
the first two equations can generally be satisfied through a readjustment of the 
thermodynamic parameters at freeze-out (Keranen et al. 2003), the entropy of the 
hadron gas will in general tend to always be much lower than that of the QGP due 
to the fact that HG degrees of freedom are much more massive. In a slow freeze-out 
scenario entropy will be generated through flow. If freeze-out is sudden, however, 
entropy conservation across the freeze-out surface in equilibrium is only possible by 
either choosing carefully the freeze-out parameters or post-hadronization reheating 
(release of latent heat) (Keranen et al. 2003). 

A more general freeze-out scenario is provided by dropping the assumption 
that hadronization happens under flavor chemical equilibrium. Chemical equilib- 
rium requires that particle and antiparticle fugacities are inverse of each other 

A, = A^^ (2.80) 

which, by the law of mass action, corresponds to the equilibration of creation and 
annihilation processes 

qq ^ 99 ^ 99 ^ qq (2.8I) 

According to the Boltzmann equation Eq. (2.2) this condition will be reached for 
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all quark flavors after infinite time. To model the approach of this condition, one 

defines the phase space occupancy 7^ so that 

Ag ^ Ag7g 7g = 7g (2.82) 

Equilibrium will then be reached when 7^ = 1. Several authors have allowed the 
possibility that strangeness is not yet in chemical equilibration when freeze-out 
happens (Becattini et al. 2003), but it is still generally assumed that fight quarks 
are equilibrated. It is not immediately clear, however, why any quantum number 
should be equilibrated at all in a phase transition during dynamic non-equilibrium. 

It turns out that if one allows light quark phase space to be overpopulated, 
i.e. grow above equilibrium (7^ > 1), the entropy of the hadron gas starts growing 
rapidly with 7^ (Letessier et al. 2000) (Fig. 2.4). This region is phase space is 
kinetically not reachable through a hadron gas evolution such as that described in 
Eq. (2.2). However, if the particles are emitted from a rapidly non-perturbatively 
hadronizing qq rich QGP, over-saturation might occur. This is the only way found 
as yet to contain the entropy created during the QGP stage without reheating or 
expansion (which requires slow hadronization) The value of 7^ is bounded above by 
the Bose- Einstein condensation condition for pions 

^-.J,,T _ 1 < °°. 7,<e"-/- (2.83) 

above which a Bose-Einstein Condensate (BEG) of pions forms, which consumes 
energy while carrying no entropy. The BEG constraint works well with the super- 
cooling hypothesis, since 7^ has more freedom to rise to the point at which Sqgp < 
Shg only as the critical temperature becomes fairly small ~ 140MeV. 

This picture is phenomenologically in accord with the Gribov hadronization 
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Figure 2.4: Entropy, Energy and number density as 7^ > 1 increases . 

scenario described in section 1.3.1. As energy density reaches a critical value, the 
color potential between quarks and gluons grows up to the point where it is ener- 
getically convenient to fragment gluons into pairs to maintain large-scale color 
neutrality. The qq phase space is then overpopulated as the gluons are depleted by 
g ^ qq reactions. Tables 2.1 and 2.2 (Rafelski and Letessier 2003a) and (Letessier 
and Rafelski 2000) make it clear that fits where 7^ is a variable parameter consis- 
tently prefer values of 7^ just below the BEC limit. In this minimum, /x^ = Tln(A5) 
goes to a value very close to 0, in accordance to a scenario where strangeness is 
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Table 2.1: The chemical freeze-out statistical parameters found for nonequilibrium 
(left) and semi equilibrium (right) fits to SPS results. We show ^/snn, the temper- 
ature T, light quark fugacity A^, strange quark fugacity A^, the quark occupancy 
parameters 7^ and '^s/lq- Bottom line presents the statistical significance of the fit. 
The star (*) indicates for A^ that it is a value resulting from strangeness conserva- 
tion constraint. For 7^ that there is an upper limit to which the value converged, 
7^ < e"*''/^ (on left), or that the value of 7, = 1 is set (on right). 





17.2 


12.3 


8.75 


17.2 


12.3 


8.75 


T [MeV] 


135 ±3 


135 ±3 


133 ±2 


157 ±4 


156 ±4 


154 ±3 




1.69(5) 


1.98(6) 


2.56(6) 


1.74(5) 


2.03(7) 


2.69(8) 


K 


1.23* 


1.27* 


1.31* 


1.20* 


1.24* 


1.24* 




1.68* 


1.68* 


1.69* 


1* 


1* 


1* 


Is/lq 


0.91(6) 


0.83(4) 


0.85(6) 


0.66(4) 


0.60(4) 


0.67(5) 


xVdof 


11.4/6 


4.3/2 


2.3/4 


23/7 


8.9/3 


4.0/5 



produced through annihilation and gluon fusion reactions. These fit's statistical 
significance is markedly better than fits based on the equilibrium hypothesis. 



2.6 How to falsify freeze-out models 



Sudden non-equilibrium freeze-out has therefore passed the elementary comparison 
with experimental data. It can fit observed hadron yields to an acceptable statistical 
significance. It also seems that super-cooling can explain the small volume and fast 
emission times inferred by HBT observations (Csernai et al. 2003), though this has 
yet to undergo a quantitative corroboration. 

However, the sudden freeze-out picture allows for much more stringent fal- 
sifiabihty tests to be made, since it demands that all soft variables be described by 
the same set of statistical mechanics parameters. The rest of this thesis Performs 
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Table 2.2: freeze-out statistical parameters found for non-equilibrium (left) and 
7g = 1 ,7s fitted (right) fits to RHIC results. We show ^/snn, the temperature 
T, baryochemical potential /x^, strange quark chemical potential /is, strangeness 
chemical potential ^s, the quark occupancy parameters 7^ and ^s/jq, and in the 
bottom line the statistical significance of the fit. The star (*) indicates that there 
is an upper limit on the value of 7^ < e^'^l'^ (on left), and/or that the value is set 
(on right). 



^Satat [GeVj 


200 


130 


200 


130 


T [MeV] 


143 ±7 


144 ±3 


160 ±8 


160 ±4 


lib [MeV] 


21.5 ±31 


29.2 ±4.5 


24.5 ± 3 


31.4 ±4.5 


lis [MeV] 


2.5 ±0.2 


3.1 ±0.2 


2.9 ±0.2 


3.6 ±0.2 


m [MeV] 


4.7 ±0.4 


6.6 ± 0.4 


5.3 ±0.4 


6.9 ±0.4 




1.6 ±0.3* 


1.6 ±0.2* 


1* 


1* 




1.2 ±0.15 


1.3 ±0.1 


1.0 ±0.1 


1.13 ±0.06 


xVdof 


2.9/6 


15.8/24 


4.5/7 


32.2/25 


-Ptrue 


90%+ 


95%+ 


65% 


15% 



some of these falsifiability tests, which can be broadly grouped in three categories. 



Pcirticle spectra In the sudden hadronization picture, particles form from a trans- 
versely expanding system and undergo little or no re-interaction after forma- 
tion. Therefore, the same temperature which describes hadron abundances 
should, with the addition of transverse flow, also be able to describe the shape 
of hadron spectra. Moreover, because of the absence of the hadron gas phase, 
the number of hadrons in each particle species is unchanged since formation, 
except for the decay of short-lived resonances. 

Therefore, spectra should be normalized through hadronic chemical potentials, 
and the resonance admixture into the spectra should be prominently present. 
The non-equilibrium parameters 7^ fitted in table 2.2 should also be necessary 
to describe particle spectra. These predictions are tested in detail within the 
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next two chapters. 

Short-lived resonances Short hved resonances (strong excited states) provide a 
direct test of post-hadronization interactions, since their final yield will be 
affected by them. The sudden freeze-out model demands, therefore, that their 
abundance be determined by the same temperature and chemical potentials 
that determine long-lived particles. Any deviation from this can be attributed 
to post-hadronization dynamics. This is explored quantitatively in chapter 5. 

Exotica (Pentaquarks etc.) Pentaquarks should be strongly enhanced in a non- 
equilibrium model: Their mass seems to be similar to that of other short-lived 
resonances (Diakonov et al. 1997), but, of course, they have a different number 
of quarks (5 in the pentaquark case). Hence, in non-equilibrium models they 
should be enhanced with respect to similar baryons by a factor of ~ 7^. This 
is a large difference for 7^ ~ 1.5. This is also explored in chapter 5. 

Freeze-out dynamics If hadronization happens through a vacuum instability at 

the QGP-HG contact surface, the freeze-out hypersurface d^H^ acquires a 
physical significance connected to the confining phase transition. Therefore 

• Experimental data should strongly prefer a particular choice of (fT,^ 

• The preferred o?^E'^ should be localized on a surface moving with time 
across the fireball. 

Chapter 4 will explore whether the d^T,'^ choice can be made through an 
analysis of particle spectra alone. Chapters 5 and 6 will provide ways to 
constrain d^H^ further. 
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CHAPTER 3 
Particle spectra at SPS: a spherical ansatz 

3.1 Sudden freeze-out spectra 

The WA97 and NA57 experiments have, for the first time, measured the normahzed 
spectra of a wide range of strange hadrons, ranging in mass from the Ks to the Q 
(Bruno et al. 2003). The relative yield of these particles has been one of the main 
pieces of evidence pointing to the formation of quark gluon plasma at SPS energies 
(Heinz and Jacob 2000). Correspondingly the spectra's variety can provide a very 
stringent constraint on the system's freeze-out conditions. 

In particular, when sudden QGP breakup occurs, the spectra of hadrons 
are not formed at a range of stages in fireball evolution, but arise rather suddenly. 
Most importantly, particles of very different properties are produced by the same 
mechanism and thus are expected to have similar m^-spectra as is indeed observed 
(Bruno et al. 2003). The reported symmetry of the strange baryon and antibaryon 
spectra is strongly suggesting that the same reaction mechanism produces A and 
A and H and H. This is a surprising, but rather clean experimental fact, which we 
will interpret quantitatively in this paper. 

When the momentum distributions of final state particles stop evolving 
during the fireball evolution, we speak of thermal freeze-out. Because a spectrum of 
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strange hadrons includes directly produced, and heavy resonance decay products, 
one can determine the freeze-out temperature and dynamical velocities of fireball 
evolution solely from the study of precisely known shape of the particle spectra. 
We demonstrate this in some detail in section 3.5. The physical mechanism is that 
the freeze-out temperature determines the relative contribution of each decaying 
resonance while the shape of each decay contribution differs from the thermal shape, 
see section 3.3. 

We note that we make in our analysis the tacit assumption that practically 
all decay products of resonances are thermally not re-equilibrated, which is equiva- 
lent to the assumption of sudden freeze out. This is consistent with our finding that 
the m_L strange baryon and antibaryon distributions of A, A, S, S froze out near to 
the condition at which the chemical particle yields were established. 

One of the key objectives of this work is to present a comparison between 
thermal and chemical freeze-out analysis results for temperature, (explosion) col- 
lective velocity and other chemical and dynamical parameters. It is important to 
realize that particle spectra and yields are sensitive to magnitude of collective mat- 
ter flow, in which produced particles arc born, for somewhat different reasons: 1) in 
thermal analysis the collective flow combines with thermal freeze-out temperature 
to fix the shape of each particle spectrum, and temperature is also controlling the 
relative yield of contributing resonances - see previous paragraph for a here rele- 
vant tacit assumption - thus both T and v are fixed by the shape of 'm± data; 2) 
in chemical analysis the particle yields required are obtained integrating spectral 
yields, with experimental acceptance in p±,y implemented. Since many particles 
have a too small particle momentum to be usually observed, the acceptance-cut 
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yields used in chemical analysis depend quite sensitively on parameters which de- 
form the soft part of the spectra without changing the number of produced particles, 
such as is the flow velocity. For this reason precise particle spectra and yields are al- 
lowing to draw conclusions about the proximity of thermal and chemical freeze-out 
conditions. 



3.2 Spherically symmetric freeze-out 



Collisions at the SPS are characterized by a high degree of stopping power: Most of 
the primordial particles lose their longitudinal momentum in the initial collision, and 
the mid-rapidity region has a high particle/antiparticle imbalance (Appelshauser 
et al. 1999) and hence baryochemical potential. Hence, it is reasonable to assume 
that the system will not have a longitudinal rapidity structure at freeze-out. One 
such emission surface which has the advantage of analjd;ical simplicity is spherically 
symmetric freeze-out. Such an ansatz is also favoured by hydrodynamic evolution, 
since the hydrostatic force per unit volume is minimized. 

We proceed to construct a spherically symmetric freeze-out hypersurface, 
in which the freeze-out time depends on the radius alone 

/ tf{r) \ 
r sin(^) cos(0) 
r sin(^) sin(0) 
y rcos(^) J 

Parameterizing this surface in terms of r,9,4> we get, according to Eq. (2.15), an 



(3.1) 
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emission element of this form: 



d^E" = rdrde 



( 1 \ 

^ sin(^) cos(0) 
^ sin(^) sin(0) 



(3.2) 



combining this emission shape with a spherically symmetric flow we get an explicit 
form from the Cooper-Frye formula (where Q is the configuration space emission 
coordinate and the particle momentum coordinate). 

= ^W^Xai^ r^dr J sm{9)d9 (^E-pT^cos{9 - 0)^ e-T(^'-''^'r™«('^-^)) 

(3.3) 

Using modified Bessel functions, the integral over 9 can be done analytically 



and the Normalization N includes the degeneracy and chemical potentials 

N^QiUXai (3.5) 

This approach can be generalized to B-E and F-D statistics, relevant for tt, as, for 
realistic chemical potentials, these distributions can be represented as a converging 
series of Boltzmann-like therms 

oo 

Ae-^ ^ j:(±l)S,,,,A»«e-C«)^ (3.6) 



n=0 



For a real hydrodynamical system the integration over r will not be trivial, since 
flow will be a non-trivial, equation of state dependent function of r. For our fit, we 
shall average the system using one flow. 

The consequences of such averaging will be explored in the next chapter in detail 
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3.3 Thermal freeze-out analysis of SPS data 

In recent months experiment WA97 determined the absolute normahzation of the 
pubhshed ■m± distribution (Bruno et al. 2003), and we took the opportunity to 
perform the spectral shape analysis and will compare our results to those obtained 
in chemical yield analysis (Letessier and Rafelski 2000) in order to check if the 
thermal and chemical freeze-out conditions are the same. Our analysis continues 
and this report gives its current status. 

We report here a simultaneous analysis of absolute yield and shape of WA97 
results of six m_L-spectra of A, A, S, S, O + Q, Kg = {K^ + K^)/2 in four centrality 
bins. If thermal and chemical freeze-outs are identical, our present results must be 
consistent with earlier chemical analysis of hadron yields. Since the experimental 
data we here study is dominated by the shape of m_L-spectra and not by relative 
particle yields, our analysis is de facto comparing thermal and chemical freeze-outs. 

We have found, as is generally believed and expected, that all hadron m±_- 
spectra are strongly influenced by resonance decays. Thus we apply standard pro- 
cedure to allow for this effect (SoUfrank et al. 1990) and described in the previous 
chapter. 

Since particle spectra we consider have a good relative normalization, only 
one parameter is required for each centrality in order to describe the absolute nor- 
malization of all six hadron spectra. This is for two reasons important: 

a) we can check if the volume from which strange hadrons are emitted grows with 
centrality of the collision as we expect; 

b) we can determine which region in m±_ produces the excess of Q noted in the 
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chemical fit (Letessier and Rafelski 2000) is coming from. 

However, since the normahzation Vqgp common for all particles at given 
centrality comprises additional experimental acceptance normalization, we have not 
determined the value of the fireball emission volume at each centrality. Hence we 
will be presenting the volume parameter as function of centrality in arbitrary units. 

The best thermal and chemical parameters are found by minimizing the 
total relative error ^ as defined in Appendix B, where the data points are given by 
the WA97 (Bruno et al. 2003) spectra for K°, A, A, S, S, 17 + n. We have checked 
the validity of the statistical analysis by the usual method, i.e. omission of some 
data in the fit. 

Only in case of Kaons we find any impact of such a procedure. Noting 
that the statistical error of kaon spectra is the smallest, we have estabhshed how 
a a systematic error which could be for Kaons greater than statistical error would 
influence our result. For this purpose we assign to K'^ experimental results in most 
of our analysis an 'error' which we arbitrarily have chosen to be 5 times greater 
than the statistical error. In this way the weight of the kaon spectra in the analysis 
is greatly reduced. In the first result figure below (Fig. 3.1) we present both results, 
standard K° error and enlarged error. We see that while in individual result some 
change can occur, overall the physical result of both analysis are consistent. Thus 
we can trust in the combined study of hyperon and kaon data. This conclusion 
is reaffirmed in section 3.5, where we will see that the minimization of involve 
more or less pronounced minima, depending on the error size of the kaon spectra, 
see Fig. 3.12. In most calculations we present in this paper we will be using, unless 
otherwise said, hyperon results combined with the Kaon data with 5 times enlarged 
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• • SXStatistical K° error 

• • Statistical K° error 



chemical fit 
T limits 



12 3 4 

Multiplicity bin 

Figure 3.1: (color on-line) Thermal freeze-out temperature T for different central- 
ity bins compared to chemical freeze-out analysis shown by horizontal solid lines. 
Original statistical error is used in the dotted results, five times statistical error for 
kaon data is used in solid vertical lines. 

statistical error. We believe that in this way we will err on the conservative side in 
our physical conclusions. 

3.4 Overview of the results 

We show here a slate of results obtained within the approach outlined above. First 
we address the parameters determining the shape of the distributions, that is 

As function of the centrality bin we show in Fig. 3.1 the freeze-out tem- 
perature T of the m_L spectra. The horizontal fines defineate range of result of the 
most recent chemical freeze-out analysis , see Ref. (Letessier and Rafelski 2000). 
It is reassuring that we find a result consistent with the purely chemical analysis 
of data that included non-strange hadrons (Letessier and Rafelski 2000). There is 
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no indication of a significant or systematic change of T with centrahty. This is 
consistent with the beheve that the formation of the new state of matter at CERN 
is occurring in all centrality bins explored by the experiment WA97. Only most pe- 
ripheral interactions produce a change in the pattern of strange hadron production 
(Ambrosini et al. 1999). The (unweighted) average of all results shown in Fig. 3.1 
produces a freeze-out temperature at the upper boundary of the the pure chemical 
freeze-out analysis result, T ~ 145 MeV. It should be noted that in chemical anal- 
ysis = 'V (Letessier and Rafelski 2000), which may be the cause of this slight 
difference between current analysis average and the earlier purely chemical analysis 
result. 

The magnitudes of the collective expansion velocity v and the break-up 
(hadronization) speed parameter ^ are presented in Fig. 3.2. For v (lower part of 
the figure) we again see consistency with earlier chemical freeze-out analysis results, 
and there is no confirmed systematic trend in the behavior of this parameter as 
function of centrality. 

Though within the experimental error, one could argue inspecting Fig. 3.2 
that there is systematic increase in transverse flow velocity v with centrality and 
thus size of the system. This is expected, since the more central events comprise 
greater volume of matter, which allows more time for development of the flow. 
Interestingly, it is in v and not T that we find the slight change of spectral slopes 
noted in the presentation of the experimental data (Bruno et al. 2003). 

The value of the break-up (hadronization) speed parameter ^[= 
l/{drf/dtf)] shown in the top portion of Fig. 3.2 is near to velocity of light which 
is highly consistent with the picture of a sudden breakup of the fireball. This 
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Figure 3.2: (color on-line) Thermal freeze-out flow velocity v (top) and break up 
(hadronization) velocity ^ for different centrality bins. Upper limit ^ = 1 
(dashed hne) and chemical freeze-out analysis limits for v (solid lines) are also 
shown. 

hadronization surface velocity ^ was in the earlier chemical fit set to be equal to 
f , as there was not enough sensitivity in purely chemical fit to determine the value 

Unlike the temperature and two velocities, the overall normalization of 
hadron yields, must be, and is strongly centrality dependent, as is seen in 
Fig. 3.3. This confirms in quantitative way the believe that the entire available 
fireball volume is available for hadron production. The strong increase in the vol- 
ume by factor 6 is qualitatively consistent with a geometric interpretation of the 
collision centrality effect. Not shown is the error propagating from the experimental 
data which is strongly correlated to the chemical parameters discussed next. This 
systematic uncertainty is another reason we do not attempt an absolute unit volume 
normalization. 
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The 4 chemical parameters A^, Ag, 7^, 75/7^ are shown in the following Fig- 
ures 3.4,3.5. These parameters determine along with the final particle yield. 
Since we have 5 parameters determining normalization of 6 strange hadron spectra, 
and as discussed we reduce the statistical weight of Kaons, there is obviously a lot 
of correlation between these 4 quantities, and thus the error bar which refiects this 
correlation, is significant. 

The chemical fugacities A^ and A^ shown in Fig. 3.4 do not exhibit a sys- 
tematic centrality dependence. This is consistent with the result we found for T 
in that the freeze-out properties of the fireball are seen to be for the temperature 
and chemical potential values independent of the size of the fireball. Comparing 
to the earher chemical freeze-out result in Fig. 3.4 one may argue that there is a 
systematic downward deviation in A^. However, this could be caused by the fact 
that the chemical freeze-out analysis allowed for isospin-asymmetric S~((iss) yield 
(Letessier and Rafelski 2000), while out present analysis is not yet distinguishing 
light quarks. 

The ratio 75/7^ shown in bottom portion of Fig. 3.5 is systematically smaller 
than unity, consistent with many years of prior analysis: when 7^ = 1 is tacitly 
chosen, this ratio is the value of 7^ in analysis of strange baryons. We have not 
imposed a constrain on the range of 7^ (top of Fig. 3.5) and thus values greater 
than the pion condensation point 7* = e^'^l'^'^ ~ 1.65 (thick line) can be expected, 
but in fact do not arise. 

It is important to explicitly check how well the particle m^-spectra are re- 
produced. We group all bins in one figure and show in Figs. 3.6, 3.7, 3.8, 3.9 in 
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Figure 3.3: Hadronization volume (arbitrary units) for different centrality bins. 

sequence A, A, S, S. It is important to note that there are some significant devia- 
tions which appear to be falling outside of the trend set by the other measurements, 
-this occurs for A as well but remains invisible in the figure due to the smallness of 
the experimental error bar. Overall, the description of the shape of the spectra is 
very satisfactory. 

We also describe the Kodata extremely well, especially in the m± range 
which is the same as that for hyperons considered earlier, as is seen in Fig. 3.10. 
We recall that these results were obtained reducing the statistical significance of 
Kaon data, and thus the conclusion is that hyperons predict both the abundance 
and shape of kaon spectra. Moreover, all the strange hadron spectra can be well 
described within the model we have adopted. 
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Figure 3.4: Thermal analysis chemical quark fugacity (top) and strange quark 
fugacity (bottom) for different centrality bins compared to the chemical freeze-out 
analysis results. 

3.5 Statistical significance of the results presented 

We have analysed the validity and consistency of our data analysis by exploring 
X^/DoF profiles. These are obtained by fixing the value of one of the parameters 
(we consider Tf, v, dtf/drf — 1/|^), and computing the related x^/^oF, the total 
error divided by degrees of freedom. These are the number of measurements minus 
number of parameters, DoF is typically 33 in this data analysis. All curves must 
have the same x^/DoF at the minimum and this minimum must point to the value 
of parameters we report. For the temperature T we produced two results shown 
in Fig. 3.12, in the bottom section for experimental (statistical) measurement 
error, and in the top part for the five times enlarged error. We recall that both 
results are presented in Fig. 3.1. We note that there is a pronounced Xt/^oF 
minimum shown on logarithmic scale) for all 8 results of which the average value 
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Figure 3.5: Thermal analysis chemical quark abundance parameter 7^ (top) and 
Is/lq (bottom) for different centrality bins compared to the chemical freeze-out 
analysis. Thick line: upper limit due to pion condensation. 

is at T = 145 MeV. We show in Fig. 3.13 the profile of Xt/DoF for the collective 
flow velocity v (top) and the freeze-out surface dtjjdrj = motion (bottom) 

being fixed. These minima can be shown on linear scale. We note a mild secondary 
minimum in the region v ~ 0.25-0.35. However, the minima we find ai v — 0.5-0.58 
are by far more significant, dtf/drf = is converging to a sharp minimum seen 

in bottom portion of Fig. 3.13 at at a value consistent with the sudden breakup 
scenario. It is necessary to include dtf/drf = along with v in the analysis to 

find this result, which was not always done in other studies of particle spectra. 



3.6 Omega spectra 

In Fig. 3.14 all four centrality bins for the sum -|- O are shown. We see that we 
systematically under predict the two lowest m_L data points. Some deviation at 
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Figure 3.6: (color online) Thermal analysis tut spectra: A . 

high may be attributable to acceptance uncertainties, also seen in the the S 
result presented earlier in Fig. 3.9. We recall that there is a disagreement with the 
Omega yields in the chemical analysis, which thus does not include in the analysis 
the production of Q. In the here presented analysis we see that this disagreement 
is arising at low momentum. 

The low-m_L anomaly also explains why the inverse m± slopes for Q, Q are 
smaller than the values seen in all other strange (anti)hyperons. One can presently 
only speculate about the processes which contribute to this anomaly. We note that 
the l-2s.d. deviations in the low m_L-bins of the + i7 spectrum translates into 
3s. d. deviations from the prediction of the chemical analysis. This is mainly a 
consequence of the fact that after summing over centrality and m±, the statistical 
error which dominates + spectra becomes relatively small, and as can be seen 
the low m_L excess practically doubles the Q yield. 
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Figure 3.7: (color online) Thermal analysis rriT spectra: A. 
3.7 Discussion of SPS spectra 



Our thermal freeze-out analysis confirms that CERN-SPS results decisively show in- 
teresting and new physics, and confirms the reaction picture of a suddenly hadroniz- 
ing QGP-fireball with both chemical and thermal freeze-out being the same. In our 
view driven by internal pressure, a quark-gluon fireball expands and ultimately a 
sudden breakup (hadronization) into final state particles occurs which reach detec- 
tors without much, if any, further rescattering. The required sudden fireball breakup 
arises as the fireball super-cools, and in this state encounters a strong mechanical 
instability (Letessier and Rafelski 2000). Note that deep super cooling requires a 
first order phase transition. 

The remarkable similarity of m± spectra reported by the WA97 experiment 
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Figure 3.8: (color online) Thermal analysis rriT spectra: S. 

is interpreted by a set of freeze out parameters, and we see that production mech- 
anism of A, A, and S, S is the same. This symmetry, including matter-antimatter 
production is an important cornerstone of the claim that the strange antibaryon 
data can only be interpreted in terms of direct particle emission from a deconfined 
phase. 

The reader must remember that in presence of conventional hadron collision 
based physics, the production mechanism of antibaryons is quite different from that 
of baryons and a similarity of the m± spectra is not expected. Moreover, even if 
QGP is formed, but a equal phase of confined particles is present, the annihilation of 
antibaryons in the baryon rich medium created at CERN-SPS energy would deplete 
more strongly antibaryon yields, in particular so at small particle momentum, with 
the more abundant baryons remaining less influenced. This effect is not observed 
(Bruno et al. 2003). 
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Figure 3.9: (color online) Thermal analysis ttit spectra: S. 

Similarity of m^-spectra does not at all imply in our argument a similarity 
of particle rapidity spectra. As hyperon are formed at the fireball breakup, any 
remaining longitudinal flow present among fireball constituents will be imposed 
on the product particle, thus A-spectra containing potentially two original valence 
quarks are stretched in y, which A-y-spectra are not, as they are made from newly 
formed particles. All told, one would expect that anti-hyperons can appear with 
a thermal rapidity distribution, but hyperons will not. But both have the same 
thermal-explosive collective flow controlled shape of m^-spectra. 

We have shown that thermal freeze-out condition for strange hadrons 
(K^, A, A, S, H) agrees within error with chemical freeze-out and we have confirmed 
the freeze-out temperature T ~ 145 MeV. These findings about the similarity of 
thermal and chemical freeze-out were controversial, when the experimental single 
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Figure 3.10: (color online) Thermal analysis ttit spectra: Kg. 

particle spectra were lacking precision, since pion spectra and two particle corre- 
lation analysis did not yield this result. However, this paper studies the precise 
hyperon and kaon spectra which reach to relatively low p± and compares with 
definitive chemical analysis of SPS data. The two particle correlation analysis in- 
volves pions, which unlike strange hadrons here considered, are potentially witnesses 
to other physics than the properties of dense and hot quark-gluon phase. 

We were able to determine the freeze-out surface l/§7^ = dtf/drj dynamics 

dt 

and have shown that the break-up velocity -^^ is nearly the velocity of light, as would 
be expected in a sudden breakup of a QGP fireball. A study with dtfjdrf has not 
been previously considered, and only collective fiow is included in the description 
of the particle source. In our analysis we find a slight increase of the transverse 
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Figure 3.11: (color online) x^/Dof per particle data point. As can be seen, Ks (error 
given by the experiment) gives the dominant discrepancy . 

expansion velocity with the size of the fireball volume, but consistently v < l/Vs. 

We have reproduced the strange particle spectra in all centrality bins. Our 
findings rely strongly on results obtained by WA97 at smallest accessible particle 
momentum, and this stresses the need to reach to smallest possible p± in order 
to be able to explore the physics of particle freeze-out from the deconfined region. 
Moreover, we demonstrated that the experimental production data of + has a 
noticeable systematic low p± enhancement anomaly present in all centrality bins. 
This result shows that it is not a different temperature of freeze-out of + f2 
that leads to more enhanced yield, but a soft momentum secondary source which 
contributes almost equal number of soft i7 + compared to the systematic yield 
predicted by the other strange hadrons. 
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Figure 3.12: The total error divided by degrees of freedom for different centrality 
bins, shown as function of (fixed) freeze-out temperature T, bottom for the exper- 
imental value of the (statistical) K'^ error, top for the 5 times enlarged kaon data 
statistical errors. 
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Figure 3.13: The total error divided by degrees of freedom for different centrality 
bins, shown as function of (fixed) flow velocity v on top and for (fixed) freeze-out 
surface dtf/drf — l/§^ dynamics on the bottom. 
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Figure 3.14: (color online) Thermal analysis rriT spectra: f2 + fi. 
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CHAPTER 4 
Particle spectra at RHIC: A boost-invariant ansatz 

4.1 Boost-invariant statistical hadronization 
4.1.1 General remarks 

The advent of RHIC data has coincided with the consensus that a model based 
on statistical hadronization, combined with transverse expansion, can explain both 
the abundancies and the transverse momentum distributions of hadrons produced 
in heavy ion collisions. The fitted parameters, and in particular the temperature, 
however, have varied considerably, ranging from as low as 110 MeV (Burward-Hoy 
2003; Afanasiev et al. 2003) to 140 MeV (Rafelski and Letessier 2003d) to as high 
as 160 and 170 MeV (Baran et al. 2004; Braun-Munzinger et al. 2001). 

Such discrepancies are not very surprising, since the models differ consider- 
ably. However, this means that before we can say that the freeze-out temperature 
has been determined, we must understand precisely the origins of these differences, 
and try to ascertain which model is "the best" , both in terms of physical reason- 
ableness and as an ansatz to fit the data. We shall proceed to give an overview of 
the ways in which these models differ. 

Firstly, spectra are normalized very differently in different models. Some 
recent work fits the particle slopes only (Burward-Hoy 2003; Hohne et al. 2003; 
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Sandor et al. 2004; Adams et al. 2003b), treating the normalization of each par- 
ticle as a free parameter. This approach can be argued for assuming a long-lived 
post-hadronization "interacting hadron gas phase" in which individual hadron abun- 
dances subject to inelastic interactions evolve away from chemical equilibrium. This 
particular reaction picture clashes with e.g. the fact that short-lived resonance ratios 
can be described within the statistical hadronization model using the chemical (sta- 
tistical hadronization) freeze-out temperature obtained in stable particle studies ( 
(Braun-Munzinger et al. 2001) and next chapter). This implies that in principle the 
relative normalization of the particle spectra should be derived from a hadroniza- 
tion scenario involving flavor chemical potentials. In fact a study of RHIC spectra 
finds that the normalization can be accounted for (Baran et al. 2004), and that 
the chemical equilibration temperature also describes particle spectra well. This is 
suggesting that any post-hadronization re-interaction phase is short and has minor 
influence on the particle yields. 

The problem is that the different ways to derive hadronization particle 
distributions have a profound effect on the resulting fitted temperature. Tempera- 
ture affects the absolute number of particles through several mechanisms and anti- 
correlates with the phase space occupancy parameters 7^, i = u, d, s (Letessier and 
Rafelski 2000). It has been found that the introduction of these parameters, mo- 
tivated by the need to conserve entropy at hadronization (Letessier and Rafelski 
2000) decrease the x^/ per degree of freedom considerably and lowers the freeze-out 
temperature by 30 MeV (Letessier and Rafelski 2000). Other workers assume the 
light flavors are in chemical equilibrium (Becattini et al. 2003; Braun-Munzinger 
et al. 1999; Braun-Munzinger et al. 2001; Baran et al. 2004). 
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Additionally, when fitting the particle spectra, the system's spatial shape 
and the way the freeze-out progresses in time have a considerable effect on the form 
of particle distributions, and hence on the fitted temperature and matter flow. The 
impact of freeze-out geometry and dynamics on particle spectra were examined well 
before RHIC data became available (SoUfrank et al. 1990) and it was realized that 
an understanding of freeze-out is essential for the statistical analysis of the fireball 
(Heinz and Kolb 2002). Even though this matter has been clearly recognized, a 
systematic analysis of how freeze-out geometry affects particle distributions is for 
the first time attempted here. In fact, each of the models used in the study of particle 
spectra (Hohne et al. 2003; Adams et al. 2003b; Burward-Hoy 2003; Baran et al. 
2004) employs a different choice of freeze-out geometry, based on different, often 
tacitly assumed, hadronization scenarios. Thus an understanding for the influence 
of hadronization mechanism is impossible to deduce from this diversity. 

4.1.2 Freeze-out geometry 

At RHIC collision energies the measured dN/drj (Olszewski et al. 2002; Staszel et al. 
2002) indicates that around mid-rapidity the system conditions can be approximated 
by the Bjorken picture (Bjorken 1983). 

To describe particle spectra measured around mid-rapidity, therefore, boost 
invariance becomes the dominant symmetry on which freeze-out geometry should 
be based. We shall construct a general hadronization scenario, combining the cut 
Cooper-Frye formula (Eq. (2.18)) with boost-invariance. 

A particle's energy at rest with respect to the collective flow is given by 
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combining the flow fleld (parametrized by the transverse and longitudinal rapidi- 
ties yL,T as well as the emission angle 0) with the observed particle's momentum 
(parametrized by the rapidity y, transverse momentum px and angle 0) 



^ cosh(i/i) cosh(i/T) ^ 
sinh(|/7^) cos(6') 
sinh(yr) sin(^) 
y sinh(yL) cosh(yT) ) 

The rest energy will then be 



^ rriT cosh(|/) ^ 
Pt cos(0) 
PTSin(0) 
\^ mT sinh(y) J 



(4.1) 



Pi^u'^ = rriT cosh(|/T) cosh(|/ - y^) - sinh(|/T) cos(6' - (j))pT. 



(4.2) 



The requirement for the Bjorken picture is that the emission volume element has 
the same dependence: 



p^d^E^ ~ A cosh(7/ - vl) + B. 



This constrains the freeze-out hypersurface to be of the form 



(4.3) 



= {tfCOsh{yL),x,y,tfSmh{yL)). 



(4.4) 



Here is a parameter invariant under boosts in the z direction, whose physical 
significance depends on the model considered. 

For central collisions, a further simplifying constraint is provided by the 
cylindrical symmetry, which forces tf, as well as yz and yr, to be independent of 
the angles 9 and (p. The freeze-out hypersurface can be parametrized, in this case, 
as 

W — (iy(r) cosh(7/i),rsin(^),rcos(^),i/(r) sinh(yi)) , (4.5) 

^ (4.6) 



dt dt 

d^'Z^' = tjrdrdOdyL ( cosh(|/i)-^ cos(6'), sin(6'), sinh(|/£,) 



dr 
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And the emission element takes the form 



rriT cosh(y - y^) - cos{9 - 0) 



tfrdrdOdyL, 



(4.7) 



with the same dependence on the angle as Eq. (4.2). 



Putting Eq. (4.2) and Eq. (4.7) together, the particle's momentum distri- 
bution function in the Boltzmann approximation is given by 



E— 
d^p 



dtf 

rriT cosh(y - i/l) - Pt^ cos(6' - 0) 



(4.8) 



exp 



rriT cosh(|/T) cosh(|/ — yi) — sinh(|/T) cos(6' — 0)pTn 



T 



t frdrdOdyi 



(the quantum statistics case, relevant for tt, can be obtained through the prescription 
in Eq. (3.6)). For an approximately boost-invariant cylindrical system we can use 
Bessel functions to perform the integrals over and 9 analytically 



Kn{z] 



oo 

cosh(ny)e-^°^^(2')dy 



I^(z)^— [ cos(n^)e-^°^(^)d^ 
27r Jo 



and obtain 



dN 



dt, 



^'¥p ^ '^Th{oiPT)Ki{(3mT) - pT-Q^Ii(apT)Ko((3mT) 



(4.9) 
(4.10) 

(4.11) 



What distinguishes the models currently considered is the time component 
of the freeze-out surface. The most general freeze-out hypersurface compatible with 
cylindrical symmetry is provided by Eq. (4.5). Generally, (a generic function of 
r) represents the time, in a frame co-moving with the longitudinal flow, at which 
the surface at distance r freezes out. 

The fits in Refs. (Burward-Hoy 2003; ?; Adams et al. 2003b) are based on 
a particular case of such a freeze-out surface, in which tf is completely independent 
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of r (dtf/dr — O). Such a picture's physical reasonableness can be questioned, e.g., 
why should spatially distant volume elements, presumably with different densities 
and moving at different transverse velocities, all freeze out simultaneously in a 
longitudinally co-moving frame?. However, such a simple model can perhaps serve 
as an approximation. 

More generally, the "burning log" model, discussed in the last chapter for 
the spherically symmetric case, assumes that the emission occurs through a three- 
dimensional hadronization surface which is moving at a constant "velocity" {dtf/dr 
throughout the fireball. Both boost-invariant and spherically symmetric versions 
of burning log model were considered. Even if the hadronization velocity encom- 
passes an extra parameter, it is worth considering since it is based on a physically 
motivated hadronization picture. Moreover, the burning log picture is a suitable 
framework in the study of sudden hadronization. Sudden hadronization occurs when 
the fireball encounters a mechanical instability, which combined with the fireball's 
high transverse fiow ensures that the emission surface spreads to the interior of the 
fireball with dtf/dr ~ c. All of the indications suggested for such a picture seem 
to be borne out by both SPS and RHIC data (Letessier and Rafelski 2000; Rafelski 
and Letessier 2003d). 

An approach based on the hypothesis of initial state "synchronization" by 
the primary instant of collision and the following independent but equivalent evo- 
lution of all volume elements assumes that each element of the system undergoes 
freeze-out at the same proper time r. In this framework each fireball element ex- 
pands and cools down independently, hadronizing when its temperature and density 
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Figure 4.1: (Color on-line) While boost-invariance fixes the longitudinal freeze-out 
structure (left), several scenarios exist for the transverse dependence of freeze-out 
(right). For spherical freeze-out, only plot on the right applies . 



Table 4.1: (Color on-line) Freeze-out hypersurfaces at contours of constant radii. 



Surface Ej^ i 


Constant tr ( tf \ i n \ t i \ 
Q^^jQ^^i^ \r) mTKi{(3mT)Io{apT) 


Hubble 
(constant r/) ''"■^ 


f cosh(?/i) cosh(yT) \ 
sinh(|/T) cos(6') 
smh^yx) sin(6') 

^ sinh(?/i) cosh(yT) J 


rriT cosh.{yT)Io{apT)Ki{PmT) — 
Pt smh{yT)Ii{apT)Ko{pmT) 


burning log 
(boost invariant) 


/ tf{r) coshlyi) \ 

r cos{9) 

r sin(^) 
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PT^h («Pt) ^0 iPrriT) 
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reach the critical value. This model was successfully used to describe RHIC rriT- 
spectra (Baran et al. 2004). In this approach t/ in Eq. (4.5) is equal to rcoshd/y) 
and the hadronization hypersurface in Eq. (4.6) becomes proportional to the flow 
vector: 

E'* = Tu^ (4.12) 
d^Ef" = TrdrdedyLuf" = dVu^ (4.13) 
r = Tsinh(yT)- (4-14) 

In this hadronization model the heavy ion fireball behaves similarly to the expand- 
ing Hubble universe. In the 'Hubble' scenario, the Cooper-Prye formula reduces 
to the Touscheck Covariant Boltzmann distribution (Letessier and Rafelski 2002; 
Touschek 1968; Hagedorn et al. ). 

||fe-^/^ - ^,d'p2S,{p'-m')e-^''-^/^ (4.15) 
V^" = Vou^ (4.16) 

(Where V is the co moving fireball's volume element in the local rest frame.) 

To summarize and illustrate the diversity of distinct hadronization geome- 
tries we present in Table 4.1 and Fig. 4.1 the freeze-out scenarios examined here. As 
we shall see the choice of freeze-out geometry produces in a fit of experimental data 
a non trivial effect capable of altering significantly the understanding of statistical 
hadronization parameters. 

4.1.3 Flow profile 



Hydrodynamical expansion of the fireball implies in general that each volume ele- 
ment will have a different density and transverse expansion rate. For this reason, 
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Figure 4.2: (Color on-line) 7r,K,A and S m^- distributions obtained with different 
freeze-out models and flow profiles. For this and subsequent figures, a uniform 
density profile was assumed . 

the integral over d^Tj can span a range of flows, weighted by density. In first ap- 
proximation one can fit data using just an "average" flow velocity throughout the 
entire fireball (Sandor et al. 2004; Torrieri and Rafelski 2001a): 

= l^dr{E-PT^)fiT,yAr),X) 

oc (i?-p^Jl)/(T,(yT),A). (4.17) 
dr 

However, if one wants to properly identify the flow profile should be 
taken into account. Hydrodynamic simulations (Teaney et al. 2001) accompanied by 
assumption that freeze-out happens when a volume element reaches a critical energy 
density indicate that the transverse rapidity will depend linearly with the radius i.e. 
vt ~ tanh(r). This condition, however, is appropriate for a static freeze-out and will 
not in general hold if the freeze-out is sudden. Other flow profiles have been tried 
in the literature, arising from dynamical hypothesis. For example, the assumption 
that the freeze-out occurs at the same time t/ results in a quadratic {v oc r^) flow 
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profile (Schnedermann and Heinz 1993), which has also been used recently in fits to 
data (Adams et al. 2003b). In the Hubble fireball (Baran et al. 2004) the freeze-out 
conditions will also result in a distinctive flow profile. Specifically with T,^ oc u^, 
we have ^yv oar. 

Density profiles also depend on the assumed initial condition and the equa- 
tion of state of the expanding QGP. It has been shown (Aggarwal et al. 1999) that 
different density choices have a considerable effect on both the temperature and 
flow fits at SPS energies. 

Fig. 4.2 shows how the choice in hadronization dynamics and flow profiles 
at same given freeze-out temperature and transverse flow can result in a range of 
inverse spectral slopes. Here the density profiles were assumed to be uniform. It 
is clear that the same freeze-out parameters give rise to a variety of substantially 
different particle spectra. Conversely, fits to experimental data will only produce 
reliable information on the freeze-out conditions if and when we have a prior knowl- 
edge of the hadronization geometry and dynamics. Therefore, conclusions about 
statistical model fits, as well as arguments whether freeze-out occurs simultane- 
ously for different particles or not, cannot be answered while the models used to fit 
the data are plagued by such uncertainties. 

We now proceed to use a Monte-Carlo method to calculate the effect these 
uncertainties have on fitted data. Afterwards, we will apply a range of fits to RHIC 
models to see if the "right" hadronization model can be isolated from particle spectra 
alone. 
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Figure 4.3: (color on-line) Results of fits to Monte-Carlo generated data samples I 
and II = 0): profiles for temperature (left) and fitted dtf/dr (right). The 
models used in the fits are described in sections 1.1 and 1.2 Full fit includes fitted 
chemical potentials and 7g,s, resonances and fitted dtf/dr. Last profile shows effect 
of fitting sample II using Eq. (4.17) . 

4.2 Sensitivity to model choice 

The ambiguities presented in the previous section mean that it is important to study 
their effect on the statistical model's fitted parameters. One way to do this is to use 
a Monte-Carlo to generate data according to a particular freeze-out model, and to 
see what happens if the "wrong" model is used to perform the fit. We have written a 
Monte-Carlo program which can be used for this purpose. An acceptance /reject ion 
algorithm is used to generate particles in a statistical distribution in the volume 
element's rest frame. The accepted particles are then Lorentz transformed to the 
lab frame. (Any flow and density profile, as well as any freeze-out surface can 
be accommodated). Resonance decays are handled through Eq. (2.34), using the 
MAMBO algorithm (Kleiss and Stirling 1992) to generate points in phase space. 
Output can be used to generate spectra or fed into a microscopic model such as 
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Figure 4.4: (color on-line) dNtot/d{y,v) arising from freeze-out of a Bjorken fluid, 
before (open circles) and after (solid circles) resonance decays. The normalization 
is arbitrary, chosen to coincide with PHOBOS data (triangles) . 



uRQMD (Bass et al. 1998). 

The Monte-Carlo output was used to produce the data points in Fig. 4.4. It 
can be seen that a boost-invariant statistical hadronization can explain the global 
properties of the system such as dNtot/drj. It can also be verified that the role 
of resonances is absolutely crucial. We proceeded to generate three datasets of 
particles. Each data set had a temperature of 140 MeV, a maximum transverse 
flow of 0.55, and out of equilibrium chemistry (7^ = 1.4, 7^/75 = 0.8). Generated 
particles include vr, i^'",/),/). A, S, f2 and their resonances. The three samples 
differ in their choice of freeze-out geometry (specifically dtf/dr) and flow profile: 
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Sample I dtf/dr = and no flow profile, as fitted in (Sandor et al. 2004) 

Sample II dtf/dr = and a quadratic flow proflle, as fltted in (Burward-Hoy 
2003; Hohne et al. 2003; Adams et al. 2003b) 

Sample III dtf/dr = 1, the boost-invariant analogue of (Torrieri and Rafelski 



Figure 4.5: (color on-line) Results of fits to Monte-Carlo generated data samples III 
(-^ = 1): profiles for temperature (left) and fitted dtf/dr (right). The models 
used in the fits are described in sections 1.1 and 1.2 Full fit includes fitted chemical 
potentials and •jg^s, resonances and fitted dtf/dr. 

We have fitted the three samples to a variety of models, producing x^/DoF profiles 
for freeze-out temperature and the dtf/dr parameter. Fig. 4.5 shows the profiles 
resulting in the fit to sample III while samples I and II are shown in Fig. 4.3. A 
full chemistry model with resonances (solid black line) seems to be equivalent (as 
far as the position of the minimum and the value of x^/DoF) to a fit in which 
normalization is particle-specific (dashed red line). However, the chemical potential 
fit has greater statistical significance since considerably more degrees of freedom are 
required for arbitrary normalization. 
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If chemical potentials are included resonances become essential since a fit 
with chemical potentials but no resonances (dotted blue line) loses all statistical 
significance. Similarly, the physical presence of non-equilibrium (7^ ,5 ^ 1) means 
chemical potentials have to include the non-equilibrium parameter for the fit to 
be meaningful (black dot-dashed line). The freeze-out geometry does not seem 
to impact the temperature minimum that much. However, the correct freeze-out 
geometry can be picked out by a comparison of fits to different models by choosing 
the model with the lower x^/DoF. Moreover, data sensitivity to temperature is 
strongly affected by freeze-out geometry: Comparing the temperature profiles for 
different choices of dtf/dr (Figs. 4.3 and 4.5) it is apparent that the temperature 
minimum is more definite in the dtf/dr = case. In the case of explosive freeze-out, 
the correlation between temperature and other parameters in the fit (notably flow) 
increases, resulting in a shallow increase at larger than minimum temperatures. 

Finally it should be noted that flow proflle, freeze-out geometry and tem- 
perature appear to be strongly correlated. If data sample II is fltted with a distri- 
bution with no flow proflle (such as Sample I) there is a non-negligible shift in both 
the fltted temperature and dtf/dr, and a small rise in x^/DoF (Fig. 4.3, brown 
dot-dashed fine). 

As flts including both flow proflle and resonances are computationally very 
intensive, we shall have to limit ourselves to the observation that the uncertainty 
in flow proflle is a source of systematic error here. In the following chapter, a probe 
capable of disentangling flow proflle from freeze-out will be developed. 
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Figure 4.6: (color on-line)Fit to RHIC data profiles for temperature (left) and fitted 
dtf/dr (right). The models used in the fits are described in sections 1.1 and 1.2 
Full fit includes fitted chemical potentials and 'jq^s, resonances and fitted dtf/dr. . 

4.3 Fit to RHIC data = 130GeV) 

Finally, we have performed a fit to the available RHIC data. The data sample we 
used is the same as the one used for the Monte-Carlo, but, since the STAR and 
PHENIX spectra had different trigger requirements (notably centrality) and accep- 
tance regions, we have used two different system volumes, one for STAR particles 
and another one for PHENIX. 

The results are similar to the Monte Carlo data in many ways. The x^/DoF 
was only slightly larger. 

A fit with particle-specific normalization gives a very similar ^ind fitted tem- 
perature to a fit including chemistry and resonances. If chemistry is included in 
particle spectra analysis, resonances and non-equilibrium are essential. The fit to 
freeze-out geometry weakly points to dtj /dr = 1, a. picture that is supported by the 
temperature profile, virtually identical to data sample HI. However, we can not 
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claim our study to be complete in this respect, since we have not yet investigated 
the effect of including flow profile in the models. As Monte-Carlo simulations have 
shown, the conclusion can differ once these are taken into account. 

Fig. 4.7 shows the hyperon and pion spectra from the global fit of Fig. 4.6. 
A comparison of the fits on the left panel confirms that a model with no chemistry 
is about as good at fitting particle spectra as a model with resonances and chemical 
potentials. However, the second fit also has predictive power: We have used the 
fitted parameters to predict the ttit spectrum for the E*. Unsurprisingly, we found 
that the E* should have roughly the same slope as the Hs, but its total multiplicity 
should be about three times as big. We therefore suggest that a greater sample 
of spectra, in particular more spectra of heavy resonances taken within the same 
centrality bin as light particles (to make sure both flow and emission volume are 
the same for each particle) would help in establishing whether chemical potentials 
are a good way to normalize hadron spectra or not. 

The only spectrum which presents a significant systematic deviation for 
most models is the n". As Fig. 4.7 (right panel) shows, most models fail to catch 
the upward dip of the low momentum pions, and indeed simple transverse expansion 
predicts the reverse trend (Burward-Hoy 2003; Hohne et al. 2003; Sandor et al. 
2004; Adams et al. 2003b). Including resonances, and allowing for 7^ > 1 helps 
(the latter is equivalent to postulating a pion "chemical potential" (Tomasik et al. 
2003)) However, to fully account for the lowest momentum pions, even addition of 
resonances and 7g > 1 are not enough. One has to add a power-law component to 
the pion spectrum 




(4.18) 
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This contribution (roughly 6% of the total pion yield in the best fit) also accounts 
for the highest px pions. Such a parametrization has been justified (Hagedorn 1984) 
and successfully used (Peitzmann 2003) before. 
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Figure 4.7: (color online)Left: S and Q uit distributions, within a global fit includ- 
ing resonances and chemical potentials (Fig. 4.6, solid black) as well as no resonances 
or chemistry ( Fig. 4. 6, dashed red). The model with chemical potentials was used 
to predict the S* itit distribution using the Monte-Carlo. Right: PHENIX vr" px 
distribution, within the global fit from Fig. 4.6. As can be seen, both resonances 
and 7g help, but are not sufficient to explain the pion distribution fully. 



Finally, Fig. 4.8 (left) shows the Blast wave fit to the K* momentum distri- 
bution. The fact K* can be fitted to the same temperature and flow as the stable 
hadrons is, in itself, remarkable, since it is a short lived resonance detected by in- 
variant mass reconstruction, which potentially decays before freeze-out and whose 
decay products can be rescattered (and made undetectable) by other particles in the 
medium. We shall not dwell on this for too long, since the next chapter is mainly 
devoted to a discussion on directly detectable short-lived resonances. However, it 
should be mentioned that, as Fig. 4.8 (right) shows, microscopic models predict the 
rescattering of decay products to be non-trivial and strongly momentum dependant 
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Figure 4.8: (color online) (left) K* ^ K~ and (f) normalized spectra within the global 
fit using non-equilibrium sudden freeze-out ( Fig. 4.6, black), (right) momentum 
dependence of K* suppression through rescattering of decay products, calculated 
in the microscopic hadronic model uRQMD . 

(Bleicher and Stocker 2004). Given this, the relative goodness of the fit to K* data 
is a strong boost for the sudden freeze-out picture. 

The acceptable K* description within a fit which includes and K~ 
is encouraging for another reason: Resonances and absolute normalization in fits 
with many particles remove the correlation between temperature and fiow: It is 
not possible to adjust the slopes by shifting the temperature and transverse fiow 
along the same contour, since in this case normalization is also affected. However, 
temperature is still strongly correlated with chemical potentials (expecially 7^ and 
7g). However, if the fit includes K*^ and Ks, this correlation is removed, since 
both spectra share a common chemical potential. Shifting the temperature from 
the minimum will result in a shift of one of these spectra away from the data, no 
matter what chemical potential is used. For this reason, a common tempeature 



131 

minimum in a fit which includes K* and Ks is an important test the sudden freeze- 
out model successfully passed. We hope further tests, in the form of more resonance 
spectra, are forthcoming. 

In conclusion, we have given an overview of different hadronization models, 
and described how they arise from different freeze-out scenarios of a system forming 
from a thermalized quark gluon plasma. We have also shown, in the boost invariant 
limit, how the various hadronization ansatzes give rise to quantitative differences in 
observed particle spectra. We have then used MonteCarlo simulated data to study 
the sensitivity to model choice of presently available experimental data, and have 
evaluated different models ability to fit presently existing RHIC data. While data 
slightly favors a chemical non-equilibrium explosive freeze-out, there is not enough 
evidence to make a definitive conclusion about this issue. We hope the forthcoming 
200 GeV results will clarify this further. 
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CHAPTER 5 
Direct detection of short-lived resonances 

5.0.1 Introduction and motivation 

The experimental measurement of short-hved hadron resonances can potentiaUy be 
very useful in clarifying some of the least understood aspects of heavy ion collisions. 
In general, evolution of a hot hadronic system proceeds according to Fig. 5.1. When 
mesons and baryons emerge from a pre-hadronic state, presumably quark gluon 
plasma, their abundances are expected to be fixed by hadronization temperature and 
chemical fugacities. This stage of fireball evolution is commonly known as chemical 
freeze-out. After initial hadronization, the system may evolve as an interacting 
hadron gas. At a certain point (which can vary according to particle species), 
thermal freeze-out, where hadrons stop interacting, is reached. 

A quantitative understanding of the above picture is crucial for any mean- 
ingful analysis of the final state particles. Many probes of deconfinemcnt arc most 
sensitive when the dense hadron matter fireball breakup is sudden and re- interaction 
time short or non-existent. Final state particles could, however also emerge re- 
membering relatively little about their primordial source, having been subject to 
re-scattering in purely hadronic gas phase. In fact, theoretical arguments have been 
advanced in support of both a sudden reaction picture (Csorgo and Csernai 1994; 
Rafclski and Letessier 2000) and a long re-interaction timescale (Braun-Munzinger 
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Position 



Figure 5.1: Stages of the space-time evolution of a heavy ion coUision. At a certain 
moment in proper time (known as chemical freeze-out), hadrons emerge. The system 
then evolves as an interacting hadron gas, until thermal freeze-out, the point at 
which all elastic interactions cease. 

et al. 1999; Heinz and Jacob 2000). Both pictures have been applied to phe- 
nomenological fits of hyperon abundances and distributions (Akkelin et al. 2002; 
Arbex et al. 2002; Torrieri and Rafelski 2001a). 

It is apparent that hyperon resonances can be crucial in resolving this am- 
biguity: their initial abundance, compared to stable particles with the same quark 
composition, depends primarily to the temperature at hadron formation. However, 
the observed abundance will be potentially quite different, and will strongly depend 
on the lifetime of the interacting hadron gas phase. Resonances can only be ob- 
served via invariant mass reconstruction, and their short lifetime means that they 
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can decay within the interacting hadron gas (Fig. 5.2). In this case, the decay 
products can undergo re-scattering, and emerge from the fireball with no memory 
about the parent resonance. Thus, the observable resonance abundance is sensitive 
to precisely those parameters needed to distinguish between the sudden and gradual 
freeze-out models. 

In this chapter, we start with a review of presently available experimental 
data, and the open questions which arise. We then proceed to describe how to cal- 
culate the initial resonance abundance and the effect of re-scattering in terms of the 
hadronization temperature and the lifetime of the interacting hadron gas. We show 
how these two parameters can be extracted from the experimental observations. 
Finally, we discuss possible answers to experimental challenges raised in the first 
section, and suggest ways by which these questions could be resolved by further 
measurements. Most of this write-up is based on recently published experimental 
measurements (Friese 2002; Fachini 2002; Adler et al. 2002b; Markert et al. 2002a; 
Afanasev et al. 2002; Zhang 2004; Salur 2004) and theoretical papers (Torrieri and 
Rafelski 2002; Rafelski et al. 2001; Torrieri and Rafelski 2001b) 
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Figure 5.2: How re-scattering can inhibit resonance reconstruction. If the resonance 
decay products reach the detector without further interactions, their invariant mass 
distribution should yield a clear peak at the resonance mass. However, if these 
decay products undergo re-scattering before reaching the detector, the signal may 
be indistinguishable from the background caused by unrelated particle pairs. 



5.1 Can short-lived resonances be explained by a statistical model? 



Since its beginning, the RHIC heavy ion program, and in particular the STAR 
collaboration, has managed to measure an impressive amount of resonance data. 
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Figure 5.3: (Color on-line) Fits to particle ratios containing several short-lived 
resonances, using an equilibrium (left) and a non-equilibrium (right) formalism. 
The fit in (right), was obtained using the SHARE program (See Appendix D). 

In addition to the K* and the A(1520), the p, /o(980). A, E* and S* have been 
measured in p-p, and sometimes Au-Au (Fachini 2004; Zhang 2004). Further results 
in Au-Au collisions are expected soon. Including these particles in a global fit, of 
the type described in section 2.2.1 places a tight constraint on the temperature, 
since many of the resonances have the same chemical composition but a different 
mass from observed "light" particles. Fig. 5.3 shows what happens when such a fit is 
done. A chemical equilibrium model (Braun-Munzinger et al. 2001) fails to describe 
short-lived resonances. Under chemical non-equilibrium, however (Fig. 5.3 right), 
both short-lived and long-lived resonances are described with equal acceptability. 
In a sense, the proponents of both a long freeze-out stage (which would pull short- 
lived resonances out of equilibrium ) and sudden freeze-out can be satisfied by this 
result: Proponents of hadronization at T=170 MeV can rely on the hypothesis 
that in-medium interactions will make short-lived resonances unobservable, while 
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the short-lived data points are well described by the non-equilibrium model, at the 
same temperature and jg^s as found in earlier fits (Rafelski and Letessier 2003b). We 
can try to see which of the pictures is right by concentrating on a few particularly 
simple resonances, and attempting to quantify the effect of a post-hadronization 
hadron gas phase. This topic will be addressed in the next section. 

5.2 Modeling resonances in heavy ion collisions 
5.2.1 Direct production at hadronization 

We assume that at hadronization, a volume element will be at thermal and chemical 
equilibrium in its local rest frame. As shown in section 2.2.1, this means that details 
of flow and freeze-out hypersurface cancel out when the global (4 tt) ratio, or a mid- 
rapidity ratio in a Bjorken-type fireball, are considered. 

In the Boltzmann approximation (applicable here since XparUcie ~ 0(1) and 
m >> T), the chemical potential will then cancel out when two particles of the same 
quark composition are compared. The ratio of the heavier resonance N* — > Nn to 
the lighter resonance will then be given by Eq. (2.24)) 

A^ g*m*^K2 (f^) + gm^K2 (f ) ^ ' ' 

(Here g is the statistical degeneracy). Ratios where the Boltzmann approximation 
is applicable and feed-down from particles of different quark constituents is small 
include E7A and A(1520)/A. K*°/K+ and 'K*^/K- has a larger correction due to 

KK^, while in other, recently measured ratios (A/p, p/tt (Zhang 2004)) the 

^Of course, the K*'^/K~ ratio can be corrected for (p K^K~ feed-down, since the NA49 
and STAR experiments have published yields for both. We recommend that such a correction 
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feed-down from particles with a different chemical potential is simply too large to 
be neglected. 

If less than 4 tt acceptance is observed, we can obtain the equivalent of 
Eq. (5.1) by integrating Eq. (3.4) and Eq. (4.11) across a finite y—pr range. As will 
be shown, the finite acceptance effects will cancel out to a very good approximation 
if particles of a comparable mass are considered. Hence, Eq. (5.1) can still be used 
for ratios with a finite acceptance detector. 

Table 5.1 summarizes the decay processes considered in our analysis and 
their parameters (Clebsh- Gordon coefficients have been used to estimate decays 
such as (N*^ N+7r-)/(N*^ A^°7r°)). 

Table 5.1: Resonances contributing to A and K production, with their degeneracies, 
rest-frame momentum (p*) and possibility for experimental reconstruction. 
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Reaction 


P* 


branching 


visible? 


^ 4 


E*°(1385) ^ E°7r° 


127 


^4% 


No 


8 


E*±(1385) ^ Ati^ 


208 


88% 


Yes 


4 


E*0(1385) ^ AttO 


208 


88% 


No 


2 


EO ^ A7 


74 


100% 


No 


4 


A(1520) ^ NK 


244 


45% 


Yes 


3 


0(892) ^ K+n- 


291 


67% 


Yes 



In Fig. 5.4 we show the relative thermal production ratios at chemical freeze- 
out over the entire spectrum of rapidity and ttit (solid lines) and central rapidity 
range (dashed lines). The sensitivity of resonance yields to hadronization tem- 
perature is apparent for all resonances under consideration. In particular, the E* 

emerges as a very promising candidate for further study. For example, at the lowest 

be done when presenting future K* results, as their interpretation becomes considerably more 
model-independent . 
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current estimates (T ~ 100 MeV) of the final break up temperature in 158^4 GeV 
SPS collisions 33% of A s are actually primary E* s, the percentage rises to slightly 
more than 50% if chemical freeze-out occurs at T = 190 MeV. For all of these cases, 
the E* primary yield should be considerably bigger than the S (if the chemical po- 
tentials in (Rafelski and Letessier 2003b) are used for the prediction), and the angle 
between the decay products should be about the same as in the Q — > AK decay. 
The NA49 detector in the SPS and the STAR detector at RHIC are both capable 
to measure hyperon and resonance yields well within the precision required to make 
a good measurement of the hadronization temperature. 

We also note that while the experimentally observed K* yield is compati- 
ble with the produced ratio in the case of non-equilibrium hadronization (Rafelski 
and Letessier 2003b), the A(1520) seems very suppressed. To discuss this further, 
however, an estimation of the effect of re-scattering on resonance abundances is 
necessary. 

5.3 Rescattering 

As explained in the introduction, direct observation of resonances relies on invariant 
mass reconstruction. Therefore, to calculate the observed resonance abundances, 
rescattering after hadronization will have to be taken into account. This can be 
looked at within a microscopic model of hadronic matter such as uRQMD. Here, 
we present such a study using a "back of the envelope" model which nevertheless 
seems to provide an acceptable quantitative description of how the propagation 
of resonances and their decay products through opaque matter deviates from an 
equilibrium prediction. 
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Figure 5.4: Temperature dependence of ratios of E*, K*^ and A(1520) to the total 
number of observed K A s and S s. Branching ratios are included. Dashed lines 
show the result for a measurement at central rapidity Ay = ±0.5. The experimental 
measurements included in the diagram are discussed in the first part of this paper. It 
should be noted that the K* /K ratio is actually K*/ K~ for the NA49 measurement, 
and an average {K* + 'K*)/{K+ + R-) for STAR. 

We first note in Fig. 5.4, that the relative signal is remarkably inde- 

pendent (within 5 %) of Temperature. This is because the S*(1530) contribution 
cancels nearly exactly the thermal suppression of the S originating in the S — S* 
mass difference. This effect could be used for a direct estimate of the S* lost 
through rescattering, even without knowledge of the freeze-out temperature, should 
the chemical parameters (Aj and 7i) be independently known (They can be obtained 
from a fit of the stable particle ratios). A simple test of sudden hadronization model 
consists in measurement of the ratio If it is significantly smaller than unity, 

we should expect a re-equilibration mechanism to be present. Otherwise sudden 
hadronization probably applies, since S* emerge from chemical freeze-out without 
undergoing interactions. 



We can, however, go further and use the suppression of the considered 
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resonances as a tool capable of estimating conditions at particle freeze-out. We 
consider the decay of a generic resonance N* 



in a gas of pions and nucleons. We shall assume that one interaction of either Y or tt 
is sufficient for that resonance to be undetectable, and that the decay products travel 
through the medium with speed Vi (where i can mean either y or tt) . The interaction 
probability is proportional to Vi, the interaction cross-section of the decay product 
with each particle in the hadronic medium ((Tij(t'i), where j can refer to either pions, 
Kaons, nucleons or antinucleons. Note that the cross-section itself depends, in a 
generally complicated way, on the incident momentum, and hence on Vi), and the 
particle density in the fireball pj. pj is increased by a factor 7i = 1/ -^1 — due to 
Lorentz-contr action, and decreases as time passes because of the fireball's collective 
expansion (parametrized by the flow velocity Vfiow, assumed to be of the order of the 
relativistic sound speed c/ \/3.) The time dependence of the densities will thereflore 



and poj, the density of j at hadronization, can be calculated from the chemical 
freezeout temperatures and chemical potentials. Putting everything together, the 
rescattering reaction rate is 



N* ^Ytt, 



(5.2) 



be 




(5.3) 




[o-m(Vi)p07r + 0-iK(Vi)poK + 0-iN(Vi)poN + CriN(Vi)PoN)] (7v)i 




(5.4) 



If we use the average. 




(5.5) 
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(where p and m are the resonance's momentum and mass) eq. (5.4) becomes, 

Pi = [{<7in) POtt + (cTiK) POK + {(TiN) PON + {(7^) Pon)] ^ ( ^ t] ' i^'^) 

rrii \K + Vflowt / 

Neglecting in-medium resonance regeneration and particle escape from the fireball, 
the population equation describing the scattering loss abundance (A^^i) is: 

^ = -N^,-N,P,, i^l,2 (5.7) 
at T 

= --N,. (6.8) 

at T 

The required nucleon and antinucleon density pn is obtained through Eq. (??) 

We consider the nucleons to have a mass of ~1 GeV, and a degeneracy of 6, to take 
the p,n and the thermally suppressed but higher degeneracy A contributions into 
account, the pion density is computed in the massless particle limit, leading to 

^- = E(^4-Jt,)*«'^A-.(^) (5.10) 

The model presented here is remarkably insensitive to the individual cross-sectional 
areas. The values we used in the calculation are given in table 5.2,but order-of- 
magnitude variations of the more uncertain cross-sections did not produce varia- 
tions of more than 30%. Similarly, the value of the initial fireball radius Rq (which is 
constrained by the entropy per baryon) does not significantly affect the final ratios. 
This reassures us that had we used a more exact approach than the approxima- 
tions in Eq. (5.5), the qualitative features of our model would not have changed. 
The results, however, exhibit a very strong dependence on both the Temperature 
(which fixes the initial resonance yield as well as the hadron density of the fire- 
ball) and fireball lifetime (in a short-lived fireball not many resonances decay, so 
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Figure 5.5: Produced (dashed line) and observable (solid lines) ratios A(1520) / (total 
A),E*/(total A), K*^/K and S/(total A). The solid lines correspond to evolution 
after chemical freeze-out of 1,2,3,4,5,7,10,15,20 fm/c, respectively. The values at 
time zero (chemical freezeout) were taken from Fig. 5.4. See Fig. 5.4 caption for 
the meaning of K, K* . 

their products do not get a chance to rescatter). Fig. 5.5 shows the dependence 
of the A(1520)/A, S7A and K*^{%92)/K on the temperature and lifetime of the 
interacting phase. It is clear that, given a determination of the respective signals 
to a reasonable precision, a qualitative distinction between the high temperature 
chemical freeze-out scenario followed by a rescattering phase and the low temper- 
ature sudden hadronization scenario can be made. We also note that despite the 
shorter lifetime of the S* and higher pion interaction cross section, more E* decay 
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products should be reconstructible than in the A(1520) case, at all but the highest 
temperatures under consideration. This reinforces our proposal that the S* is a 
very good candidate for further measurement. 

Diagrams such as those in Fig. 5.5 still contain an ambiguity between tem- 
perature and lifetime of the interacting hadron gas phase. A low observed ratio 
can either mean a low freeze-out temperature or a lot of rescattering in a long 
re-interaction phase. However, this ambiguity can be resolved by looking at a se- 
lection of resonances, with different masses and lifetimes. Fig. 5.6 shows how the 
initial temperature and the lifetime of the re-interaction phase decouple when two 
resonance ratios are measured simultaneusly. A data point on diagrams such as 
those in Fig. 5.6 is enough to measure both the hadronization temperature and to 
distinguish between the sudden freeze-out scenario and a long re-interaction phase. 

The plots in Fig. 5.6 can also be used as consistency checks for the model: 
For example, the near independence of E*/S on temperature means that the equal 
temperature lines in the E*/S vs E*/A diagram are nearly insensitive to the details 
of the rescattering model. Moreover the mass differences and lifetimes of the E*, K* 
combine in such a way as to make the E*/ (all A) vs i^*°(892)/(all K-) diagram fold 
into a very narrow band. Any serious shortcoming within our rescattering model 
would be revealed if the observed particle ratios stray from this band. 

To understand the low A(1520) multiphcity measured by NA49 and STAR 
it should be realized that the A(1520) is a very unusual particle (Watson et al. 1963). 
Unlike most other hadronic resonances (E*, K* , A, p ecc), its high spin is due not to 
valence quark spin configuration but to the fact that the A(1520) s valence quarks 
are (using the constituent quark picture) in an L = 1~ state. The dominant decay 
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Table 5.2: Scattering model parameters 



c^ttAt (mb) axN 


Ctttt CTttK CTATTV O'lVAT 


24 20 


40 20 24 50 




rA(1520) r_fs'*0(892) 


35 MeV 


15.6 MeV 50 MeV 


escape rate (fm" 


negligible 


\' 


0.5 


R(fm) 


8 145/T [MeV] 




220 MeV 



(N K), however, has to go through a relative momentum h—2 (d-wave) process, 
through a channel which is very close to threshold production. All this conspires to 
reduce the A(1520)'s width: Isospin conservation reduces the number of channels 
it can decay (most notably A(1520) Att is not allowed), while the high relative 
angular momentum (and negative space parity) and threshold suppress the decay 
phase space for both the dominant [NK) and additional (Ett, Atttt) A(1520) decay 
modes. 

Therefore, measurement of the A(1520) presents both experimental chal- 
lenges (The high partial wave component introduces an angular/spin dependence 
of the decay products which is difficoult to measure in an unpolarized high multi- 
plicity experiment) and theoretical uncertainities (perhaps new in-medium physics 
is at work in L 7^ resonances (Markert et al. 2002b; Zimanyi and Levai 2004)). 
It should be noticed that, as Fig. 5.7 shows, a 50% suppression of the A(1520) sig- 
nal at hadronization would mean the data is perfectly compatible with the sudden 
freeze-out model described in (Rafelski and Letessier 2000). This strongly suggests 
that other resonances should be used to constrain the fireball freeze-out properties. 
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Figure 5.6: How temperature and fireball lifetime decouple when two resonances of 
different masses and widths are compared. A point on any of the above diagrams is 
potentially sufficient to fix both of these quantities. The experimental measurements 
are discussed in the first section. See Fig. 5.4 caption for the meaning of K* . 



The non-suppression of the K* makes it likely that other hadron resonances 
can also be detectable. The diagrams shown in the previous section make it clear 
that the S* should be abundantly produced, and it has many characteristics which 
would make it a logical next step in the study of resonances produced in heavy ion 
collisions. will also be a candidate for this kind of study. 

Focusing on these resonances will also make regeneration less likely. In 
particular, regeneration is less likely to affet the S* detectability: 95% of S* decay 
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Figure 5.7: K* / K vs A(1520)/A assuming half of the A(1520) s are suppressed 
in-medium. See Fig. 5.4 caption for the meaning of K* . 

through the p-wave S* Att channeL However, regenerating S* s in a gas of A s 
and vr s is considerably more difficoult, since Avr scattering will be dominated by the 
s-wave Att S^. For these reasons, a measurement of the E*/A and (which 
is also unlikely to regenerate due to its small width) should be able to distinguish 
a sudden hadronization model from a slow sequential freeze-out with certainity. 



5.4 Production of pentaquarks 

The discovery by the NA49 collaboration (Alt et al. 2004)of a new S~(1862) 
I = 3/2 narrow F < 5 MeV resonance in their pp background, (^snn = 17.2 
GeV) demonstrates that heavy ion detectors are capable of discovering states missed 
at previous energy scans. Statistical hadronization and strangeness enhancement 
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mean that these states should be produced abundantly in heavy ion collisions. The 
enhanced production rate, together with experimental capabilities, makes the Heavy 
ion experiments promising probes for the discovery of new hadronic states. 

This newly discovered hadron resonance has, given the mass and charge, 
an exceedingly narrow width. This feature is common with ©+(1540), another 
recently reported resonance (Nakano et al. 2003; Barmin et al. 2003; Stepanyan 
et al. 2003), which decays into the channel with quark content uudds and / = 0. 
This is beheved to be the predicted (Diakonov et al. 1997), lowest mass, pentaquark 
state. The 5*(1862) can be interpreted as its most massive isospin quartet member 
ssddu, ssudu, ssudd, ssuud with electrical charge varying, respectively, from —2 to 
+1, in units of |e|. 

Appearance of these new resonances can have many consequences in the 
field of heavy ion collisions. We at first explore how the introduction into the 
family of hadronic particles of these two new resonances. ©+(1540) and S*(1862), 
influence the results of statistical hadronization fit to relativistic heavy ion hadron 
production experimental results. We use the same data set as has been employed 
in Ref (Rafelski and Letessier 2003b) and obtain predictions of how the relative 
abundances of these new resonant states vary as function of the heavy ion collision 
energy. 

Importantly, only the two already identified states with / = 0, and 1 — 3/2 
of the anti decuplet, which also includes the / = 1/2, and 1 = 1 states are of 
relevance in the study of the statistical hadronization fits. Thus, in our analysis, 
we do not depend on the unknown masses of / = 1/2, and 7 = 1 states. However, 
the interpretation of the newly discovered narrow states as pentaquarks enters our 
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considerations decisively. In our approach (Rafelski and Letessier 2003b), as in other 
recent work (Becattini et al. 2003), the chemical equilibrium and non-equilibrium 
is considered. Accordingly, we allow quark pair phase space occupancies, for light 
quarks 7^ 7^ 1, and/or strange quarks 7s 7^ 1. 

The pentaquark (if this is indeed a pentaquark) valance quark content en- 
ters the assigned chemical fugacities and phase space occupancies in a different way 
than either mesons or baryons. Hence, predictions for the pentaquark yield are 
likely to be very different whether equilibrium is assumed or not. In particular, if 
the phase space is saturated above equilibrium, this will yield to a large (~ 7'^) 
enhancement of pentaquarks with respect to hadrons of a similar mass. And nature 
has been generous enough to provide us with detectable hadrons of a comparable 
mass and width: The A(1520) and S(1530), with mass and width nearly identical 

to e+. 

Since we study at SPS the total particle multiplicities, and at RHIC the 
central yields which can be considered produced by rapidity-localized fireballs of 
matter, we require in our fits balance in the strange and antistrange quark content. 

There are two independent fit parameters when we assume complete chemi- 
cal equilibrium, the chemical freeze-out temperature T and the light quark fugacity 
\ = V-^u-^d = e^''^^^'^\ The baryochemical potential /Xf, is the physical param- 
eter controlling baryon density. Strangeness conservation fixes the strange quark 
fugacity As (equivalently, strangeness chemical potential, for more details see, e.g., 
(Rafelski and Letessier 2003b)). Adding the possibility that the number of strange 
quark pairs is not in chemical equilibrium, 7s 7^ 1, we have 3 parameters, and al- 
lowing also that light quark pair number is not in chemical equilibrium, we have 4 
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parameters. These three alternatives will be coded as open triangles, open squares 
and filled squares, respectively, in all results we present graphically. 

We find that the new resonance ©+(1540) infiuences significantly the sta- 
tistical hadronization fit to particle production at the lowest SPS energies. In a 
baryon rich environment the introduction into the fit of ©"^(1540), a, b — 1 baryon 
with 'wrong' strangeness infiuences the strangeness balance condition, and thus in- 
directly the individual yields of all strange hadrons. This leads to a reduction in 
the statistical fit error for our hadronization study of the AOA GeV Pb-Pb reactions 
where we see a significant change in the relative yield of kaons and A. We also find 
changes in the details of the statistical fit parameters. In comparison to (Rafelski 
and Letessier 2003c), aside of the introduction of the new resonances, we also have 
harmonized our hadron decay table with those used by the Krakow group (Baran 
et al. 2004). The improvement of the particle yield fit is both, a theoretical confir- 
mation of the validity of the statistical hadronization model of particle production, 
and its apphcability at low SPS energies. We show how the fit error evolves in 
Fig. 5.8, which is also presented in the bottom lines of tables 2.1 and 2.2 along 
with the number of data points and resulting degrees of freedom. Considering the 
small number of degrees of freedom at SPS, we need x^/dof < 1 to have good sig- 
nificance of the fit^. The errors seen in Fig. 5.8 are, for the chemical nonequilibrium 
case (filled squares) , sufficiently small to allow us to conclude that the introduction 
of ©+(1540) assures that the statistical hadronization works well down to the low- 
est SPS energies. To compare with earlier results on x^/dof, obtained prior to the 

discovery of these new resonances, see Ref. (Rafelski and Letessier 2003c), figure 16. 
^See appendix B 
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Figure 5.8: (Color online) x^/dof for statistical hadronization fits at SPS and RHIC: 
results are shown for 40, 80, 158^4 GeV Pb on stationary Pb target collisions and 
at RHIC for 65+65 and lOO+lOOA GeV Au-Au head on interactions. 

An interesting point, seen in Fig. 5.8, is that the chemical equilibrium 
fit 7s = 1, 7g = 1 is rendered unacceptable at all SPS energies in presence of 
the new resonances. The semi-equilibrium fit, which allows a varying strangeness 
saturation, but assumes light quark equilibrium is generally resulting in twice as 
large compared to the full non-equilibrium approach. In a study of profile 
as function of 7^ we find a clear and strong minimum for 7^ ^max ^ gm^/(2T)_ 
Acquisition by the fit of this limiting value implies that there is no fitting error in 
the 7q presented below. 

The chemical freeze-out parameters of the fits considered play a very im- 
portant role in predicting the (relative) yield of hadronic particles, and this depen- 
dence is even stronger for many pentaquark states, due to their unusual quantum 
numbers. These fit parameters for RHIC are shown in table 2.2, and for SPS in 
table 2.1 along with the freeze-out temperature. We note that for the full chemical 
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Figure 5.9: (Color online) Yield of B^(1540) in relativistic heavy ion collisions, 
based on statistical hadronization fit to hadronization parameters at SPS and RHIC 
40, 80, 158^4 GeV Pb on stationary Pb target collisions and at RHIC for 65+65 
and 100+lOOA GeV Au-Au head on interactions. Relative yields with Ks-, A, and 
A(1520) are shown from bottom to top. 

non-equilibrium, the freeze-out temperature is found to be smaller than for semi- 
equilibrium case. This reduction is over-compensated in pentaquark yields by the 
significantly increased value of 7^. 

We now consider the relative yields of the new resonances in figures 5.9 
(left) and 5.9 (right). These yields vary strongly with collision energy for the case 
of G~*'(1540) in Fig. 5.9 (left), but are rather constant in Fig. 5.9 (right). Certainly 
our result differs greatly from expectations arising from an earlier study of the sta- 
tistical model production of the 0"*"(154O) resonance (Randrup 2003) where the 
decisive variation of the particle yield with chemical potentials was not explored. 
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Moreover, the hadron yields, presented in (Randrup 2003), did not include the con- 
tributions from decay of short lived hadron resonances. We checked that the relative 
particle yields shown in (Randrup 2003) for zero chemical potentials and varying 
temperature are mathematically correct, also as a cross check of our program. 

In Fig. 5.9 (left), we show (from top to bottom) the relative yields 
0+(154O)/A(152O), e+(1540)/A, e+{15A0)/Ks for chemical nonequilibrium (solid 
lines), semi-equilibrium (7^ = 1, dashed lines) and chemical equilibrium (dotted 
lines). The yields of A used here include 50% weak interaction cascade from S. 

The reason that the chemical nonequilibrium is leading to greater than equi- 
librium yields is that the lower hadronization temperature is overcompensated by 
the chemical factors, e.g., 6+(1540)/A(1520) ~ l/27q(Ag/As)^ ignoring the small 
mass difference. The factor 1/2 is due to the difference in spin degeneracy. The 
actually observed yield ratio ©+(1540)/A(1520) is probably going to be still greater: 
As explained in the previous section, A(1520) is seen at 50% of the expected sta- 
tistical hadronization yield in heavy ion collisions. As the ©"'"(1540) is an L=0 
resonance, we predict that its yield should more closely match the statistical model 
calculation. 

In Fig. 5.9 (left), we also recognize that the reason that there is such a 
significant impact at low SPS energies of ©+(1540) is that it is produced at the 
level of -MO-20% of A in fits at AOA GeV. This is due to the large prevaihng 
baryochemical density. Clearly, this is the environment in which one would want to 
study the properties of this new resonance in more detail. However, at all energies 
considered, we find that at chemical non-equilibrium the ©"^"(1540) is more abundant 
compared to A(1520) and thus this new resonance could become an important probe 
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of the hadronization dynamics, provided detection issues, such as different branching 
ratios and experimental acceptances for the analyzed decays, (A(1520) — > pK^ and 
— > pKs) are resolved. 

The observation of the pattern of relative yield of 6"'"(1540), seen in Fig. 5.9 
(left), would firmly confirm the 4-quark, one anti quark content of this state. 
Namely, were for example the ©"'"(1540) another tri-quark baryon state, the yield 
ratio with (strange) baryons would be quite flat as function of collision energy. We 
further note that the absolute magnitude of the relative yield, seen in Fig. 5.9 (left), 
will be of help in estabhshing the degree of chemical equihbration. 

In Fig. 5.9 (right), we show at the bottom the expected relative yield of the 
S {lSQ2)[ssddu\ relative to 'Er[ssd]. The S*(1862) adds at the percentile level to 
the yield of observed 5 and thus it is less influential in the statistical hadronization 
approach. The absence of variation of the relative yield with collision energy is due 
to cancellation of chemical factors. This relatively small relative yield at all collision 
energies here considered shows that indeed the pp environment, where it has been 
identified by the NA49 collaboration, is most suitable. The dotted hues, in Fig. 5.9 
(right), are visibly breaking the trend in some of the results, indicating that the 
large chemical equilibrium fit generates unreliable statistical model parameters. 

We also show, in Fig. 5.9 (right) on the top, the yield of the pentaquark state 
Ti{mQ)[sddu%jt\ which for purpose of this study is assumed at the mass indicated. 
Again due to cancellation of key chemical factors in ratios shown in Fig. 5.9 (right), 
both being proportional to ■jqXd/^u, the ratio is flat (except for the failed fit chemical 
equilibrium results). Considering that ~ and 7^ ~ 1.6, the magnitude of 
relative yields seen in Fig. 5.9 (right) is primarily due to the hadron mass, and 



155 



degeneracy. 

We have shown that inclusion of the pentaquark states in the study of parti- 
cle production in heavy ion collisions improves the quality of our fits to experimental 
data. We find that 6"'"(1540) state influences the low energy SPS particle yield fit 
results. It can be expected that it will be detectable, in particular at low heavy ion 
collision energies, and thus should become a new probe of hadronization dynamics. 
The other pentaquark states will be hard to observe in heavy ion collisions. 

5.5 Momentum dependence of the resonance-peirticle ratios as a freeze- 
out probe 

In section 5.4, we have shown that the measurement of resonances can probe both 
the hadronization temperature, and the lifetime of the interacting hadron gas phase 
(Torrieri and Rafelski 2002; Rafelski et al. 2001). Ratios of a generic resonance 
(henceforward called N*) to the light particle (which we will refer to as N) with 
an identical number of valence quarks are particularly sensitive probes of freezeout 
temperature because chemical dependence cancels out within the ratio. If we ex- 
amine this ratio within a given rriT > uin* range, we expect to disentangle flow and 
freeze-out conditions, since the ratio N*/N should not depend on mr for a purely 
static and thermal source (if there is no flow, N* oc e"*^ so {N* /N)mT>m* is 



We therefore take the most general Boost-invariant freeze-out hypersurface 
in the Boltzmann hmit (see Table 4.1, Eq. (4.11)) 



constant. 




(5.11) 
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T 



and use it to calculate the ratio between two particles with the same chemical 
composition. The chemical factors cancel out, and we are left with 



where g* and g refer to each particle's degeneracy and the function S{mT-,PT) is 
given by Eq. (5.12). (Note that rriT is the same for N* and N, but px varies). 

Fig. 5.10 shows the application of this procedure to the cases {K* + 
'K*)/{Ks) (top), E*(1385)/A (middle), and r}'/r] (bottom) at two freeze-out temper- 
atures and flows: T = lAQMeV ,Vmax/ c = 0.55 on left and T = 1 70 Me V, fmaa;/c = 
0.3 on the right. Significant deviations from simple constant values are observed, 
showing the sensitivity of the ratio to freeze-out geometry and dynamics. The an- 
alytically simple result in Eq. (5.14) is valid only if the light particle Y has been 
corrected for feed down from resonances, including N* . In other words, Eq. (5.14) as 
well as Fig. 5.10 require that decay products from reconstructed N* do not appear 
on the bottom of the ratio. Experiments usually do not do such feed down correc- 
tions (Fachini 2002; Markert et al. 2002a; Friese 2002), since this would increase 
both statistical and systematic error on the ratio, and it is not always possible to 
do such corrections at all (undetected decays) or in the full range of experimental 
sensitivity. 




(5.14) 
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Introducing the feed down corrections into Eq. (5.14), we obtain 

Kbserved ^ g*S{mT,P*T) .^ -^^x 

^observed gSimr, Pt) + Yli 9ibN*^NR{mT, PTi)' 

Here, S{mT,PT) describes the directly produced particles and has the form given 
by Eq. (5.11) and each term i?(mT,Pri) describes a feed down contribution, in the 
form given in section 2.3.1 (In particular Eq. (2.34) and following) 

Fig. 5.11 shows the ratios, including feed down of resonances, for the same 
particles and statistical hadronization conditions as were studied in Fig. 5.10. In 
the S*/(all A) case we omitted the feed down from H to A which is usually corrected 
for (if this is not done the ratio E*/(all A) would depend strongly on the chemical 
potentials). We did allow for the (f) — > KsKl feed down, since it is a strong decay 
that cannot so easily be corrected for^. We note that the feed down from particles 
with a different chemical composition cannot always be corrected for, and thus some 
resonances ratios will also acquire a (mild) dependence on the chemical potentials. 
This is even true for ratios such as r/'/(all r/), given different ss content. In this 
chapter, these type chemical corrections were set equal to unity. 

To further study the sensitivity of resonance-particle mr-ratio to freeze- 
out dynamics, we also present the (feed down corrected) case as a function of pt 
rather than m-r in Fig. 5.12. Unsurprisingly, we see grossly different behaviors, 
with many of the results coalescing. This of course is an expression of the fact 
that N* and N have dramatically different pt at the same rriT and vice versa. We 
believe that the rriT ratio will in general be more sensitive to freeze-out dynamics, 
since its dependence on itlt is dominantly due to freeze-out geometry and dynamics. 

However, thepr dependence seen in Fig. 5.12 provides an important self-consistency 
^See footnote in subsectfon 5.2.1 
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check for our previous results. We have found that the rriT ratios are often greater 
than unity even though there must be more ground state particles than resonances. 
Now it can be seen in the px ratio, that this requirement is satisfied. 

5.6 Discussion 

In general the the mr and pr dependence of the ratios in Fig 3 and, respectively. 
Fig 4 depends on freeze-out geometry, temperature and flow velocity. The intro- 
duction of a steeper flow proflle will further raise all of the considered ratios, since 
a considerable fraction of particles will be produced in regions that do not flow as 
much. The influence of freeze-out dynamics will generally go in the same direction 
as freeze-out approaches the explosive limit {dtf/dr —>■ 1). However, both the mag- 
nitude and the qualitative features of the two effects (flow and freeze-out velocity) 
will be considerably different. Especially, when more than one ratio is measured, 
it would appear that we will be able to determine the freeze-out condition. This is 
in contrast to the mr distributions in Fig. 4.2, where the effects discussed in this 
paper result in linear corrections, which tend to compete, making the task of ex- 
tracting the freeze-out dynamics much more ambiguous. Thus, there is considerable 
potential of resonance-particle mr-ratios as a freeze-out probe. 

The presence of a long living hadronic gas rescattering phase can distort 
our freeze-out probe. In particular, the apparent N*/N ratio will be altered due to 
the depletion of the detectable resonances through the rescattering of their decay 
products. Its dependence on mr will be affected in a non-trivial way, since faster 
(higher px) resonances will have a greater chance to escape the flreball without de- 
caying, thus avoiding the rescattering phase altogether. Regeneration of resonances 
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in hadron scattering may add another dependence which is different for the 
E*/A and the K* / K ratios (Bleicher and Stocker 2004). While, as discussed in 
the last section, evidence that a rescattering phase plays a great role in particle 
distributions is lacking, it would seem that the "safest" probes for freeze-out are 
the particles and resonances most unlikely to rescatter. 

For this reason we have included the r/'/?] ratio in our considerations. 
77 ^ 77 and 77' 77 have very different branching ratios, but have the same 
degeneracies and similar but rather small partial widths. The electromagnetic de- 
cay mode is practically insensitive to posthadronization dynamics. Regeneration 
effects are suppressed since the hadronic two body decay channel is suppressed. 
All these features make these particles interesting probes, allowing for the analysis 
considered here. 77, 77' mesons have been measured at SPS energies in the 77 decay 
channel (Albrecht et al. 1995; Lebedev et al. 1994), and detectors such as PHENIX 
are capable of reconstructing the same decays at RHIC. 

In summary, we have presented an overview of the different statistical 
freeze-out models used to fit heavy ion data. We have shown how the freeze-out 
geometry and freeze-out dynamics influences the hadron spectra. Our primary re- 
sult is the finding that the rriT dependence of the resonance-particle ratios is a 
probe of freeze-out. We have presented these ratios for three particle species and 
two freeze-out conditions and have considered how our results could be altered by 
posthadronization phenomena. 
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Figure 5.10: (Color online) Dependence of the K*/K, S*/A and rj'/eta on the 
Freeze-out model . 
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Figure 5.11: (Color online) [K* + ir*)/(all Ks), S*(1385)/(all A) and V/(all r/) 
ratios, including feed down from resonances. 
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Figure 5.12: (Color online) Vt dependence of (K* +;^)/(all Kq), S*(1385)/(all A) 
and vl I (all r\) ratios, including feed down from resonances. 
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CHAPTER 6 



Azimuthal anisotropy 



6.1 Introduction 

Matter flow azimuthal anisotropy has long been considered a promising soft hadron 
observable in the study of heay^ ion collisions. This anisotropy is considered an 
important evidence for collective matter flow (OUitrault 1992; Voloshin and Zhang 
1996), as it indicates thermalization early in the system's evolution (Sorge 1997). 

Generally one simplifles the complex transport problem and considers hy- 
drodynamic model, i.e. dynamics of flrst moment of the momentum distribution 
of particles. Assuming in the local reference frame an asymmetric thermal mo- 
mentum distribution we evaluate below the actual flnal state particle distributions. 
Azimuthal anisotropy in hydrodynamics arises from anisotropic gradient of pressure 
(force fleld) when the collision between nuclei is not exactly head-on central. In this 
way hydrodynamic model leads to the azimuthal momentum dependent asymmetry 
of flnal state particle distributions. 

Anisotropy in the flnal state can be quantifled in terms of a Fourier decom- 
position of the momentum distribution (Poskanzer and Voloshin 1998), 



where the angle is deflned with respect to the reaction plane (OUitrault 1993) in 




(6.1) 
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an event-by-event analysis (Fig. 6.1). 

The second coefficient V2, usually called elliptic flow, has been subject to a 
considerable amount of experimental investigation at AGS, SPS and RHIC energies 
(Pinkenburg et al. 1999; Appelshauser et al. 1998; Adler et al. 2002a; Adler 
et al. 2003). The fourth coefficient V4 has recently been measured at RHIC (Adams 
et al. 2004). (odd components, such as the directed flow vi, disappear in the 
Boost-Invariant limit) 

Generally, Vn is a function of momentum. It has been found that V2{pt) 
at low pt and central rapidity agrees with hydrodynamical predictions (Kolb and 
Heinz 2003; Hirano and Tsuda 2003). However, RHIC experiments have shown that 
at ~ 1-5 GeV V2 rise saturates to a particle-dependent limit (Adler et al. 2002a; 
Adler et al. 2003), something so far hydrodynamic approach could not explain Since 
this Pt value is also seen as being too soft for V2{pt) to be determined by "hard" 
perturbative processes, this behavior has not been fully understood. 

However, coalescence of partons obeying perturbative dynamics has been 
shown to lead to saturation of t'2(pr) (Molnar and Voloshin 2003). While this 
approach does not as yet cover consistently the low pt range (Pries et al. 2003) and 
hence is yet to undergo detailed quantitative testing, its success suggests that, more 
generally, the process of quark-gluon hadronization can influence significantly the 
resulting V2{pt)- Here we address this possibility systematically, but qualitatively in 
that we do not introduce resonance decays and flow proflle of matter. Our objective 
is to identify physical mechanisms rather than to explain the experimental data. 

Our study will show that there are V2{pt) saturation mechanisms based 
solely on rapid hadronization of a hydro-dynamically evolving opaque system, such 
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as quark- gluon plasma. We shall use the Cooper- Prye approach of Eq. (2.17) 

(Cooper and Frye 1974). 



As we have argued in the previous chapters, this formula is particularly relevant 
if most observed particles are emitted when the system undergoes a phase transi- 
tion to a gas of particles having a much larger mean free path (QGP — > expanding 
Hadron Gas). In this case, the Cooper- Frye formula, in particular the hadronization 
hypersurface ci^S^, ceases to be just a computational prescription to generate par- 
ticles from a continuum, but acquires physical significance as the representation of 
the hadronizing QGP in position space. Thus a more thorough exploration of how 
it can modulate the shape of V2{pt) is required . The form of ci^E^ is determined 
by the dynamics of the QGP — > HG phase transition. 

Several theoretical studies and also general behavior of experimental results 
suggest that hadronization happens rapidly, perhaps through viscous fingering of 
the vacuum (Csorgo and Csernai 1994), driven by a mechanical instability at the 
point when the pressures (including that of flow) of the two vacua balance , but 
the velocity still point outwards (Rafelski and Letessier 2000). In this situation it is 
natural to expect a "burning log" type emission (Torrieri and Rafelski 2003b), with 
a fast (~ c) moving emission surface rapidly consuming the fireball from outside 
in. In contrast to this dynamical hadronization picture, a "blast-wave" freeze-out 
model in which the freeze-out time does not depend on the transverse freeze-out 
radius is presently popular among experimentalists, perhaps due to the simplicity 
of its use in a fitting procedure (Burward-Hoy 2003; Adams et al. 2003b; Afanasiev 
et al. 2003; Sandor et al. 2004). 




(6.2) 
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Particle emission 



Anisotropic expansion 



Figure 6.1: A general freeze-out surface in configuration space . 

As we will show these two scenarios (burning log and non-dynamic blast 
wave) give considerably different V2{pt) for the same initial matter flow anisotropy. 

In particular, we find that a burning-log type freeze-out dynamics leads in qualita- 
tive terms to the observed V2 saturation. 

6.2 Vn in the Cooper-Frye approach 



We shall use the approach taken in chapter 4, using boost-invariance to construct 
our flow and E^. Hence, in the most general boost- invariant case 

/ tfCOsh{yL) \ ( cosh(yi) \ 



Tf cos{9) 
Tf sin(^) 
V ^/sinh(?/L) ) 



vtCos{9) 
i'Tsin(^) 
y sinh(?/L) ) 



(6.3) 



Where 9 is the angle of the emission point in configuration space ( Fig. 6.1, not to be 
confused with the emission angle 0). However, we forgo cylindrical symmetry, and 
Fourier expand both the fiow and hadronization hypersurface, keeping only even 
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terms given the symmetry, 

oo 

tf(r^^) = J]A2nt/(r) cos{2ne), 

n=0 

oo 

'^/l^)^) = A2nr/(r) cos(2n^), 

n=0 

oo 

VT{r, 9) = A2nv{r) cos{2n9). 



(6.4) 
(6.5) 
(6.6) 



n=0 



We assume a purely elliptic fireball (n = 0, 1 only), and moreover a "fast" freeze- 
out so as to neglect tf ^-dependence. If the fireball is approximately circular in 
transverse space, the elliptical freeze-out hypersurface becomes an ellipse (with a 
correction of 0(5r^)), 



It = 



( tf{r)cosh{yL) \ 
(r + 6r) cos{9) 
rsm{9) 
V *f(^) sinh(yL) J 

( cosh(|/L) ^ 
[v -t- bv cos(6')) cos(6') 
[v -|- bv cos(^)) sin(^) 
y sinh(yi) ) 
1 

- (^; + 5^;cos(e))2■ 



(6.7) 



(6.8) 



(6.9) 
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We can now parametrize the freeze-out hypersurface element in terms oi r,9,yL — 
(radius, angle and longitudinal rapidity yi) 

dr 89 dyi 

^ cos]i{yL){r + 5rs\v?{9)) ^ 
cos{9) 

^{r + 5r)sm(9) 

\^ sinh(yL)(r + 5r sin^(6')) J 



dyLdrd9. 



(6.10) 



If we combine the obtained u'^ and ci^E^ with the usual parametrization for the 



emitted particle's 4-momentum 



rri'X' cos 



(6.11) 



Pt sin(0) 
Pt cos(0) 
y mTsinh(y) ) 

we obtain an expression of the boost-invariant momentum distribution in terms of 

4>,PT,mT and y: 

E% = J^""^' ^dr JZ^ dy, ^ dO K cosh(y - y,) x 
X (1 + f sin2(^)) -p^^(cos(^ - 0) + f^sin(^))" 

/(^[mT cosh(y — yi) — pt{v + Svcos{29)) cos{9 — 0)]) 
Where f{x) — in the Boltzmann approximation (good for all particles except 
pions) and for pions: 

oo 

/(^) = J](-l)"Ar^e-"^ 



(6.12) 



n=0 



Each Vn can now be calculated using 
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the integrals over and yi can be done analytically in the boost-invariant limit, 
using the modified Bessel Functions 

Jo 

/oo 
dvL cosh(n|/L)e-^"°^'^(2/i) 
■oo 

and the following result 



2n 
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^sin(n0)e»^°^('^)+^^^°('^) = 
Jo 27r 



where T„(x) and Sn{x) are defined in terms of Chebyshev polynomials 

Tn{x) — cos(ncos~"^(x)) 

Sn{x) = sin(ncos~^(a;)) = ^/l - T^{x) 



(6.14) 
(6.15) 



(6.16) 
(6.17) 



Putting everything together, we get 



VniPr) = 



JrdrJddJn 
2jrdrJdeJo' 



where (a = Pt{v + 6v cos(6'))) 



>o 



i^i(2^)(i + ^sin(^)2) _p^|firo(2i^) T„(cos[^])4(a) -p^^^iro(^) 



(6.18) 



dt, 



[I cos[9] (T„+i(cos[6'])/„+i(a) + T„_i(cos[6'])/„_i(a)) ■ 
isin(6') {Sn+i{cos[9])In+i{a) - Sn-i{cos[9])In-i{a)) ] 



(6.19) 



and 



Xi(^)(l + ^sin(^)2) -p^^JlKoi'-xmr^ Io{a) - pt^Xo(^)/o(«) 



dtf 



(6.20) 
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For a quantitative fit to be physically meaningful, resonance contributions (~ 50% 
of the total Ks,A., more for tt (Torrieri and Rafelski 2003a)) need to be taken into 
account. As we have seen in the previous chapters, calculating these contributions is 
numerically intensive (Baran et al. 2004), and work to perform such a quantitative 
fit, together with analysis of particle spectra (Torrieri et al. 2004), is currently in 
progress. 

6.3 Results and discussion 

We have used Eq. (6.18) to quahtatively explore the dependence of V2{pt) for Ks, A 
and TT on the parameters of statistical freeze-out. It is apparent that the freeze-out 
dynamics parameter dtf/dr has a non- negligible effect on V2{pt) (Fig- 6.2 right). 
In particular, by changing the freeze-out hypersurface towards the "burning log" 
{dtf/dr = —1) limit, it is possible to raise and lower the saturation scale and pr 
at which the saturation occurs. Varying dtf/dr also controls the mass dependence 
of the observed i'2(pt), with an increased dtf/dr leading to a different saturation 
scare for the and A, as observed by experiment. Comparing the right panels of 
Fig. 6.2 left and right, it becomes clear that a simultaneous fit for the V2{pt) of 
several different particles would constrain both the flow anisotropy and the freeze- 
out dynamics. 

Of course, a physical fireball is not characterized by a single flow, but by 
a spatially varying flow profile. Fig. 6.3 (right) explores the effect of a range 
of physically reasonable profiles on the observed V2- Unsurprisingly, adding a flow 
profile diminishes the observed V2, while maintaining the saturation scale determined 
by dv/v and dtf/dr. 
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Finally, we note that corrections due to azimuthal anisotropy of the freeze- 
out hypersurface and chemical non-equilibrium are also very small, even as the 
chemical potential approaches the Bose- Einstein condensation hmit.(Fig. 6.3 left. 
Fig. 6.4 left). Temperature corrections are more noticeable (Fig. 6.3 left). 

Eq. (6.19) makes it apparent that the interference between space and flow 
anisotropy means that even a purely elliptic fireball will acquire f „>2 7^ for all even 
numbers Fig. 6.4 (right) shows the contributions of v„ arising from the freeze-out 
of a purely elliptical fireball. It can be seen that, for n = 4, this contribution is not 
negligible, and is in fact comparable to that arising from Hydro (Kolb 2003). 

This result, together with the recent experimental observation of directed 
flow (Adams et al. 2004), suggests that and away from central rapidity could 
also be promising observables: Our model finds 1)3,5 = since we assume boost- 
invariance. However, a full 3-fluid calculation, together with Eq. (6.19), would yield 
these quantities easily. Given the significant observed directed fiow (Adams et al. 
2004), as well as the V4 calculated in Fig. 6.4, odd Vn>i s could be within experimen- 
tal observation. These coefficients make for particularly interesting probes, since 
they arise solely from freeze-out effects and deviations from boost-invariance. 

In conclusion, we have shown that observed V2{pt) depends strongly on 
the freeze-out hypersurface, and that a burning-log type freeze-out can explain the 
observed i'2(pt) saturation. We have further shown that i'n(pr) is an extremely 
sensitive freeze-out probe, capable of constraining statistical freeze-out models very 
tightly. As experimental data in this area becomes richer and more precise, more 
realistic numerical studies of both spectra and Vnipr) might establish themselves as 
crucial analysis tools. 
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Figure 6.2: Left: Varying flow anisotropy. Right: Varying freeze-out dynamics . 
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Figure 6.3: Left: Varying geometrical anisotropy. Right: Varying flow spread . 
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Figure 6.4: Left: Varying flow profile. Right: Higher Vn s from an elliptic fireball 
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CHAPTER 7 
Summary and outlook 

7.1 What chapters 3-6 have shown 

The statistical hadronization model, together with hydrodynamic flow, are capable 
of describing a large amount and variety of experimental data with comparatively 
few fit parameters. However, the variety of statistical models on "the market", 
which manage to describe the data acceptably (to at least part of the heavy ion 
community) somewhat invalidates statistical hadronization as a physical picture. 

A "slow" staged equihbrium freeze-out starting at T=170 MeV (Heinz and 
Jacob 2000), A sudden equilibrium freeze-out at T=165 MeV (Baran et al. 2004), 
and a fast explosive freeze-out starting at T=140 MeV (Rafelski and Letessier 
2003b) are not "the same model" . They are profoundly different physical pictures, 
which lead to different conclusions about what the QGP-HG phase transition should 
look like. 

That all three have been successfully employed to describe experimental 
data simply means that two (or possibly all three!) of these pictures are the heavy 
ion equivalent of epicycles: Models with no connection with reality, which are lucky 
enough to fit some of the experimental data. 

In chapter 2, we have presented several candidates for a probe capable of 
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distinguishing the epicycles from the real physics Chapters 3-6 examined each of 
the candidates in detail, in light of the experimental evidence. To what extent have 
we succeeded at isolating the true physical picture? 

7.1.1 Does non-equilibrium sudden hadronization work? 

Sudden hadronization scenarios are considerably easier to falsify than a complicated 
staged freeze-out phase. In both the equilibrium and non-equilibrium scenarios, par- 
ticle yields place strong constraints on temperature and chemical potential. Particle 
spectra place equally strong constraints on temperature and flow. If the temper- 
ature which fixes the yields ("chemical freeze-out") differs from the temperature 
which fixes the spectra ("thermal freeze-out"), a fit to normalized yields should 
conclusively fail. 

That such a fit does not fail is, perhaps, the most firm conclusion to be 
drawn from this thesis. In chapters 3 and 4, we have shown that fits based on the 
same temperature for yields and spectra succeed , or at least are acceptable, both 
at SPS and RHIC energies, if the non-equilibrium sudden hadronization ansatz is 
used. 

If, however, hadronization happens in equilibrium, we find that the x^/DoF 
increases beyond meaningfulness, a conclusion enforced by the lack of definite min- 
ima in equilibrium profiles. Fig. 5.3 shows that the equilibrium scenario fails 
to describe short-lived resonance yields (which it does not aim to describe any- 
ways) : The scarcity of these resonances drives the fit to a temperature at which its 
impossible to describe ratios such as A/S without introducing 7^ and 7^. 
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The checks performed in this thesis are not, however, the only ones possible. 
In particular, in the following two years experiments should publish normalized 
spectra of unstable resonances. These spectra pose much tighter constraints because 
such resonances have the same quark composition as light particles. Hence, their 
normalization, with respect to these light particles, is controlled by temperature 
alone. 

Resonances and absolute normalization remove the correlation between 
temperature and flow: It is not possible to adjust the slopes by shifting the temper- 
ature and transverse flow along the same contour, since in this case normalization 
suffers. However, temperature is still strongly correlated with chemical potentials 
(especially 7^ and 7^). Using spectra of short-lived resonances will remove this cor- 
relation, and show once and for all if there is a meaningful temperature minimum. 

It should be noted that spectra of cf) and K* were already modeled at RHIC, 
and are well described by the single freeze-out ansatz. Much more, however, is yet 
to come. 

7.1.2 Have we ruled out staged hadronization? 

Ruling out Staged hadronization is considerably more difficult, since this model is 
not very quantitative beyond the "Two temperatures, one for chemical and one for 
thermal freeze-out" requirement. 

One can argue that a good description of both abundancies and spectra 
of particles with very different interaction channels within the same temperature 
and flow is, by itself, a piece of evidence against staged hadronization. If staged 
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hadronization occurs, why should the Q, with no s-channel interaction with the tt, 
have the same freeze-out temperature and flow as the tt, which interacts via the 
s-channel with pretty much every other particle in the hadron gas. 

Why should the A and A, or H and S be described by the same temperature 
and flow at the SPS, given that in a baryon-rich medium particles and antiparticles 
have very different reaction channels (Ap Ap. Ap nK, ...). This is even more 
true in the case of spectra of short lived resonances: As shown in Fig. 4.8, micro- 
scopic models predict a very px dependent re-scattering of decay products, which 
should bring the spectra of short-lived resonances away from hydro-dynamically 
predicted values. 

However, this evidence is not really enough to rule out the model, consid- 
ering that, as proven in chapter 4, a fit for quite a large range of particles can be 
made with either a "sudden freeze-out -|-resonances" ansatz, or a "no resonances and 
arbitrary normalization" description. The ^^'^ the fitted temperature comes out 
remarkably similar in these two approaches. 

Direct detection has made clear that the presence of resonances in spectra 
is an experimental fact rather than a model assumption. The analysis of resonance 
yields performed in chapter 5 points to a short interacting hadron gas phase, as 
well as a low chemical freeze-out temperature. Of course, proponents of staged 
hadronization would argue that the effect of out-of-equilibrium regeneration would 
invahdate this conclusion (Bleicher and Stocker 2004). 
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7.1.3 Have we constrained freeze-out dynamics? 

Here, the answer must be "not yet": Our fits to spectra have shown that freeze- 
out dynamics, parametrized for central and near-central collisions by is an 
important factor in determining inverse slopes, and can be extracted if enough 
particles are fitted. We have presented attempts at such extraction at both SPS 
and RHIC energies. A minimum in the profile was found, but it is not fully 
convincing. One way of constraining this minimum further is through fitting more 
low Pt particles, since including the low region makes the fit considerably more 
sensitive to the variation of ^ (The phase space volume of the emitted hadron 
~ E —pT^)- PHOBOS recently lowered the observed pr range by two orders of 
magnitude (Back et al. 2004), and more particles (certainly K, hopefully protons) 
will be coming. 

However, a deeper problem remains: As we have shown, the choice of fiow 
profile correlates with the choice of freeze-out dynamics. This is not surprising, 
since both quantities are determined by the condition chosen for freeze-out, as 
can be seen by comparing different hydrodynamic models (Teaney et al. 2001; 
Schnedermann and Heinz 1993). However, it means that extracting conclusively 
E'* from fits becomes considerably more computationally intensive, which makes 
exploring other probes, more capable of extricating freeze-out dynamics separately 
from fiow profile, necessary. 

Some probes were proposed in chapters 5 and 6. Resonances ratios as a 
function of are a promising analysis tool, but experimental results have yet to 
be presented. V2, described in chapter 7, is also a promising probe, but once again 
resonance contributions are needed to describe existing experimental data (Dong 
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et al. 2004). 

Finally, there is the traditional probe for exploring freeze-out in configura- 
tion space: Hanbury-Brown-Twiss interferometry (HBT) (Wiedemann and Heinz 
1999). Since this thesis did not include HBT calculations, we have not discussed the 
current state of affairs in any detail beyond the problems faced by statistical models 
in fitting HBT (Heinz and Kolb 2002). The solution to these problems might be 
tied to sudden freeze-out (Csernai et al. 2003), or to the fact that assumptions 
going into HBT models might be, for many reasons, overly simplistic (Gastineau 
and Aichehn 2000; Bolz et al. 1993; Wong 2004). 

7.2 Outlook on the heavy ion experimental program 

This thesis was written well before all of the RHIC data has been analyzed. Indeed, 
the most interesting and surprising experimental results might well still be in the 
future. As noted elsewhere, while the resonance program has collected a convincing 
signal for a huge array of particles (with more to come), detailed quantitative obser- 
vations, such as yields, ratios and spectra for all detected resonances and exotica, 
are yet to be published. 

Charm and bottom production is another topic which awaits most of its 
experimental data. Heavy quarks will probably never be chemically equilibrated at 
any energy, and the great majority of their yields will always be produced during 
the initial high Q"^ collisions (Thews et al. 2001; Schroedter et al. 2000). How- 
ever, a quark gluon plasma will equilibrate them thermally, and hadronization will 
distribute them in a radically different way than microscopic hadronic processes. 
For sudden freeze-out to be consistent, it should describe quarks independently of 
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their mass. This means that charmed particles {t/j^t/j'^D etc), beauty states and 
mixtures (Be) should be produced according to their statistical weights, and their 
momentum distribution should be modeled by the same flow as that of light quarks. 
We shall see if this is really the the predictions of statistical models can 

differ drastically (~ order of magnitude) from the purely hadronic ones (Schroedter 
et al. 2000). 

Beyond the current RHIC program, we can look forward to the TeV energy 
range accessible to the LHC. While hadronization conditions will probably not be 
that different from RHIC (the chemical potentials at mid-rapidity should both be 
close to zero), the initial conditions and equilibrated evolution are likely to be 
markedly different. In particular, initial collisions are likely to generate a number 
of heavy quarks comparable to the strange quarks created at SPS and RHIC. If 
the system thermalizes as rapidly as the RHIC system seems to be, its temperature 
during the initial stage of the hydrodynamic evolution is likely to be well in the 
perturbative QGP range. Perhaps, if this phase lasts a long time, we will also see 
interesting electroweak effects (another phase transition accessible to experimental 
study?) . At present it is not clear what the effect of this initial evolution will be, 
and whether any hadronic signals will survive subsequent cooling. 

The most interesting physics, however, might not come so much from push- 
ing the energy boundary higher and higher, as from comparing as many different 
collision energies, system sizes and chemical conditions as possible. The SIS, AGS 
and SPS have explored conditions in a wide range of nuclei at energies ranging from 
1 to 20 GcV. RHIC has pushed this boundary to the hundreds of GeV, the LHC 
will push it into the TeV range, and the future SIS facility will look for low energy 
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but high-density systems. Taking all these experiments together results in a large 
region in temperature, chemical potential and system size which needs to be sys- 
tematically analyzed. At what point do qualitatively new features (eg equilibration, 
or 7q^s 7^ 1 ) come in? What is the critical system size or collision energy needed for 
"most of the system" to become thermalized? Where exactly is the phase transition 
boundary, and are there any critical points? 

We do not yet know the answer, but can say convincingly that soft physics 
will play an important role in obtaining it. As remarked in chapter 1, hard probes 
(such as jet quenching) are useful to tell us some bulk properties of the matter 
produced in highly energetic heavy ion collisions. If, however, one wants to explore 
how these bulk properties change with energy, or even relate the measured values to 
thermodynamics, hard processes, on their own become much less useful. It is only 
soft physics which will be able to determine at what point a qualitative change in 
microscopic degrees of freedom has occurred. To accomplish this, however, requires 
more than the new experimental data. Statistical hadronization itself needs to be 
understood better, both at the phenomenological and fundamental level. 

7.3 Outlook on statistical hadronization 

As the previous sections have shown, there is a lot more data which need to be 
analyzed in the statistical model framework to be able to tell convincingly which 
statistical model reflects physics. Fortunately, as the previous chapters have men- 
tioned, the soft physics data taking at both SPS and RHIC is far from over. The 
range of particles analyzed in this theses, and many more, need to be systemati- 
cally examined at the large range of collision energies and system size accessible to 
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modern day and future experiments. 

For this reason, the immediate next step as far as the work described in 
this thesis is concerned is the development of a statistical hadronization code to 
be made available open-source to experimental collaborations. This will provide 
experimentalists with an open and universal "standard statistical model" to which 
to compare all their data, as well as a guarantee against errors (which, as we found, 
are very easy to make, given the thousands of particles and decay modes which 
need to be calculated). We have been developing such a code (named SHARE) 
together with the Krakow group (Baran et al. 2004). The first part of SHARE, 
which calculates An yields, has recently been released (Torrieri et al. 2004). Further 
modules (covering spectra, V2 and HBT) are currently in development. 

Another possible direction for statistical hadronization research is inves- 
tigating its mathematical consistency, and improving its formalism. At present, 
there is no indication that cither sudden freeze-out or staged freeze-out are self- 
consistent. Would inelastic interactions of a hadron gas really freeze-out at T=170 
MeV, as staged freeze-out assumes? Would elastic interactions of an explosively 
hadronizing plasma at T=140 MeV really be negligible? If the answer to both 
these questions is no, to what extent does departure from equilibration impact 
observables? What is the quantitative systematic error such effects give to fitted 
thermal and flow parameters? 

One way to answer these questions is to feed statistical model output into 
microscopic models, such as uRQMD. At present, uRQMD is incapable of fully 
describing any of the RHIC observables (strange particle yields, collective flow, 
and resonance production). However, this does not preclude statistical production 
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followed by hydrodynamic evolution. The Monte-Carlo module described in chap- 
ter 3 is being developed and incorporated into SHARE, as a "plug" for uRQMD. 
Hopefully, this will lead to some of the above questions being answered. 

On a more fundamental level, it is clear that the statistical hadronization 
formalism needs developing. Reducing all particles to Lorentz-scalar objects emitted 
from a 3D locus in space, as the Cooper-Prye formula does, is a rough approxima- 
tion which might miss important physical features of the hadronizing system. In 
particular, interaction of the two phases of matter at hadronization might have non- 
negligible consequences (Bugaev 2004). Enforcing conservation laws and entropy 
non-decrease locally on each hadronizing volume element might lead to non-trivial 
changes of the statistical parameters across the phase transition. 7^ > 1, discussed 
extensively in this dissertation, is a simple example. It might not be the only param- 
eter which changes significantly as the system undergoes a rapid phase transition, 
(other workers have proposed a change in T (Keranen et al. 2003)). 

More ambitiously, the above questions might be answered by tying the 
statistical model to a more fundamental picture of hadronization, possibly involving 
a microscopic effective theory. Perhaps insights in this direction might be obtained 
by examining at what energy range (px and rapidity) the statistical model ceases to 
provide a good physical description of the data, and what takes its place. Recently, 
a lot of interest has been aroused by the fact that several intermediate energy 
observables {px ~ IGeV, below the scale described by pQCD) can be described in 
terms of parton coalescence (Pries et al. 2003; Molnar and Voloshin 2003). While 
this model's ability to describe the data is compelling, its failings on a fundamental 
level (energy-momentum conservation and entropy non-decrease) mean that is a 
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limiting case of a larger process. In particular, it should be possible to combine 
coalescence and statistical hadronization as two limits of a single framework (this 
is especially true for non-equilibrium statistical hadronization: After all, it can be 
argued that statistical hadronization with 7^ ^ 7^ 1 is another way of viewing parton 
coalescence into massive states, with probabilities given by thermal weights). 

It should be reiterated, however, that these advances in theoretical un- 
derstanding are very unlikely unless the right statistical hadronization model is 
experimentally selected. Max Planck found that black body radiation fits a 
distribution years before he found a fundamental description of this phenomenon. 
Similarly, one hopes our rough phenomenological fits will help elucidate one of the 
greatest unsolved mysteries in modern physics: The structure of that strange, com- 
plex and unpredictable medium which we have naively been calling empty space. 
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APPENDIX A 
Coordinates and units 



Throughout this dissertation we use the so called "natural" system of units, in 
which 

h^c^kB^l (A.l) 

Given a unit of energy, the conversion (he) will make it an inverse of the distance. 
Thus, 

he = 197MeVfm, IGeV = 5.4fm-^ (A.2) 

When dealing with a collision at ultra-relativistic energies, it is convenient to split 
the momentum into components parallel and perpendicular to the collision and 
express the parallel part as an energy-momentum rapidity 

it is easy to show that y transforms linearly under a Lorentz boost 

y^y + yo (A.4) 

it is also easy to show that 

E — niT cosh.{y) pl — mTsmh.{y) (A. 5) 
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where we defined the "transverse mass" mr — \/Vt + Hence, a particle's 4- 
momentum can be parametrized as 

/ my cosh (i/) ^ 

'Pt cos(0) 
"Pt sin(0) 
y mTsinh(y) j 

Experimentally measuring rapidity can be problematic, as it requires particle iden- 
tification. Hence, many experiments measure pseudo-rapidity 



(A.6) 



1 , fP + PL 

2 \P-PL 



-In 



tan 



(A.7) 



where 9 = tan(]9r/pL)- in the ultra-relativistic limit m << p —>■ E 

2p p2 _ 

so the approximation works much better in the central rapidity {pi «p) region. 
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APPENDIX B 
Fitting 



Suppose we have a theory which relates variables y and x through a parameter a 

y = f{x,a) (B.l) 

and a set of uncorrelatcd experimental data points {yi,Xi) with an experimental 
error cTj. The probability that the "true" value y gives the experimental result i/i is 
given by some probability density function P{y, yi, (Ti). Hence, the probability that 
N data points y were measured as y^ is given by the likelihood function 

N 

Ptot = \{P{f{x,,a),y,,ai) (B.2) 

i=l 

by maximizing Piot, or (for mathematical simplicity) 

N 

\og{Ptot) = Yl log {P {f{xi, a),y^, ai)) (B.3) 

i=l 

we can find an estimate of the parameter a. 

In particular, if the error function is a Gaussian 

P{x, Xi, ai) = -^e"^^ (B.4) 

the Maximum likelihood corresponds to the minimum where 
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By Taylor-expanding in the region of the minimum 

X^~Xmm+-^(«-«mm)^= (B.6) 

2 

and remembering that Ptot ~ the error on the parameter a can be found by 

Aa = (B.7) 

V 

This procedure can be generahzed from one parameter a to a set of parameters cci 
(In this thesis, minimizations with 13 parameters arc performed). The error then 
becomes the determinant of a matrix of correlation coefficients 

^a,, = ; ^ . (B.8) 



However, in this case the structure of the minimum in parameter space becomes 
non-trivial, with contours, plateaus, saddle points etc. making finding of the "true 
minimum" potentially very hard. In this thesis, we rely on the MINUIT numerical 
minimization package (James and Roos 1975) to minimize but we test the 
dependence on the parameters by explicitly calculating profiles (x^ as a function 
of one parameter) and contours (x^ as a function of two parameters) . 

The x^ fit is also a test of how well does the theory y — f{x, a) describe the 
experimental data. After a is found, the minimized x^ can be used to obtain the 
"confidence level" of the fit, ie the probability, given f(x, y, amin)-, that the total x^ 
obtained in an experiment is of that value or higher. 

/oo ^ poo 
^ Y[P{x,x,,a,)^ J ^ e-^'lHz (B.9) 

the quantitative dependence of CL(x^) as a function of x^ and the degrees of free- 
dom is shown in Fig. B (Hagiwara et al. 2002). As can be seen, the confidence 
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Degrees of freedom n 
Figure B.l: Confidence levels as a function of degrees of freedom. 

levels tend to an asymptotic limit as the degrees of freedom increase, corresponding 
to an expectation of — DoF. This hmit is vahd in the fits to particle spectra 
performed in chapters 3 and 4. However, if the degrees of freedom in the fit are 
limited (as in the fits to particle yields discussed in chapter 2), the confidence level 
depends strongly on the number of degrees of freedom, and a x^/DoF considerably 
less than 1 is required for the fit to be statistically significant. 

It should be remembered that the confidence level plot is only significant 
when all errors are Gaussian. This is generally not true if large systematic errors 
are present. In this case, when these errors can be estimated precisely they should 
be taken into account by summing systematic and statistical errors 

x' = y( /("y")-^; \\ (B.io) 

Z-^ \ ^statistical I systematic I V / 

1=1 ^^i / 

If, however, this can not be done statistical stops being a quantitative measure of 
the fit's validity. In chapter 3 we encounter precisely this problem with Kg. 
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APPENDIX C 
Relativistic phase space 



The relativistic phase space for a physical system is the volume in 4-momentum 
space corresponding to the system's invariant mass (which we call m. For many- 
particle systems, it corresponds to \/s) 

n ^ j d^pSiPi^p" - m^) ^ J d^p dE5 {^/ - p^ - m?^ (C.l) 



This integral can be simplified keeping in mind that 

1 

\E\ 

Idxi x=xo 

differentiating the delta-function in Eq.( C.2) by E it can be seen, immediately, that 

d'pj^ ^6 [e - V^^^) = j ^5 {e - V^^^) (C.2) 

It can be seen that this expression is relativistically invariant by Lorentz- 
transforming a volume of phase space explicitly (without loss of generality, in the z 
direction) 

dp'^dp'ydp'^ = dp:,dpy-f{dpz - vdE) = dp^dpydp^-f ^1 - 

Since 

dE d — - pz 



we get 

E^ 

'e 



dp^dpydpz ^1 - = dp^dPydPzl - ^'^j = dp^dpydp^ 
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Hence, we have proven that 

E' ~ E 

as required. 

We have therefore derived the relativistic phase space for a particle of mo- 
mentum p, energy E and mass m. For a system of N particles with momentum 
p and energy E, such as those discussed in Chapter 2, it is easy to generalize our 
result to 
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APPENDIX D 
Publications 

Several parts of this thesis have been published in peer- Reviewed Journals. We shall 
give a listing of publications contained in this thesis, together with the chapter which 
follows the publication most closely. 

• (Torrieri and Rafelski 2003a; Torrieri and Rafelski 2001a; Letessier et al. 2001) 
Chapter 3 

Paper (Torrieri and Rafelski 2001a) was chosen for the loP select article 
collection. See http://www.iop.org/Select/ 

• (Torrieri and Rafelski 2002; Rafelski et al. 2001; Torrieri and Rafelski 2001b; 
Letessier et al. 2003) Chapter 4 

• (Torrieri and Rafelski 2004) Chapter 5 

• Part of chapter 2 (the inclusion of width in particle yields) (Torrieri et al. 
2004)has recently been submitted to Communications in Physical Computing. 

• Chapter 6 (Torrieri and Rafelski 2004) is part of a continued research effort 
which will eventually be published. 
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